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ABSTRACT. A survey is given of some of the more outstanding developments in 
magnetism during the post-war years, including the great extension of our knowledge of the 
ferromagnetic domain through the theoretical work of Néel and others and by the systematic 
use of the Bitter powder figure techniques; the main ferromagnetic, paramagnetic and 
nuclear resonance phenomena and experiments on the diffraction of neutrons by anti- 
ferromagnetic crystals are also described. 


§1. INTRODUCTION 


N choosing the subject of this Presidential Address I had in mind the fact 
| that in 1939 an ‘International Réunion d’Etude’ was-held at Strasbourg, 

when practically every branch of magnetism was surveyed in a then 
up-to-date manner. In 1950 a somewhat similar international colloquium was 
held at Grenoble, and on this occasion the subjects of ferromagnetism and 
antiferromagnetism were widely discussed. I had the good fortune to be present 
at both conferences, and so I want to talk to you today about the developments 
in certain branches of ferromagnetism which have taken place in the past ten 
or twelve years, and which have more particularly interested me. ‘lhe fact that 
time will not permit me to deal with all important branches today must not be 
taken to indicate that I have not been interested in them. 


§2. THE STUDY OF DOMAINS 


I should like first of all to deal with some rather beautiful and striking 
extensions of our knowledge concerning ferromagnetic domains. We all know 
that Pierre Weiss put forward the idea of a ferromagnetic domain in 1907 in 
order to explain the fact that we may easily magnetize a specimen of fairly pure 
iron to a high degree of technical saturation by the application of a very small 
external field. In 1935 Landau and Lifshitz* published a remarkable paper 
which was almost completely ignored for ten years, and has only recently found 
its way into English papers and treatises on magnetism. ‘They asked the following 
question. Why should several domains instead of one be formed in a small 


* See note at top of References at end of paper. 
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specimen of a single crystal which has once been magnetized to satura tion” 
The answer they gave was that a single domain is extremely improbable because ~ 


of the huge magnetostatic energy which would in such a case be associated with — 
free poles on the exterior of the specimen. By a series of simple diagrams they | 


proved that a single domain must inevitably break down into many domains, 


magnetized consecutively in opposite, or as we say in antiparallel directions. 
Consequently, the positive and negative magnetizations at the ends of many 
such domains neutralize one another, and Landau and Lifshitz showed how the 
effects of these magnetizations might be still further reduced by the establishment 
of closure domains which ensure that little or no magnetic flux leaves the surface 
of a specimen. 


It might at first sight appear that any specimen would therefore contain an 


infinite number of domains, but we have to take account of the fact that 
neighbouring antiparallel domains must be separated by a wall. Now, energy 
is necessary for the formation of such a wall, a Bloch wall as it is called, because 
the wall must consist of magnetic material in which the directions of the electron 
spins in neighbouring atoms must change continuously as the wall is traversed 
from one domain to the other. We know that the energy of interaction between 


two neighbouring electrons whose spin vectors make a small angle ¢ with one — 


another is given by the expression E,S?¢ where E, is the exchange integral. 
Therefore it is easy to see that a gradual change in direction of magnetization 
throughout a long line of atoms requires less energy than that required for the 
same change in direction of magnetization over a small number of atoms. 

But the wall cannot have an infinite thickness, because the energy of magnetic 
anisotropy plays a part. It is well known that we have to supply different quantities 
of energy per unit volume to magnetize a crystal of a ferromagnetic metal in 
different directions, as one can see from a glance at the (J, H) curves for a single 
crystal of iron. Now, the majority of the electron spins contained within a Bloch 
wall must be directed away from the directions of easy magnetization, and so 
the wall must possess magnetic anisotropy energy approximately proportional 
to the thickness of the wall. We therefore see that the thickness of a wall is 
settled by the conflicting claims of the exchange energy between the electrons 
within it and their magnetic anisotropy energy. In the case of iron it is calculated 
that a wall between two antiparallel domains is about 10->cm thick, or consists 
of about 300 atoms in one line. 

But I have been talking as if it were always necessary that the domains should 
have their magnetic vectors directed antiparallel or in opposite directions on 
the two sides of a domain boundary. ‘This picture is indeed correct for the case 
of a single crystal of cobalt and also for a single crystal of the compound MnBi. 
A moment’s reflection will show that this cannot be the normal arrangement for 
single crystals of iron or of nickel. We know that in the body-centred cubic 
iron crystal there are six directions of easy magnetization, viz. the [100] 
directions, along the cube edges. In the face-centred cubic crystal of nickel at 
room tempcratures the eight easy directions of magnetization are along the cube 
diagonals, the [111] directions. Yet, the idea that there must be no free surface 
magnetization on the faces of a domain wall is all-important and leads us to new 
concepts both of the nature of the domain structure of a demagnetized iron 


crystal, for example, and of the way in which that structure changes when a 
magnetic field is applied. 
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‘Thus, the old familiar simple picture, fig. 1 (a), of a demagnetized iron crystal 
as made up of little cubes, each magnetized parallel to a cube edge and arranged 
anyhow within the crystal, must be abandoned and replaced by a picture of the 
kind shown in fig. 1(d). In the new picture it is obvious that the specimen is 
magnetically neutral and that the intensity of magnetization is continuous across 
every domain boundary, a feature which I shall be at pains to emphasize in 
many of the diagrams and photographs which will be shown today. Figure 1 
presupposes that the crystal surface is an exact crystal plane; but we shall see 
later that when the surface is inclined at a very small angle to a true crystal plane, 
the main domain boundaries shown by experiment are usually somewhat 
complicated and often very beautiful. In addition, of course, we have to take into 
account the closure domains which are found at the edges and sides of a crystal. 

Let us therefore consider how domain boundaries may be rendered visible. 
This may be done by what is known as the Bitter figure technique. We place a 
small quantity of a finely divided ferromagnetic powder, or a drop of a colloidal 
solution of magnetite, upon the surface of the crystal. In such experiments the 
particles of powder etc. simply replace the iron filings which we used in our 
schooldays. It is easy to understand that on a crystal surface where there is no 
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Fig. 1. Diagrammatic representation of the domain structure of a demagnetized iron crystal. 


free surface magnetization the line or narrow region in which a domain boundary 
intersects the surface must be manifest by a stray magnetic field. Consequently, 
the finely divided and permanently magnetized particles will be acted upon by 
such stray fields and form deposits in the manner discussed by Kittel (1949 a, b) 
and Mee (1950). 

Of course, it is imperative that the crystal surface shall be most carefully 
polished and free of surface strain, so that after preliminary polishing with 
emery down to the finest grade, it is usual to polish the surface electrolytically : 
for this purpose a useful bridge circuit has been devised (Bates and Mee 1950). 
Finally, if necessary, the sample is annealed in vacuo at 900°c. 

Perhaps a word about electrolytic polishing will not be out of place here. 
The specimen is mounted in a polystyrene mould with its edges well protected 
with polystyrene glue. It is then mounted on a wheel which can be rotated 
at a chosen rate to cause the surface to be uniformly polished. I imagine that 
most of you know that electrolytic polishing is caused by a film of high resistance 
which covers the anode. Above any tiny surface projections the film will be less 
thick than above flat regions of the surface, and therefore the protuberances are 
more vigorously attacked. As it is essential that the viscosity of the solution 
used in the polishing bath should be everywhere the same so that ionic mobilities 
may be constant, we now rotate both anode and cathode on a common spindle 
(Bates 1952), and we use the special bridge circuit already mentioned to maintain 
the polishing current at the optimum value. 

The well-prepared crystal surface is then mounted horizontally in a horizontal 
or a vertical magnetic field and observed through a metallurgical microscope. 

38-2 
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A drop of the colloidal solution is placed on the surface and a microscope cover 
glass placed on it, when the Bitter figures quickly form. If one uses a vertical 
field one measures its intensity by means of a flat search coil, but in the case of 
a horizontal field one uses a small magnetic potentiometer. 

I shall now review, very briefly, a few of the patterns which one obtains by ~ 
this technique. If one obtains a maze pattern, it is generally because the surface 
is not strain-free, and it must be repolished and annealed. It is easy to prove, 
however, even by the use of imperfectly prepared specimens that the main or 
principal domains are definitely situated below the crystal surface. 

A very distinctive feature of the Bitter figures on a well-polished surface of an 
iron single crystal which is almost parallel to a true crystal plane is the appearance 
for the demagnetized state of a ‘fir-tree’ pattern. Such ‘fir-tree’ patterns have 
been examined and explained by Bozorth, Williams and Shockley. Their 
explanation can readily be found in at least two modern treatises on magnetism 
(Bates 1951, Bozorth 1951), but it may help you if I place their interesting 
diagrams on the screen. If the crystal surface were ideally a true crystal plane 
we should get few, if any, Bitter figures upon it, but the further the surface 
departs from the true crystal plane the heavier become the Bitter figure deposits. 
‘There is no doubt about the essential correctness of the beautiful explanation 
given by Bozorth and his collaborators. 

Let me point out to you that the trunks of the ‘fir-trees’ are formed by the 
Bloch walls which separate alternate domains which are magnetized respectively 
parallel and antiparallel on either side of a trunk, and that the branches are 
formed by closure domains, magnetized in a direction roughly perpendicular 
to the trunks, whose function is to take care of a great quantity of magnetic flux 
which would otherwise emerge from the crystal surface. Indeed, from the 
changes in appearance of the ‘fir-trees’ from place to place over the crystal 
surface we can now deduce exactly how that surface is shaped, and it almost 
seems as if Bitter figures, or rather absence thereof, is the best evidence of a truly 
ground and annealed iron (100) surface that we can possibly supply, for the 
branches always point downhill. 

There are at least three separate methods by which it can be conclusively 
ascertained how the principal and closure domains are magnetized. We may 
scratch the surface with a fine glass fibre and note the deposits in the stray field 
at the scratch, we may use a tiny probe of permanently magnetized metal, or we 
may note the way in which deposits form on the sides of any little depressions in 
the surface because of the free poles on the sides of the depressions. 

I want later to deal with some quantitative aspects of powder figures, but 
first I should like to consider the application of the newer domain concepts to 
one or two well-known problems of the magnetic properties of iron single crystals. 
Let us follow Néel’s description of magnetization processes. He considers that 
in the demagnetized iron crystal each domain vector is aligned parallel to one of 
the six directions represented by three cube edges and their extensions in the 
negative direction, and he designates this condition a state of six-phase 
magnetization, or Mode I magnetization. 

If we now apply a small external field to the crystal, supposing the latter for 
clarity to be of ellipsoidal shape, then before it can produce an effective field 
inside the crystal three magnetization phases must disappear, as they are entirely 
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contrary to the direction of the external field. When they have completely 
disappeared we have left three-phase or Mode II magnetization. It is now possible 
by further increasing the applied field to produce an effective field within the 
crystal itself and thus cause one of the phases to be so much less favoured than the 
others that it disappears, when we are left with two-phase or Mode III 
magnetization. Eventually, on further increasing the field only single-phase 
or Mode IV magnetization can persist. 

If the single crystal of iron is shaped as an oblate spheroid whose equatorial 
plane is a (100) crystal plane and we apply an external field H, in this plane at 
any angle @ to the nearest [100] direction, we can deduce by a commonsense 
argument that the resultant of H, and the demagnetization field of the spheroid 
must lie symmetrically between the two nearest [100] directions, i.e. the 
resultant field must lie along the [110] direction in the equatorial plane. But 
what must then be the shape of the domains? They must be leaf-shaped (French, 
feuillet) with their boundaries perpendicular to the [110] direction. 

I need not deal further with this problem; it has been considered in extenso 
by Néel and also by Lawton and Stewart. I will merely place on the screen 
one or two slides which illustrate how completely the components of magnetization 
of the ellipsoidal specimen parallel and normal to the applied field should 
theoretically behave and the excellent accord with experimental results obtained 
many years before the theory I have just outlined was given. But I should like 
to refer to what is known as the Kaya relation, which obtains for a long rod or 
ellipsoidal iron single crystal whose axes of easy magnetization make angles of 
direction cosines «,, %) and «, with the long axis of the ellipsoid. Obviously, the 
component of magnetization normal to an applied field will now always be very 
small, and a commonsense argument shows that when a very small field eventually 
acts inside the specimen, Mode II magnetization must reign and each of the 
domain vectors be directed parallel to one of the three directions of easy 
magnetization. Consequently, we must have the parallel component of 
magnetization equal to I,/(a,+a,+ 3), a relation which Kaya deduced for the 
retentivity of such a specimen. 

The time has now come for me to deal with quantitative experiments on 
powder figures, and I must first mention an exceedingly fundamental experiment 
made by Williams and Shockley who cut a single crystal of silicon—iron in the 
shape of a picture frame. ‘They were able to cause a 180° boundary to form 
approximately along the middle of the sides of the picture frame, and to cause 
it to move by applying a magnetizing field by means of a coil wound around 
one side of the frame. ‘They thus made some very striking quantitative 
measurements of the sideways displacement of the boundary as a function of the 
magnetization. 

While I am dealing with picture frame measurements, I should like to 
mention those of Stewart (1951) who examined the velocity V of the displacement 
of the 180° wall and found that it obeyed the relation, first given by Sixtus and 
Tonks, V = A(H — H,); here H is the applied field and H, is a field analogous to, 
if not actually equal to, the coercivity of the material. But, in these experiments 
A was found over 1000 times smaller than the values usually found with 
polycrystalline materials. ‘The suggested explanation is the enormous braking 
effect of eddy currents in the single crystal case. In the course of this work 
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it was found that the downward sweep of the hysteresis curve in the 
demagnetization quadrant was re-entrant, and the two phenomena are possibly 
linked together. 

The effect of eddy currents was the explanation advanced by Snoek to account 
for the fact that the virgin magnetization curve often lies partly outside the 
closed hysteresis loop when measurements are made on rod or strip specimens 
on open magnetic circuit. It is supposed that eddy currents prevent the 
magnetization from rising rapidly to its proper value when the field is changed, 
so that the true field acting on the outside of the specimen is not H—DI but 
H-—DI', where I>I', whereupon the effective field as it were ‘overshoots’, 
and the outer portions of the specimen assume an intensity of magnetization 


(@) No field applied 


(b) Field applied 


Fig. 2. Domain structure proposed by Néel. 


greater than they ought, particularly on the steeply ascending portions of the 
loop. Incidentally, Kayser has shown that even if we use an anchor ring specimen 
with no free poles we may get the same Snoek phenomenon if the material of 
the ring possesses strong magnetic anisotropy. 

In Nottingham we have made a series of quantitative experiments to test 
the correctness of certain magnetization processes put forward by Néel. 
Although the most complete results were only very recently published in our 
Proceedings (Bates and Mee 1952a,b) I should like to show you some of 
them today together with many pictures which could not be published because 
of space (and cost). In fig. 2 we have domain structures which Néel showed 
to be appropriate to a long single crystal strip of iron, with its main upper and 
lower surfaces very nearly parallel to true (100) crystal planes, and the 
[011] direction running down the middle as depicted on the upper surface of 
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the strip. When a very weak, indeed practically zero, field acts upon the strip 
parallel to this [011] direction, the principal domains form slices magnetized as 
shown in fig. 2(a) with their boundaries perpendicular to the [011] direction. 
When the field is increased we get marked changes in the closure domains as 
in fig. 2(4). Now, Néel has given a formula which allows us to calculate the 
period d of the domains, which is equal to the total thickness of two successive 
domains, in terms of the true field acting on the specimen. _ His calculation of 
_ the period is based on the particular arrangement of the principal domains and 
their closure domains as given in fig. 2(b). I might mention that we are here 
dealing with a transition from Mode II to Mode III magnetization. 

We have made many measurements of the period d both on the (100) surface 
and on the (110) planes formed by the sides of the specimen. I now place some of 
them on the screen, together with graphs showing the variation of d with H, 
and you will see that there is no doubt whatever about the essential correctness 
of Néel’s ideas, although more profound study of the experimental results shows 
that these ideas need to be modified in some particulars. For example, we feel 
fairly certain that we do get some free pole distributions on the (110) faces of our 
specimens. Again, recent work with single crystals of very pure iron cut in the 
_ manner described seems to show that the boundaries on the (110) faces are not 

vertical as drawn, but are definitely inclined at an angle of about 15° to the 
vertical. In any case, it is a very strange thought that the period of a principal 
domain structure, which is in effect the whole state of magnetization of a specimen, 
can be so markedly influenced by what occurs in a set of closure domains on the 
sides of the specimen! This is a feature to which I shall return later. 

Perhaps this is a convenient point at which the effects of impurities or 
inclusions in the crystal may be introduced. We all instinctively feel that if we 
had a perfect crystal it would be possible to magnetize it to technical saturation 
in an extremely small field, and the obvious mode of increase of magnetization 
in a particular direction is by boundary displacements. So we at once seize 
upon the idea that the coercive field is a limiting field which must be overcome 
before the boundary will move from one stable position to another. It also seems 
reasonable to suppose that the initial permeability p, is a measure of the internal 
restoring forces which act upon a domain boundary subjected to small reversible 
displacements only, the coercive force then being the maximum restoring force 
which a boundary can stand. 

Now, it is well known that Kersten put forward an explanation of coercivity 
in terms of the movement of a domain boundary from one set of inclusions to 
another. But Néel showed that this theory was incorrect because no account 
was taken of the huge quantities of magnetostatic energy which would be 
associated with a wall passing across or through an inclusion; indeed, he showed 
that a factor of about 100 was omitted, and he thinks that the real effect of 
inclusions is to produce large fluctuations in the intensity of magnetization or 
its direction inside the domains themselves. 

Let us look for a few moments at some pictures of the Bitter figures which 
were recorded in the neighbourhood of inclusions (Bates and Mee 1952b). 
In many of these pictures you see that a complicated closure domain structure 
exists round an inclusion. In fig. 3 are given some of the types predicted by 
Néel, and these are fully confirmed in many of the photographs which I now 
project on the screen. 
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I should like to complete this somewhat minor incursion into the realms of 
coercivity by reference to an important piece of work by Stoner and Wohlfarth. 
It seems obvious to me that for some time three sets of people, Koenigsberger 
in Germany, Néel in France, and Stoner and Wohlfarth in Leeds, were grappling 


with the same question. ‘This question may be put in the following way. 


If we have a particle of a ferromagnetic element so small that it never contains 
more than one single domain, how does it behave in a magnetic field; and what 
field would be required to reverse its direction of magnetization? 

The most complete picture of what may happen was given by Stoner and 
Wohlfarth and I am sorry that I can only give a brief outline of it here. Neglecting 


™~ 


any energy of magnetic anisotropy or of internal strain, the total energy of an — 


ellipsoid of revolution, with longitudinal and transverse demagnetization 
coefficients D;, and Dg respectively, is given by $(D,—Dg,)1,2V«?, where V is 
the volume of the ellipsoid and « is the angle which J, makes with the long axis. 
When fields are applied parallel to the long axis, we obtain a rectangular 
hysteresis loop with the coercivity H,=(D,—Ds)[, A graphical method of 


Fig. 3. Closure domains around an inclusion (after Néel). 


calculation shows that for an assembly of such ellipsoids which do not interact 
with one another the average coercivity is about 0-96 of this value; and, if the 
ratio of the long and short axes is 1+¢«, where « is small, then (D, — Dg) = 87/5, 
so that H,=4 100 in the case of iron. 

Let us apply these ideas to the case of a magnet made of compressed finely 
powdered iron, where we must take account of mutual interactions between 
domains. Assuming that a number of cavities of relative volume wv are uniformly 
distributed in 1cm* of uniformly magnetized material, then the energy of the 
whole is $7(1 —v)v1,?, which is 1—v times the energy which would be possessed 
by the material without cavities. H, is therefore correspondingly decreased. 
Taking the case of simple cylindrical particles for which the difference in the 
two demagnetization coefficients is 27, and assuming that these particles of 
density d) are compressed together to form a structure of density d, we find 
H,=0:26n(1—d/d))I,, giving a value of about 500 oersteds for iron.’ Now, 
Weil investigated powder magnets made of 70% iron and 30% cobalt and proved 
that H, 1s proportional to 1—v, a result in complete agreement with these views. 

We have just seen how important a cavity or an inclusion in the surface of a 
metal can be; it is also of interest to discuss the effects of what might be lightly 
termed an exclusion. In fig. 4(b) is shown a part of the (100) surface of a single 
crystal of pure iron. The lowest corner of the crystal has been cut off along a 
(111) plane as shown. We know that if the crystal had not been thus maltreated, 
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the main and closure domains for the demagnetized state would follow the 
pattern of fig. 4(c). What was actually photographed by Mr. D. H. Martin 
at Nottingham is reproduced in fig. 4(a). A partial interpretation of this 
somewhat beautiful pattern is given in fig. 4(b) where the dotted square 
indicates the extent of the area photographed in fig. 4(a); the word partial is 
used to suggest that the closure domain system which must exist at the lowest 
edge of the (100) surface shown in fig. 4 (5) has not yet been photographed or 
examined. ‘The tree patterns found in the photograph do not particularly 
interest us at this stage; they are only the outward sign that the surface is not 
ideally cut parallel to the (100) plane. But what does interest us is the injection 
of a new domain system cutting across the simple parallel set of boundaries which 
would exist if the corner of the crystal had not been removed. ‘The somewhat 
startling idea, to which I referred earlier, that according to Neéel’s theory the 
magnetization of a strip is determined by what happens at its edges, is now 
supported by this direct experiment. 


(a) (6) 
Fig. 4. The effect of an exclusion from a surface on domain structure. 


I said earlier on that we had to abandon the old cubic block picture of domains 
shown in fig. 1 (a). But this statement may not be entirely correct, for I now 
show you some pictures, taken in my laboratory a few days ago by Mr. D. H. 
Martin, which indicate that two neighbouring domains in a silicon-iron crystal 
can have their vectors directed head-on to one another, and you see that the 
boundary between them is a very tortuous one. 

To conclude this portion of my address, I should like to show you a film of the 
formation of domains on the surface of a single crystal of cobalt which is particularly simple 
from a magnetic point of view because it possesses only two directions of easy magnetization 
parallel to the hexagonal axis. As a prelude to the film we will first examine a series of slides 
of many of the features which will appear in the film (Bates 1951). 


§3. GYROMAGNETIC MEASUREMENTS 


In 1922 there appeared the first of a series of researches on the measurement 
of the gyromagnetic ratio for ferromagnetic metals which culminated in the 
brilliant measurement of the gyromagnetic ratios for paramagnetic ions of the 
iron and rare earth groups by Sucksmith; I should like to spend a little time 
in dealing with the present state of our knowledge of the value of the ratio for 
the ferromagnetic elements. You all know that when such an element is 
magnetized there is a change in the measured intensity of magnetization which 
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is accompanied by a corresponding change in angular momentum, and — 
experiments have been designed either to measure the change in angular — 
momentum which accompanies magnetization or the change in magnetic | 
moment which accompanies a change in rotation of a ferromagnetic specimen. 

The results obtained in England and America prior to 1939 are summarized — 
in the following tables. 


Table 1. Collected Values of pe/m for Iron 


Authors Ratio Material Method 
Chattock and Bates 1:00+1% Tron Direct rotation by 
magnetization 
Sucksmith and Bates 1:0041% Iron Modified resonance 
Coeterier 1-:00+ (1%?) Iron Controlled resonance 
Ray-Chaudhuri 1-004 1°/ Iron oxides Resonance 
Barnett 1-03841% Iron Magnetization by 
rotation 
Barnett 1:03240-3% Iron Modified resonance 


Table 2. Collected Values of pe/m for Ferromagnetics other than Iron 


Authors Ratio Material Method 

Sucksmith and Bates 1-002 + 0:004 Nickel Modified resonance 

Barnett 1:046+ 0-004 Modified resonance 
Sucksmith and Bates 1-002 + 0-006 Heusler Alloy Modified resonance 
Sucksmith 1-030 + 0:064 Cobalt Modified resonance 

Barnett 1-078 + 0-007 Modified resonance : 
Sucksmith 0-990 + 0-048 Magnetite Modified resonance 
Coeterier S351 7) Pyrrhotite Controlled resonance 


In these tables the results are expressed in terms of the gyromagnetic ratio ; 
it is defined by p=(change in angular momentum)/change in magnetic 
moment) =U/u. If the magnetic carrier is an electron without spin describing 
a plane orbit we should have p=2m_/e, if the carrier is a spinning electron we 
should have p = m/e, while for the general case in which the carrier is an electron 
system for which the Landé splitting factor is g, we should have p=2m/ge. 
From the tables it follows that g is very nearly equal to 2, so that we deduce that 
the elementary magnetic particle owes its magnetic properties almost entirely, 
if not entirely, to electron spin alone. 

But, it is clear from the tables that some differences exist between the several 
determinations with any given element, and these are important, because if the 
higher values of p are correct then g must be slightly less than 2, from which 
it follows that orbital motion of the electron must play a small but definite part. 
Naturally, one would like this question to be cleared up completely. 
Unfortunately, we have had only two determinations of the ratio since the last war. 

The first-was by Meyer, of Strasbourg, who used the controlled resonance 
method and made measurements with iron and nickel and with the alloys FeNi, 
(FeNi)Co, FeCo and Fe,Ni. He paid careful attention to the width of the 
slit in front of the photoelectric cell whose activation caused a magnetizing 
field to be reversed. The average value of g which he obtained for the metals 
and alloys was 2:001+0-005. The work appears to me to have been done with 
extreme care. 
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Thus far so good; but we have a second method by Scott (1951) based on the 
method, and using much of the apparatus, which was so successfully used in a 
determination of the value of the specific charge of the electron by Kettering 
and Scott in 1944. The method was a kind of controlled resonance method, in 
which changes in amplitude were caused by periodic, synchronized reversals 
of the magnetizing current in a coil wound directly on the specimen; these 
changes were measured, together with the corresponding changes in the 
magnetic moment of the specimen. The final result for the gyromagnetic ratio 
for iron was (1-:0278 + 0-0014)m/e or g=1-946. Again, the work appears to me 
to have been very thoroughly carried out. 

So we still have differences in the values for g; and perhaps the truth of the 
matter is that while we can make relative measurements of changes in 
magnetization and in angular momentum with great accuracy, we cannot yet 
make absolute measurements with the necessary degree of certainty. 


§4. FERROMAGNETIC RESONANCE 

Let us, however, turn now to some experiments on a new phenomenon 
discovered by J. H. E. Griffiths in 1946, for which he was fittingly awarded the 
Vernon Boys Prize by the Physical Society in 1951. Let us suppose that a 
demagnetized ferromagnetic metal is exposed to a field H, along the z axis and 
that this field is slowly increased from zero; then we know that the induction 
of the specimen must change according to the usual (B,H) curve. Now, let 
us further suppose that at any chosen value of H, we superimpose upon it an 
oscillating field of very small amplitude. When the two fields are applied along 
the z axis we may thus find what is usually termed the incremental permeability. 

But when the oscillatory field is, say, along the x axis we may then measure 
the transverse incremental permeability. In either case for low frequencies the 
incremental permeability falls from a value greater than unity for low values 
of the steady field to unity when the latter is high enough to cause saturation. 

When, however, the frequency of the oscillating field is very high, e.g. in the 

microwave region, a resonance effect is observed and the incremental 
permeability, instead of decreasing to unity as would be expected, actually 
rises from a low value for a small steady field to pass through a maximum and 
then decrease towards unity as the steady field H, is increased to a very high 
value. In other words, when yp, is plotted against H, we get a kind of broad 
resonance curve, as shown in the slide (Standley 1949). 

We may explain this resonance effect in the following way. ‘The 
ferromagnetic carriers may be considered as under the influence of an effective 
field H’, so that they perform Larmor precessions of frequency 

Pmeveyenynic)in Pe  Seeease (1) 

When the Larmor frequency is equal to that of the oscillating field, resonance 
is set up and energy absorbed at a maximum rate from the oscillating field. 
Alternatively, we may suppose that the absorption may occur between a series 
of Zeeman levels in the ferromagnetic material, so that hy=gu,H’, where 
fy is the Bohr magneton and g is the Landé splitting factor. Writing w=2z1, 
we have w=g(e/2mc)H'=yH’' where y is defined as the magneto-mechanical 
ratio. It isa pity that this term has extensively crept into the literature in place 
of the term gyromagnetic ratio coined by O. W. Richardson in 1908 and used 
extensively in England for over a quarter of a century. 
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Using a semi-classical model Kittel showed that the general resonance 
equation for an ellipsoidal specimen with axes x, y, s might be written as 


w= (H,+(N,—N,)M,\A,+(Ne—Ne)Ma], «se (2) 


where M, is the saturation magnetization of the material along the direction | 
of H, and N,, N,, and N, are the demagnetization coefficients along the three 
axes. ‘Two cases are frequently met in practice: firstly, the case of a sphere, 
where we find w=yH,, and secondly a very thin sheet parallel to the xz plane, 
where we have w =>(B,H,)"?, a statement derived by both Kittel and Luttinger 
(1948) and Van Vleck (1950) by quantum theory treatment of the problem. 

We therefore see that from such ferromagnetic resonance measurements we 
may obtain y and hence calculate g. Let us for a moment consider the usual 
experimental arrangements. The specimen usually forms part or the whole 
of one wall of a cavity tuned to resonate at a fixed microwave frequency v, which 
is often about 24000 Mc/s, and the Q or circuit magnification of the cavity is. 
measured for a series of values of H,. It is reasonable to assume that when H, 
is changed only the losses in the ferromagnetic specimen itself are affected, and 
so we may calculate from the Q measurements the quantity (uwpp,)"?, where 
Ltn is the permeability of the specimen determined by resistive measure- 
ments of this kind and p, is its electrical resistivity. If however we retune the ~ 
cavity as H, is varied, we thereby measure pz, the inductive or dispersive 
permeability. 

Now, at high frequencies we must obviously treat the permeability as a 
complex quantity. Writing w.=,—ju, we can show that wp = (py? + Me”)? + pg, 
bey, = (My? + Me”)? — a. Hence, from simultaneous measurements of ppg and py, 
the two quantities 4, and , may be found. 

Such experiments have so far been made on about a dozen different 
ferromagnetic materials, and it is significant that the values of g so obtained 
generally lie between 2:0 and 2-3, although in Nottingham a value greater 
than 3 has been found for manganese arsenide. Unfortunately, in my opinion, 
it is customary among workers in this field to refer to their values of what is in 
effect the magneto-mechanical ratio factor as g values, and to describe the 
corresponding values obtained from the older direct measurements of the 
gyromagnetic ratio as g’ values. (It seems obvious to me that gp (Richardson) 
for gyromagnetic data and gg (Griffiths) for ferromagnetic resonance data would 
have been more appropriate, but the unfortunate terminology now seems rather 
firmly established.) 

A few values are given in the following table, using the objectionable symbols 
under protest: 


Table 3 
Specimens g go 
| (alaaapele: 
Ni BS) 1°92 2-00 
Co Di, S37) 1-94 
Fe 2:°12—2:17 Ways 2°00 (1:95)* 
Magnetite 2:20 1-93 2:02 
Permalloy 2:07—2:14 1:91 — 


(1) Measurements by Barnett. 
(2) Measurements by other workers. * Scott’s value. 
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It has been shown by Kittel (1949 b) and by Van Vleck (1950) that although 
g and g’ are both defined by ge/2mc=AM/J, where AM is the change in 
magnetic moment and AJ is the corresponding change in angular momentum, 
it may be suggested that they do in fact refer to different quantities. For, when 
we determine the gyromagnetic ratio, because angular momentum is conserved 
the change AJ/,,;,; for the whole system is zero. But in the ferromagnetic resonance 
experiments AJj,14;= A/27. 

Again it has been shown that 


gie[2me = Mevint Monit” while ge/2me = Maint Mont | 
J spin t+ orbit spin 
so that when the orbital contribution is much smaller than that due to spin, 
we have De to) gi eee a Pee ae (3) 


The last equation is only approximate and is based on a simple model, and 
while the results in the above table may appear to indicate that g’<2 and g>2, 
it is clear that the equality is not satisfied. It may be pertinent to add that eqn. (2) 
for w* holds for zero or small damping. 

The experimental results for certain manganese ferrites, manganese nitride 
(Mn,N) and Heusler alloys show that g is equal to 2-00 within the limits of 
experimental error of about 2°. In other cases, g is greater than 2-0 and in the 
case of MnAs, mentioned above, is greater than 3. The simple relation between 
g and g’ cannot account for such cases. Moreover, in the cases of nickel and 
permalloy no changes in the g values are observed when the specimens are raised 
from room temperature to above the Curie point; changes have however been 
recorded by some Japanese workers (Okamura, Torizuka and Kojima 1951) in 
the case of cobalt ferrite as the temperature approaches the Curie point, and 
Standley and Adam (unpublished) have found a similar behaviour with MnSb. 

In order to explain the broad resonance bands which are observed Stevens 
has suggested that electrons may fluctuate from one lattice site to another, 
e.g. from d?® to d? sites in the case of nickel; if the lifetime in a given site is short 
compared with the Larmor precession a very broad resonance would be expected. 


§5. PARAMAGNETIC RESONANCE 


I have said so much about ferromagnetic resonance that it may have escaped 
notice that it is in reality only a special case of electron resonance, and that, in 
general, we ought to expect similar resonances in non-ferromagnetic materials. 
However, we do not find such resonances in the majority of substances, for 
instance, we do not find them with pure water. ‘The reason for this failure of 
the electrons to resonate in such cases is that they behave as if they are coupled 
together in pairs with the two spins directed antiparallel to one another. 
Consequently, each coupled pair forms a non-magnetic system on which 
the constant and oscillating fields have no effect. 

But in the case of the strongly paramagnetic salts we have uncoupled electrons 
available and electron resonance according to eqn. (1) must take place. 
Absorption lines are found whenever the separation of the energy levels of a 
paramagnetic ion is equal to the quantum of energy corresponding to the radio 
frequency. ‘This experimental procedure is very important, for it enables us to 
make a detailed examination of these energy levels without the necessity of making 
corrections for diamagnetism, to measure the relevant g values or splitting 
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factors directly, to measure directly small splittings of the order of lcm? — 


produced by crystalline fields and: also to measure very short spin-lattice 
relaxation times of the order of 10~!°sec. . ; 
The paramagnetic absorption spectrum of a salt may be very complicated — 


and difficult to interpret, for unless the ion is perfectly free of residual spin-orbit __ 


coupling, the paramagnetic crystal may for example be magnetically anisotropic — 
and the spin levels in a magnetic field will then depend on the orientation of the 
magnetic field to the axes of the crystalline field. Again, it is possible for the 
unit cell of the crystal to contain two or more paramagnetic ions, each exposed 
to a differently orientated though otherwise identical crystalline field, so that the 
paramagnetic resonance spectrum of each ion will be different for an arbitrary 
setting of the crystal. 
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Fig. 5. Orientation of magnetic field H, electron spin g, and nuclear spin Kg when anisotropy is 
present. The figure is drawn for g) >g), and B>A. The dotted vector Kg refers to the case 
where the coefficient giving the interaction between the nuclear magnetic moment and the 
magnetic field of the unfilled electron shell is negative. 


I cannot here go into the wealth of information obtained for example by 
Bagguley, Bleaney, Griffiths, Penrose and Plumpton, in detailed experiments 
on over one hundred salts of the iron group, which were published in our 
Proceedings, but I should like to refer to the brilliant discovery of the hyperfine 
structure of paramagnetic resonances by R. P. Penrose, whose untimely death 
we all regret. Ina flight of genius he realized that with a combination of low 
temperature and magnetically dilute ions, great resolution of the absorption 
lines could be expected. He found that at low temperatures in the microwave 
absorption band of a dilute Tutton salt, in which 95% of the magnesium had 
been replaced by copper, with A=3-2cm, with H, about 3000 oersteds along 
the K, axis, four resonance lines were observed with field separations of about 
93 oersteds, as found by adjustment of H,. ‘These lines corresponded to the 
four possible orientations in the field H, of the copper nucleus, for which 
I=3/2 units of spin. The overall structure of the lines varied with the 
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orientation of H, in the crystallographic ac plane, being a minimum along 
the K, axis, for which the electronic g factor is a maximum, and a minimum 
near the K, axis for which the electronic g factor is a minimum. 

Incidentally, a kind of physical explanation of these different g factors has 
recently been given by Bleaney (1951) and by Garrett (1950). Referring to 
fig. 5, taken from Bleaney’s paper, we may see that a magnetic field H applied 
at an angle @ to an axis of crystalline symmetry has components H cos 6@ and 
H sin @ respectively along the normal to this axis as shown. The components 
of the magnetic moment must accordingly be proportional to g,H cos@ and 
&,H sin @, where g is proportional to the length of the vector, making an angle ¢ 
with the z axis, which is formed by compounding the two components. The 
vector g is parallel to the resultant magnetic moment of the electrons, but it 
is not parallel to H. Consequently, the factors g,/g and g,/g enter into the 
equations for hy as multiplying factors of cos @ and sin @ respectively. 


§6. NUCLEAR RESONANCE 


To those of our Society who are concerned with nuclear physics, perhaps 
the most interesting use of magnetic resonance techniques is in the determination 
of nuclear magnetic moments. For this purpose we use eqn. (1) in the form 
v=2uH/h where » now stands for the magnetic moment of the nucleus in 
which we are interested. As a special case we may consider the measurement 
of the magnetic moment of the proton itself, when we have v=2u,H/h. Let 
us in this case picture a sample containing protons, e.g. a tube of water, placed 
upright in a strong, constant, horizontal magnetic field H,, whose value may 
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Generator 
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Fig. 6. Block schematic diagram for determination of magnetic moment of a proton. 


conveniently be fixed at 10000 oersteds. Around this tube let us wind a coil 
to be supplied with alternating current of a radio frequency of the order of 
tens of megacycles per second. Measurements of the impedance of this coil 
show us when absorption corresponding to proton resonance occurs, for this 
is manifested by an extremely sharp peak on the curve of absorption plotted 
against frequency. ‘This experiment is in constant use in our teaching 
laboratories at Nottingham. ‘The apparatus as set up by Dr. R. Street is shown 
diagrammatically in fig. 6. 

The coils MM are supplied with 50c/s a.c. to modulate the field in the 
gap between the pole-pieces NS of the electromagnet. The specimen tube P 
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contains a weak solution of cobalt nitrate. AB is a length of cable which 
produces a phase shift of 180°. When the circuits T, and T, are very nearly _ 
identical the signal derived from the point C is almost zero because the ~ 
introduction of a 180° phase shift by means of AB results in feeding T, and T, 


in antiphase. The circuit thus behaves as a kind of bridge with out-of-balance 


voltage derived from C. When the modulated field in the gap between the 
pole-pieces NS passes through the value required to make the precessional 
frequency of protons equal to the frequency of the h.f. currents in the coil 
surrounding the specimen, both the inductance and the associated resistance 
of the coil change. Therefore at a critical value of field which occurs in general 
at two epochs of the modulating field the bridge exhibits a slight out-of-balance 
voltage at C which is amplified by the communication receiver and presented — 
visually as a cathode-ray tube display. It is primarily the ferromagnetic 
uniformity of pole-piece material which determines the homogeneity of field 
acting in the specimen, and hence the sharpness of the resonance signals. 

The apparatus is assembled easily from components generally available in 
the laboratory, and is convenient for use as an experiment for advanced 
students and also as a method of producing very accurately known magnetic 
field intensities, for the calibration of search coils etc. For a signal generator 
frequency of 18 Mc/s the field in the gap required to produce proton resonance ~ 
is 4000 oersteds approximately. 

The importance of proton resonance is that it provides an excellent meetliee 
of measuring a magnetic field. For, if we use a fixed field of 10000 oersteds in 
this experiment, then the sharp proton resonance occurs at a frequency 
v=42-6 Mc/s. When one realizes that in some of this work resonance peaks 
less than 0-2 oersted wide are obtained when the frequency is fixed and the field 
varied,'it is clear that the accuracy of field measurement is now very high indeed. 
An idea of the order of accuracy in magnetic resonance measurements may be 
gained from the fact that the ratio of the nuclear magnetic moment of the 
proton to that of the deuteron has been given by Smaller (1951) as 
3-2571999 + 0-0000012. 

Suppose now that we have a proton situated in a crystal lattice. Then it 
will be exposed to the magnetic fields of the other protons surrounding it. 
Consequently, a broadening or shift of the resonance peak should be observed, 
and should enable the distance apart of the protons to be found. Since protons 
cannot be located by x-ray diffraction experiments, these magnetic resonance 
measurements may assume considerable importance, for together with neutron 
diffraction experiments they form the only means by which we may determine 
proton positions in a crystal lattice. 

Now, we are used in physics to being able to repeat our experiments 
ad libitum. In the case of these resonance peaks, however, the height of the peak 
may vary with time. This is because every time a radio-frequency quantum is 
absorbed a proton is turned from the parallel setting into the antiparallel setting, 
and it takes time to come back. ‘This effect has now been found experimentally, 
and an application of Boltzmann’s law shows us how small is the relative 
difference in the number of protons directed parallel and antiparallel respectively 
to the applied field. For example (Darrow 1951) if we have 2000000 protons 
in thermal equilibrium at room temperature, the whole effect is occasioned by 
1000007 set parallel and 999993 set antiparallel, a difference of only 14. 
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$7. NEUTRON DIFFRACTION EXPERIMENTS 


The mention of neutron diffraction reminds me that I have not dealt today 
with another outstanding feature of magnetic research during the last ten 
years, namely the intensive study of antiferromagnetism and ferrimagnetism. 
I should normally have dealt in some detail with these phenomena, but I have 
good reason to believe that this year’s Holweck Prizeman, Professor L. Néel, 
will devote a large part of his Holweck Lecture to their consideration. 

If we now assume that between the metal ions and the oxygen ions in 
substances like CoO, MnO or NiO there are forces, of what is called a 
superexchange type, which cause the metal ions represented by the red spheres 
in the model before you to have their spin vectors directed towards you, while 
the ions represented by the blue spheres have their spin vectors directed away 
from you, then it follows that the whole system is almost magnetically neutral. 
We should, of course, expect the system to show a complicated variation of 
susceptibility with rise in temperature, and an anomaly in its specific heat 
analogous to that found with a ferromagnetic metal, and these are characteristic 
of the typical antiferromagnetic substances which I have just mentioned. 

Frem the point of view of x-ray diffraction measurements the model represents 
two unit cells because the scattering of x-rays does not depend markedly upon 
the direction of electron spin. But, suppose we have a radiation whose scattering 
does depend upon the directions of the magnetic spin vectors in the individual 
ions. ‘Then the two kinds of metal ions represented by the two kinds of circles 
are not equivalent as far as this radiation is concerned, and the unit cell of metal 
ions is doubled linearly for this radiation. Now a beam of neutrons may be 
regarded as radiation of this type, and for neutron scattering the model represents 
only one quarter of a unit cell. Consequently, in addition to the lines which we 
get with x-ray scattering, we get other lines when we employ neutron scattering. 
Alternatively, if we carry out diffraction experiments on the material at a 
temperature sufficiently high for the antiferromagnetism to be absent we get 
similarity of the neutron and x-ray diffraction patterns, but if we make the 
experiments at temperatures sufficiently low for the antiferromagnetism to be 
present we get double the number of lines in the neutron diffraction pattern. 
Conclusive experiments on these lines have been carried out by Shull and his 
collaborators (1951), and I am sure that we shall hear very much more about them 
in the near future (cf. Street 1951). 


§8. CONCLUSION 


L have tried in the short time at my disposal to set before you some of the ~ 
more important developments in magnetism during the post-war years. I realize 
that I have left untouched many important pieces of work and I can only hope 
that what you have heard today may entice you to seek out descriptions of these 


elsewhere. 
REFERENCES 


Many references are made in this Address to works which are described in Moderin 
Magnetism by L. F. Bates (1951). Such references are therefore not given in the following 
list. 

Bates, L. F., 1951 a, #7. Phys. Radium, 12, 323; 1951 b, Modern Magnetism, 3rd edn. 
(Cambridge : University Press); 1952, Rev. Metall., in the press. 
Bates, L. F., and MEE, C. D., 1950, 7. Sci. Instrum., 27,317; 1952 a, Proc. Phys. Soc. A, 65, 
129; 1952 b, Ibid., 65, 140. 
PROC. PHYS. SOC. LXV, 8—A 39 


594 D. M. S. Bagguley and }. H. E. Griffiths 


Beanery, B., 1951, Phil. Mag., 42, 441. 

BozortH, R. M., 1951, Ferromagnetism (London : Macmillan). 

Darrow, K. K., 1951, Elect. Eng., 70, 401. 

Garrett, C. G. B., 1950, Proc. Roy. Soc. A, 203, 375, 392 ; 1951, Ibid., 206, 242. 
KirtTEL, C., 1949 a, Rev. Mod. Phys., 21, 541; 1949 b, Phys. Rev., 76, 743. 
Kirrei, C., and Lurrincer, J. M., 1948, Helv. Phys. Acta, 21, 480. 

Lanpau, L., and Lirsuitz, E., 1935, Phys. Z. Sowjet, 8, 153. 

Mer, C. D., 1950, Proc. Phys. Soc. A, 63, 922. 

Oxamura, T., Torizuka, Y., and Kojima, Y., 1951, Phys. Rev., 84, 372. 
Scort, G. G., 1951, Phys. Rev., 82, 542. 

SHULL, C. G., STRAUSSER, W. A., and WoLLaN, E. O, 1951, Phys. Rev., 83, 333- 
SMALLER, B., 1951, Phys. Rev., 83, 812. 

STANDLEY, K. J., 1949, Thesis, Oxford University. 

STEWART, K. H., 1951, Proc. Phys. Soc. A, 63, 761. 

STREET, R., 1951, Sci. Prog., 39, 258. 

VaN VLECK, J. H., 1950, Phys. Rev., 78, 266. 


Paramagnetic Resonance and Spin-Lattice Relaxation 
Times in the Copper Tutton Salts 


By D. M. S. BAGGULEY anp J. H. E.-GRIFFITHS 
The Clarendon Laboratory, Oxford 


MS. reccived 14th Fanuary 1952 


ABSTRACT. Paramagnetic resonance absorption in single crystals of the copper 
‘Tutton salts has been investigated at 3:10 cm wavelength between 90°K and room 
temperature. In salts diluted with zinc Tutton salt the hyperfine lines are clearly resolved 
at 90° K but broaden and merge into one another as the temperature is raised. From 
these measurements the spin-lattice relaxation time 7 can be estimated and is found to be 
about 2x10-*sec. It is found that 1/7 varies linearly with temperature between 
150° k and 290° x. _An investigation of the crystalline field in the Tutton salts shows 
that for the concentrated salts the field has a strong rhombic component above 150° k, 
whilst in the salts diluted with zinc to less than 1° copper the crystalline field has 
tetragonal symmetry at 290° k, but the tetragonal axis lies in a plane perpendicular to that 
containing the tetragonal axis at 90° K. 


§1. INTRODUCTION - 
i experiments of Gorter and his collaborators (Gorter 1947) have 


enabled spin-lattice relaxation times to be measured for many salts of 

the iron group between 1°k and 90°, and a theory which is in broad 
agreement with the experimental results has been formulated. At present, 
however, there are few results with which to compare the theoretical predictions 
for the temperature range 90°K to 300°k. It is the purpose of this paper to give 
an estimate for spin—lattice relaxation times in this temperature interval which 
have been derived from measurements of line widths in the paramagnetic 
resonance spectra of magnetically dilute copper Tutton salts. 

In crystals containing less than 0-1°% copper the hyperfine structure is clearly 
resolved at 90°K, but as the temperature is raised the absorption lines gradually 
merge into each other until at room temperature only a single broad absorption 
line is observed. ‘l'he magnetic interaction between the copper ions is negligible 
for these dilutions, and it is improbable that the variation in line width can be due 
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to a change in position of the neighbouring protons since the minimum line width 
measured at 90°K is in good agreement with theoretical estimates based on x-ray 
measurements made at room temperature. ‘The broadening of the hyperfine 
absorption lines is interpreted therefore in terms of a variation in spin-lattice 
relaxation time when the temperature is raised from 90°K to room temperature. 

With such broad lines and dilute samples the resonance absorption at room 
temperature is very weak, and so the experimental technique is difficult. The 
measurements are further complicated by a change in the crystalline field so that 
the interpretation involves assumptions about both the hyperfine splitting and 
the line shape. However, it does not seem possible to fit our experimental 
curves if the relaxation time decreases more rapidly than 1/T between 150°K 
and 300° k, and we believe that the magnitudes we have estimated are sufficiently 
accurate to provide a reasonable check on any theory of relaxation times in this 
temperature range. 

The crystalline electric field in the Tutton salts was determined, in so far as 
Was necessary, by measurements at a wavelength of 3:10cm for temperatures. 
between 90°K and room temperature. For the concentrated salts our results. 
are in good agreement with the susceptibility measurements of Krishnan and 
Mookherji (1938), but indicate that the crystalline field has a component of 
rhombic symmetry at all temperatures. In crystals diluted with the corresponding 
zinc ‘Tutton salt the crystalline field appears to be very nearly tetragonal at room 
temperature, but the axis lies in a plane perpendicular to that containing the 
tetragonal axis at 90°K. For intermediate temperatures the crystalline field has 
rhombic symmetry. ‘These measurements are discussed further in § 4(ii). 


§2. EXPERIMENTAL ARRANGEMENT 


In order to observe paramagnetic resonance it is usual to place the sample 
in a cavity resonator situated between the poles of a magnet and measure the 
change in Q of the resonator as the applied magnetic field is varied, the incident 
frequency meanwhile remaining constant. A change in Q of less than 1°% cannot 
easily be detected by the direct method of measuring the power transmitted 
through the resonator using a silicon tungsten crystal detector (Bagguley and 
Griffiths 1950), and for improved sensitivity one of two other methods is usually 
employed. When the absorption lines are sufficiently narrow it is convenient 
to extend the technique of field modulation employed in nuclear resonance 
experiments and to sweep the magnetic field completely through the line or to 
measure its differential. ‘To apply this method to broad lines, requiring a field 
sweep of a few hundred gauss, is difficult and not easily achieved at a frequency 
which will pass through an a.c. amplifier. Alternatively a ‘magic tee’ or a hybrid 
ring circuit may be used (England and Schneider 1950). We have found it 
more convenient to use a bridge circuit in which the steady signal transmitted 
by the resonator when no absorption takes place is balanced by a similar signal 
from a reference resonator. ‘The reflector voltage of the 3cm klystron is 
modulated at 1 Mc/s and the resulting frequency modulation traverses the 
response curves of the resonators once in each half cycle of the impressed voltage. 
The power transmitted by each resonator therefore has a 2 Mc/s amplitude 
modulation, and after detection by silicon tungsten crystals the signals are fed 
in antiphase to a balanced transformer tuned to this frequency. ‘The system is 
balanced when there is no paramagnetic absorption, and the unbalance within 
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an absorption line is amplified using a narrow-band 2 Mc/s amplifier and 
displayed on a cathode-ray oscilloscope. 

With two resonators in such a bridge circuit it was hoped to balance out 
frequency and amplitude fluctuations due to the reflex klystron, and by detecting 
a 2Mc/s signal a reduction in crystal noise was expected so that a considerable 
increase in sensitivity might be obtained. In fact the sensitivity has been 
improved by a factor of about one thousand compared with the simple crystal 
detector and galvanometer method, but it appears that the method has 
limitations which prevent it from attaining sensitivities which should in principle 
be available to microwave spectroscopy. The balanced arrangement is perfect 
only when the resonators have the same Q value, and this cannot be achieved 
when the resonator containing the paramagnetic sample is cooled through a 


temperature range in which the dielectric loss of the sample undergoes a_ 


considerable change. Further, the use of two crystal detectors introduces 
additional noise voltages, and it seems that at 2 Mc/s this is still a serious factor. 
‘The method is, however, extremely flexible and, being independent of ring 
circuits and ‘magic tee’ junctions, can be used continuously over a wide range 
of wavelengths. In addition only one reflex klystron is required for the 
spectroscope described in this paper, so that if too high a sensitivity is not 
required the method may be of considerable value in microwave research. 
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Fig. 1. Block diagram of apparatus. 


A block diagram of the apparatus is shown in fig. 1. The reflex klystron K 
(C.V.223) was operated at 3-10 cm and produced about 50 mw power. Sufficient 
stability could only be obtained by supplying all the potentials of the klystron 
from batteries. ‘The reflector voltage was modulated with an amplitude of 40 v 
at 1 Mc/s from a quartz crystal controlled oscillator (Montgomery 1947). Power 
from the klystron was fed into the resonators R, and R, viaa ‘'T’ junction. Each 
resonator was excited in the H,, mode, R, was coupled to the waveguide by coaxial 
line and probe, and Ry was directly coupled by aniris fin. in diameter. The O of 
each resonator was about 5000 at 3-10cm. ‘The power transmitted by each 
resonator could be varied independently by the attenuators A, and A, and, 
after detection by the silicon tungsten crystals D, and D,, the resulting 2 Mc/s 
signals were fed in antiphase to the mutually perpendicular windings of the 
balanced transformer B. ‘The crystals were matched to the transformer by a 
tapping and the input and output circuits were tuned. The output from the 
transformer was amplified by a narrow-band crystal amplifier tuned to 2 Mc/s 
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with a gain of about 120db. During an experiment the bridge was balanced 
initially using the attenuators A, and A, for amplitude adjustment, and phase 
changes were controlled from the tuned circuits of the transformer. It was. 
possible to balance to better than 1 part in 104. Paramagnetic absorption 
decreased the power through R, and the out of balance signal was amplified 
by the 2 Mc/s amplifier. The output from this amplifier was detected and fed 
through a low-pass filter with a cut-off at about 100c/s and operated a diode 
switch which impressed the absorption curve on a 500c/s sine wave. After 
further amplification this signal was detected and applied to the Y plates of the 
oscilloscope. A 50v bias was applied to the 500 c/s switch so that no signal was 
visible until a 50 v unbalance was registered at the output of the 2 Mc/s amplifier. 
This was effective in removing the residual unbalance of the bridge and enabled 
the magnitudes of the two signals to be so adjusted that paramagnetic absorption 
further increased the output voltage applied to the oscilloscope. 

The magnet was excited from the 100 v d.c. mains and produced a maximum 
field of 6500 gauss. By means of a motor-driven rheostat a periodic field sweep 
was obtained having an amplitude variable between zero and 4000 gauss and a 
period variable between 10 and 20seconds. The field was calibrated using a 
fluxmeter and standardized search coils to an accuracy of within about 1%. 
A rotating search coil with slip rings and carbon brushes provided a voltage 
proportional to the magnetic field. The coil rotated at 3000r.p.m. and the 
output voltage after amplification and detection was applied to the X plates of the 
oscilloscope. ‘The absorption lines were photographed using Kodak recording 
paper R.P.30. 


§3. CRYSTALLOGRAPHY 


A detailed discussion of the crystallography of the copper Tutton salts in 
relation to paramagnetic resonance has been given previously (Bleaney, Penrose 
and Plumpton 1949) and we shall use the notation and designation of magnetic 
axes given by these authors. The crystals are monoclinic, the dimensions of 
unit cell for the potassium salt being a, b, c=9-1 A, 12-2 A, 6-14, whilst for other 
salts the axial lengths are slightly greater, corresponding to the increase in ionic 
volume. Unit cell contains two inequivalent Cu?" ions, the one at (4, 3,0) being 
obtained from that at (0,0,0) by translation followed by reflection in the 
(0, 1,0) plane. Each ion is surrounded by six water molecules situated at the 
corners of a distorted octahedron, the mean Cu—H,O distance being about 2A. 

Single crystals of the potassium, ammonium and rubidium Tutton salts 
were grown by slow cooling of a saturated solution. ‘These were diluted with 
varying amounts of the corresponding zinc salt so that mixed crystals having 
between 1°%, and 0-05°% of the copper salt could be used. For the determination 
of the crystalline field symmetry and the g values, plane sections were cut, 
1 or 2mm in thickness and about 0-5 cm? in area, containing the magnetic axes, 


i.e. the planes K,K,, K, Kg, K,Ks. 


§4. RESULTS AND DISCUSSION 
(i) Introduction 


Measurements were made at 3-1cm wavelength at temperatures between 
90°k and room temperature. From the values of magnetic field H at which 
resonance absorption occurred, the g values were determined using the relation 
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hy =gBH where v is the frequency and 8 the Bohr magneton. The results of ; 
the measurements at 290°K are summarized in table 1. The g values are in agree- — 


ment with those derived from the susceptibility measurements on the concentrated 
salts made by Krishnan e¢ al. (1933, 1938) at room temperature. The 
appropriate diamagnetic correction to the susceptibility measurements was 
calculated and found to be —200x10-8c.g.s. e.m.u. per gram mol. The 
accuracy of the g values is between 1°% and 2%. 


Table 1 
Bat eee a p o g Lo ca 

K,Cu(SO,),.6H,O Rhombic 1052 — 2°31 2:07 2-25 
(102°) (2-33) (2-06) (2-28) 

(NH,)sCu(SO,)».6H,O Rhomibic i ee ES 2-09 2-25 
: (77°) (2:32) (2-09) (2-28) 

Rb,Cu(SO,),.6H,O Rhombie”  1dse ) ee Se ae 2-24 
105°) (2°31) (2-06) (2°23) 

K,(Zn-+ Cu)(SO,),.6H,O* Tetragonal 1152 322 2:26 2:06 2:24 


(NH,).(Zn-+ Cu)(SO,)..6H,O* Tetragonal 167° B2a 2:28 2:09 2:24 
Rb,(Zn+ Cu)(SO,)..6H,O* ‘Tetragonal ily Sar 


* These salts were diluted with zinc Tutton salt to less than 1% copper. 
The figures in brackets are derived from Krishnan et al. (1933, 1938). 
The angles %, « are as defined by Bleaney, Penrose and Plumpton (1949). 


Table 2 
(Cu+Zn)K,(SO,),.6H,0 (Cu-+Zn)(NH,).(SO,),.6H,O (Cu+Zn)Rb,(SO,). . 6H,O 
TS) enh) (2) (3) (4) (1) (2) (3) (4) (1) (2) (3) (4) 
90 220 230 210 

125 200 18 12 4-2 

163 180 23 17 3-0 

170 195 25 19 2:7 

173 180 26 20 2°5 

210 195 26 20 2:5 

215 2 180 29 7483 2:2 

240 180 30 24 2:1 

290 150 36 30 1:7 150 33 D7; 1:9 150 35 ZT 1:9 


Half width of complete 
line observed at 290° x is 
100 gauss. 


Half width of complete 
line observed at 290° K is 
90 gauss. 


Half width of complete 
line observed at 290° k is 
90 gauss. 


(1) Hyperfine splitting (gauss); (2) half width of individual lines (gauss); (3) AHs (gauss); 
{4) 7 (sec x 10°), 


Typical photographs of the absorption spectrum of the potassium and 
rubidium salts containing less than 0-1% Cu are given in fig. 2 (Plate). These 
show the variation in width of the hyperfine lines as the temperature is raised 
from 90° to room temperature. It can be seen that the absorption lines which 
are Clearly resolved at 90° K merge into one another as the temperature increases. 
Measurements were made with the magnetic field normal to an (0, 1,0) plane 
since this direction must remain a principal magnetic axis independent of 
crystalline field changes, and hence the spectrum can contain only four hyperfine 
lines. In table 2 the hyperfine splitting and line widths obtained from such 
photographs of the absorption spectrum are given together with the derived 
values of spin-lattice relaxation time. 
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(ii) Crystalline Field 


It has been reported previously (Bleaney, Penrose and Plumpton 1949) 
that at 90°K for concentrated copper Tutton salts the crystalline field is 
predominantly of cubic symmetry with an additional tetragonal and smaller 
rhombic component. We have observed two absorption lines in both the K,K, 
and K,Ky planes of the K,, (NH,). and Rb, concentrated copper Tutton salts 
at room temperature. If the crystalline field had tetragonal symmetry we should 
have observed a single absorption line for all orientations in one of these planes, 
namely, that in which the magnetic field direction made equal angles with the 
tetragonal axes of the two ions, and the fact that we could resolve two lines in each 
plane indicates that a large rhombic component of crystalline field must be 
present. * 

In fig. 3 (Plate) we give typical absorption curves for the concentrated 
potassium salt at 45° to K, in the K,K, and K,K, planes. The maximum and 
minimum g values in these planes are (for the potassium salt): K, Ks, 2max= 2°40, 
eae tT 7 es Came 292, Lin 2 04. , 

For crystals diluted with zinc Tutton salt to less than 1°% copper the 
crystalline field is very nearly tetragonal at 290°K. ‘The absorption line at room 
temperature has a half width of about 100 gauss, so that the measurements are 
not extremely sensitive to small departures from tetragonal symmetry, but 
only a single absorption line is observed at 45° to K, in the K, K, plane (a direction 
in which any rhombic field would produce maximum separation of the absorption 
lines), and also there is very good agreement between the values of g, measured 
in the K,K, plane and along K,. These values are given in table 3. 


Table 3 
K.K, plane K, 
Salt Si 81 8 
(Zn+ Cu)K,(SO,)..6H,O 2:05 DMS 2:26 
(Zn-+Cu)(NH,)(SO,)..6H,O 2-04 2-26 2-28 
(Zn+Cu)Rb,(SO,),.6H,O 2:08 BAG PLES} 


The tetragonal axis, at room temperature, lies in the K,K, plane at an angle 
of about 30° to K, (see table 1). ‘This is approximately at 70° to the direction 
of the tetragonal axis at 90° K, which is in the K, K; plane inclined to K, at an angle 
of about 40°. There is a considerable change in crystalline field, therefore, as 
the crystal is cooled down to 90°k, but this does not result in a pure cubic field 
at any temperature, there being a strong rhombic component present at 
intermediate temperatures. 

Polder (1942) has given a theoretical account of the magnetism of the copper 
ion in the Tutton salts assuming the crystalline field to be of cubic symmetry 
with a small tetragonal perturbation. He remarked that the susceptibility 
measurements (Krishnan ef al. 1933, 1938, Hupse 1942) did not distinguish 
between two possible orientations of the tetragonal axis, and he arbitrarily chose 
that orientation which gave reasonable values for the splitting of the orbital 
levels and the susceptibility of the individual Cu?+ ions when the tetragonal 

* Arnold and Kip (1949) have reported measurements on CuK,(SO,),.6H,O at 290° K and 


conclude that the crystalline field is tetragonal. These authors, however, examined the K,K,; plane 
only, and this is insufficient to exclude the possibility of a rhombic crystalline field. 
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field had a positive sign. He showed that the g values parallel and perpendicular 
to a tetragonal axis were given by 


4X A 
n=2(1- a 4,=2(1- eh 


Here A is the spin-orbit coupling coefficient; F,, F;, F,, Ff; are the levels into 
which the 2D state is split by the crystalline electric field, and of these F is the 
ground state whilst F,, F are components of the triplet in a cubic electric field. 


ee ee ee 


ewan 


Polder assumed that the tetragonal axis was in the K,K, plane inclined — 


at 40° to K, and obtained good agreement between theory and experiment 
using the values g,=2-44, g,=2-:05, A=—852cm™4, F,—F,=15 400 cm4, 
F,—F,=26600 cm!. The experiments reported in this paper, however, show 
that whilst Polder’s assumption is valid at 90°K it is not correct at room 
temperature.. Using the values for g, and g, derived from the paramagnetic 


resonance experiments at 290°K, i.e. g,=2-05, g, =2-26, the values for P,— Fs ~ | 


and F,—F, are found to be 135000cm™ and 6500 cm respectively, so that, 
on Polder’s theory, F lies below F, and the splittings are of widely differing 
magnitude. It is usually assumed that the crystalline field arises from the 
octahedron of water molecules surrounding the copper ion, and under these 
circumstances we should cxpect the values of F,— F, and F;—F, each to be 
about 20000cm™?. ‘These values are confirmed by optical measurements on 
solutions of cupric salts which show a broad absorption band at about 12300 cm? 
and no further absorption below 32000cm™+ (Abragam and Pryce 1951b). 
A splitting of this order of magnitude can give rise to the observed g values if 
tetragonal terms opposite in sign to those given by Polder are used (Polder 1942, 
Abragam and Pryce 195la). We then obtain g,=2-0, g, =2—6A/AE, where 
AEF is the splitting of the ?D state in a cubic field and hence, for g, =2:26cm™ 
and A= —828cm™+ (Abragam and Pryce 1951b) we obtain AE=19000cm'1. 
The agreement between experiment and theory is now much more satisfactory ; 
the splitting AZ is of the expected order of magnitude, and the value for g, =2-0 
approaches the experimental result of 2-05. The small difference may be due 
to a weak rhombic component of the electric field. We conclude therefore that, 
in crystals of the copper Tutton salts diluted with zinc to less than 1° copper, 
the crystalline field at 290° is, to a good approximation, tetragonal, but with a 
negative sign, and the tetragonal axis lies in the K,K, plane making an angle of 
about 30° with Kg. 


(iii) Hyperfine Structure 


Abragam and Pryce (1951a) have given a theoretical interpretation of the 
hyperfine structure for a Cu?+ ion in a crystalline electric field of tetragonal 
symmetry. ‘he Hamiltonian (in their notation) is 


W =, PHS, + 8 BH S_+ H,S,)+ ASI, + BUS ol a+ Syly) 
+ QUP— 31+ 1)}—yByH 
In this equation the last two terms are small and may be neglected for this 
discussion. When the steady magnetic field is directed parallel and perpendicular 


to the tetragonal axis the overall separation of the hyperfine lines is 3A and 3B 
respectively, whilst for any arbitrary orientation of the crystal the separation 


i 
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will not be less than the smaller of these values. For a tetragonal field the 
constants 4 and B are given by 


2 = 2 
A= EPS (54 rt-gy—2), B= PEER ( 48+ g,—2) 


in which the notation is again that of Abragam and Pryce (195la). For the 
dilute copper Tutton salts, in which the tetragonal component is negative, the 
lower sign is to be taken so that, writing 2yBBy/r?=0-025, «=0-28, g,=2-05, 
&,=2-26, we obtain |A|=85 x 10-4cm™, |B|=78 x 10-+cm-4, from which the 
overall splitting of the hyperfine lines corresponding to 3|A| is 270 gauss, and 
that corresponding to 3|B| is 220 gauss. ‘The hyperfine structure is therefore 
almost isotropic, and in fact the calculated value of the hyperfine splitting when. 
the magnetic field is along the 6 axis of the crystal is 225 gauss. 

It will be seen from table 2 that we have given the value for the overall 
hyperfine splitting along the } axis at room temperature as 150 gauss, and this is 
considerably less than the theoretical estimate of 225 gauss. It is necessary 
therefore to indicate how the experimental value has been obtained. When 
constructing the absorption spectrum theoretically the line shape was first 
determined (see §4(iv) below) and the line width and hyperfine splitting then 
adjusted to give the best fit with the experimental curves. Below 250°K it is 
possible to distinguish the peaks of the individual absorption lines; the hyperfine 
splitting is therefore known fairly accurately and the width of the lines is adjusted 
to reproduce the observed spectrum. At room temperature, however, the 
absorption lines are completely merged together so that the choice of hyperfine 
splitting is somewhat arbitrary. We have chosen a value which gives results 
most consistent with measurements at the lower temperatures both for the 
hyperfine splitting and the line width. In any case it seems unlikely that 
the hyperfine splitting can be as great as 225 gauss at room temperature since 
the total width of the absorption line at half power is only about 200 gauss. If the 
splitting were to be as great as 225 gauss the separate peaks would be resolved 
unless there has been a considerable broadening of the absorption lines, and 
such a broadening would of course further increase the width of the absorption 
band beyond 225 gauss. It may be that the term 2y8fy/7° or « has changed, but 
these are not expected to be temperature dependent and are taken from the 
accurate measurements of the hyperfine structure described by Bleaney, Bowers 
and Ingram (1951). 


(iv) Line Widths 


In the copper Tutton salts diluted with zinc to less than 1% copper the 
contribution to the line width from dipolar coupling between the Cu?* ions is 
negligible, and so the line width arises from spin-lattice coupling and the magnetic 
interaction of the Cu?* ions with the protons in the neighbouring water molecules. 
It is usual to assume that the complete width is given by simple addition of these 
two contributions so that 6 gauss, representing the proton contribution, must 
be subtracted from the measured line width in order to derive the effect of 
spin-lattice interaction. ‘The line shape should be a combination of the gaussian 
and damped harmonic oscillator equations, but we have found sufficient 
agreement by replacing this bell-shaped curve by that for a damped harmonic 
oscillator alone. In order to construct the spectrum we assumed, therefore, 
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a line shape of the form = 1,(1-+.?), in which » is the distance from resonance. | 
measured in units of half the width between the half-power points, i.e. in units: 


of the half width of the absorption line. 


From the half widths of the absorption lines so constructed as to give best | 
fit with the experimental curves the spin—lattice relaxation time was derived | 


using the relation 7=1:14x 10-7/gAH,. Here g is given by hv=gPH at the 
centre of an absorption line, AH, is half the line width between half-power points 
minus the 6 gauss mentioned above, and 7 is the relaxation constant as defined 
by Van Vleck (1940). The relaxation times for the different copper salts are 
plotted in fig. 4. It can be seen that they lie on a single curve and extrapolate 


back reasonably well to the point at 90° k which is taken from the measurements _ 
of Bleaney, Bowers and Ingram (1951). Above 150°x the points on the (1/7, T) 


curve lie approximately on a straight line, so that in this temperature interval 
the relaxation time does not appesr to be changing more rapidly than 1/7. 


i D Bleaney et ad. (1951) 
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Fig. 4. Spin-lattice relaxation time as a function of temperature. 


The theory of spin-lattice interaction distinguishes two coupling mechanisms. 
The ‘direct’ process, in which a magnetic ion can interchange energy with the 
lattice vibrations only by a simple quantum exchange, gives rise to a relaxation 
time proportional to 1/7, whereas if Raman scattering of lattice waves can take 
place, the relaxation time should be proportional to 1/77(T’<0) or 1/T?(T>6@) 
where @ is the Debye temperature (Kronig 1939, Van Vleck 1940). 

At liquid helium temperatures the relaxation times for the concentrated 
copper ‘T'utton salts vary more slowly than 1/7? (Cooke, private communication), 
so that the ‘direct’ process appears to predominate below 4°K. The value of 7 
at 4-015" for CuK,(SO,)..6H,O is about 9-5 x 10-8sec in a field of 675 gauss 
(Bijl 1950). Measurements at 77° and 90°K have been reported by Broer and 
Klemperman (1947) for CuK,(SO,),.6H,O and Cu(NH,).(SO,)..6H,O in a 
field of 2400 gauss, which is comparable with the field strength used in a 


ey ee 
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paramagnetic resonance experiment. ‘The values of + obtained from these 
authors are 


HTK 90° k 
Cu(NH,).(SO,)2..6H,O 10-7 x 10-8 sec 6:4 x 10-8 sec 
CuK,(SO,),.6H,O 14:3 x 10-8 sec 8:6 X 1078 sec 


Another estimate of the spin—lattice relaxation time, derived from the width of 
hyperfine lines in paramagnetic resonance, has been given by Bleaney, Bowers 
and Ingram (1951). At 90°K the relaxation time is reported to be about 
5x10-8sec, in good agreement with the dispersion measurements on the 
concentrated salts. ‘Thus the relaxation time is varying rapidly with temperature 
between 4°K and 90°x, and it is thought that the Raman process is predominant 
in this temperature interval. Our measurements indicate that a striking change 
occurs about 150°K. In the 140° interval between this temperature and 290°K 
the relaxation time decreases only by a factor of two, whereas in the 60° interval 
above 90° the relaxation time has changed by a factor of twenty. It seems 
unlikely that the linear variation of 1/7 with T above 150° can be adequately 
described by the present theory of Raman scattering, which gives a relaxation 
time changing at least as rapidly as 1/7. 
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ABSTRACT. One obtains two different theories of semiconductors, depending on whether 
the free energy due to the spin of electrons in donor (or acceptor) centres is taken into account 
(approach A) or not (approach B), It is shown (i) that the former approach can be carried 
through for a degenerate semiconductor, (ii) that it applies formally to the case when the | 
interaction between impurity centres and conduction electrons is taken into account, 

provided one of the symbols is reinterpreted, and (iii) that it is to be preferred. Approach 
B has, however, been used almost universally. It is shown that it corresponds to the 
unlikely situation in which fifty per cent of the electrons from each impurity level come 
from paired states, and the remaining fifty per cent from unpaired states having the same 
activation energy. | 


| 


§1. INTRODUCTION 
T does not seem to be realized that two models are in use in current 
semiconductor theory. ‘The first 1s essentially due to Mott and Gurney ~ 
(1948, p. 157) and depends on minimizing the free energy. It will be 
referred to as A. It has not so far been formulated so as to take the Fermi—Dirac 
distribution of current carriers in the conduction band into account rigorously, 
and has, therefore, not yet been applied to degenerate semiconductors. ‘The — 
second approach, to be referred to as B, is essentially Wilson’s original approach ~ 
(Wilson 1931, Seitz 1940, p. 186). It is known to hold also for degenerate 


semiconductors, but it cannot take into account the contribution to the free 
energy due to the unpaired electrons which remain behind in impurity centres 


at any given temperature. 

We shall show that approach A can be extended to degenerate semiconductors, 
and that the equations obtained differ markedly from those obtained by 
approach B. The present form of approach A is always to be preferred to approach B. 


§2. EXTENSION OF APPROACH A 
The simple generalization of Mott and Gurney’s argument for approach A 
is as follows. Consider a crystal of an excess semiconductor having donor 
impurity levels at the energies Ej(j=1, 2, .. .), which are generally negative if 
the zero of the energy scale is taken at the bottom of the conduction band. The 
free energy is 


A= —Xn,E; + (Xjn,)RTn* —8RT 9 F'30(n*) —RT log P—g, ......(1) 
where n, is the number of electrons donated by the jth level, and the summations 
are over allimpurity levels. ‘The first term represents the energy to bring electrons 


into the conduction band. ‘The next two terms come from the theory of a gas 
of n electrons, according to which its free energy is 


RT [nn * — 39 F 319(7*)], where n= “(2m *RT)32 and 


Re 
me aK 
)=[- 1+exp (w—7*)* ee 
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7* is the Fermi level ¢ divided by RT, m* is the effective electron mass. ‘The 
last term but one gives the free energy of the vacant places from which electrons 
have been removed. Let N,; be the number of available donor electrons and 
N;—N,; the number of electron vacancies in the jth energy level at the absolute 
zero of temperature. Then there will be N;—N,;+,; electron vacancies at 
higher temperatures. P is the number of ways in which N,—N,,+, vacant 
places can be distributed among N, sites, N, —N.. +n places among N,j sites, etc. 
‘Thus, using Stirling’s approximation, 

0 log P Nex— 

ma? SEN ee) es nee (3) 
‘The last term of (1) is some function of the N,’s, N,;’s and ,;’s which is to take 
account formally of the contribution of the spins of the electrons to the free 
energy, of the effect of interactions among the electrons and impurity centres and 


of any other terms which have not been included explicitly. For equilibrium 
{0A/dn,),=0, so that by (1) to (3) 
Wi Neon rte og So we 
ee ee a, Sa Oa (a PR e ekeienels (4) 


Here we have used the results 


dF, . “ y * 
Fee =eFealn®) (c=0) and Bayp=mFialn*) 
Hence the number of electrons in the level F, at temperature T is 
N 
Nop Mp = moose ts 5 
1+ 8, exp (7 =4") ( ) 
ie 1 og 
where p-=exp aa | airs: (6) 
and 7,=£,/RT. ‘The Fermi level is determined by 
N 
x la oo PP ln*). iw ee ee 
»| Ne- eas | Fino) (7) 


If g takes account only of the spin contribution to the free energy then any 
quantity f, has the value } if none of the N., sites of the rth energy level has paired 
electrons (two electrons of opposite spin) at the absolute zero, while it has the 
value 2 if all NV, sites have paired electrons, since the argument follows then exactly 
along the lines laid down by Mott and Gurney. (More generally, it may readily 
be shown that if a fraction f, of the u, electrons donated by the level £, comes 
from paired sites, then B,,=27/-1.) Approach B is obtained if 8,=1 (for all r). 
‘This corresponds to the case in which fifty per cent of the electrons which come 
from the 7th level come at all temperatures from paired states, and fifty per cent 
from unpaired states. ‘This situation is most unlikely as electrons in paired states 
have in general a lower activation energy than electrons in corresponding 
unpaired states. It follows that approach B 1s less satisfactory than approach A. 

The value of £, has an important effect on the relation between the parameters 
of the theory. We merely show (in the figure) the degree of ionization of a single 
impurity level at room temperatures as a function of the impurity concentration 
for B,=4, 1,2. ‘The figure shows incidentally that, for a fixed activation energy, 
impurity centres with paired electrons are more readily ionized than impurity 
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centres with unpaired electrons. This is entirely due to the effect which the : 
spin of the electrons left behind in the impurity centres has on the minimum of a 
the free energy. A wrong choice of 8, in the analysis of experiments can lead al 
to an incorrect activation energy. In particular, if Z, is independent of 
temperature, a wrong value of 8, may lead one to infer an incorrect and temperature-_ 
dependent activation energy, which would agree with the true activation energy 
only at the absolute zero of temperature. 


Paired Electrons, 4 =2 


og Approach B, , =I 


c Unpaired 
i=] 

S 06 Electrons, 
N Beh 
s 

je 
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7, (cm) 
The effect of the two different approaches on the degree of ionization of the impurity centres. 


A simple excess semiconductor at room temperature is assumed, having a donor level 
0:1 ev below the conduction band and having N,= Ney. 


§ 3. SPECIAL CASES 

‘The above results can be extended to take the full band and acceptors into. 
account. We shall, however, merely indicate the modifications brought about 
by approach A in the usual formulae for some important special cases. Suppose: 
only one impurity level is present. The mass action law for the ‘reaction’: 


excess electron + ionized donor = neutral donor 
is easily seen to be 
n,(N, — Ne, +71) 
AT = Bite M ali ®) exp (nrg pce (8) 
el 1 


K being the reaction constant. Using the classical approximation 


Fig") SE 7 CkD | a ee eee (9) 
and putting N,=N,;, 8, = 4, one obtains from (8) a formula which has been 
used by Pearson and Bardeen (1949)} in the analysis of their experimental. 
results, which is therefore based on approach A. Most other work on. 
semiconductors is, however, based on approach B. 

If N.1=N,, (7) can be written in the form 


NV, =[1 +8, exp. (9% ta) ep diay). tee cee (10) 


{| When comparing their equation with (8), their value of Na must be put equal 
to E 
possibility of having acceptors present has here been excluded. ar aa ere a 
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If 8, is replaced by unity, (10) yields a result which is well known in theories 
based on B. From (10) follows for given EF, a relation between the maximum, 
-o*; Say, of »* and the temperature for which it is reached : 


no* =Byn — log (—§Pym, — 1). oS ee, 


For £, =1 one obtains a result first given by Shifrin (1944), using approach B. 
From (10) there also follows the well-known approximate expression for the 
Fermi level, provided (9) can be used and exp (n*—7,) >1: 


nae ee 
Q= puts Es Toe ak ee ee allZ) 


These modifications must be borne in mind in the usual expression for the 
electrical conductivity and the various coefficients which arise in second-order 
effects (Hall effect, thermoelectric power, etc.: for a list of formulae, see R. W. 
Wright 1951, approach B) wherever the formulae involve the Fermi level or the 
impurity level. The reason is that the former is affected by the value of B, and 
the latter occurs usually in the improved approach A multiplied by f, (cf. the 
equations of §2). We merely note here that if the mean free path / is independent 
of the energy and inversely proportional to the temperature, and we put 
l=I,/RT, then the electrical conductivity is given by (Shifrin 1944) 


o=ogsoe( ia-expy*), oy lOmm*tel 3h) css: (13) 
At low temperatures (12) and (13) lead to 


oy (2N,h?)12 E, +0g/on 
log (oT?!*) = log E a | aa ae Disilenatonnrte, 


If the term g is due only to the spin of the electrons in the impurity levels, then 
(0g/On,)/RT is independent of temperature. Thus, if the activation energy has 
been determined from a plot of log (oT**) against 1/7, the value obtained will 
be independent of f§,, and will therefore be the same whether the analysis is 
based on approach A or approach B. If g takes into account interactions, and so 
becomes a function of N, and 7, (14) shows that the experimentally determined 
activation energy must also be expected to depend on these quantities, as 
observed for silicon alloys containing boron (Pearson and Bardeen 1949, 
Castellan and Seitz 1951), cuprous iodide in equilibrium with iodine vapour 
(Vine and Maurer 1951, using approach B), and the oxides of titanium, iron, 
nickel and zinc (Mitchell 1952, using approach B). 


ACKNOWLEDGMENT 


The author is indebted to Miss L. S. Souter of the Research Laboratories 
of the General Electric Company Ltd., Wembley, Middlesex, for some 
stimulating discussions concerning the paper by Pearson and Bardeen. It is 
to these that this note owes its origin. 


+ An analytical proof is also known (Landsberg, Mackay and McRonald 1951, approach B) and 
is readily generalized to deduce (11). 

t There is some confusion in the literature with respect to this equation. It should be noted 
that in Mott and Gurney’s book (1948) the ‘(2)’ refers to the case of paired electrons on p. 158 and to 
the case of unpaired electrons on p. 175. In Wright’s book (1950, p. 51) on the other hand, ‘(2)’ 
should be replaced by ‘ 2’ on the right-hand sides of eqns. (26) and (27). ‘The term in the logarithm 
- of eqn. (28) will then be found to be multiplied by the additional factor (2)/2, and this is in agreement 


with our eqn. (12). 
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ABSTRACT. 'The quantization of the complete gravitational field is carried out by 
extending the work of an earlier paper. ‘The main obstacles in the quantization of Einstein’s 
field are overcome by expressing the field quantities in the Riemannian space as expansions 
in the flat space, and then splitting the gravitational field into the linear and the non-linear 
parts. The linear part of the gravitational field is regarded as the free gravitational field, 
while the non-linearity is treated as a direct interaction between the gravitons. This 
treatment is quite general, but it suffers from the usual limitations of the perturbation 
method. 

The gravitational self-energies of the photon and the electron are also investigated by 
using the linear approximation of the gravitational field. It is found that the photon 
self-energy vanishes unambiguously, while the electron self-energy is quadratically 
divergent. 


§1. INTRODUCTION : 

HE quantization of the approximate linear form of Einstein’s gravitational 

| field has been carried out in an earlier paper (Gupta 1952). We shall now 

consider the quantization of the exact non-linear field. Einstein’s 

gravitational field differs from other fields in two respects: it involves the use of 

the Riemannian space, and it is non-linear. We shall overcome these two 
obstacles in the following way. 

As usual, we shall put gq!” =” —x«y"", where e'” is the Minkowskian metrical 


tensor given by ‘ 
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and « is a constant. ‘Then we can regard the Minkowskian space as the zero-th 
order approximation to the Riemannian space. Moreover, by expanding all 
the field quantities in powers of x, it is possible to express them as expansions 
in the Minkowskian space. This procedure does not involve any departure 
from Einstein’s ideas. However, we shall differ slightly from Einstein in the 
choice of the coordinate conditions, and use those given by De Donder (1921), 
Fock (1939) and Papapetrou (1948). 

After passing over from the Riemannian space to the Minkowskian space, 
we shall split the gravitational field into the linear and the non-linear parts. 
Then we shall treat the linear part as a ‘free’ gravitational field, while the 
non-linearity will be treated as a direct interaction between the gravitons. It 
will also be shown in the next section that the splitting of the gravitational field 
into the linear and the non-linear parts may be given a simple physical 
interpretation. 

The present treatment is quite general, although it suffers from the usual 
limitations of the perturbation method. However, the treatment of the 
non-linearity of the gravitational field as a perturbation seems well justified, 
because we know that the linear approximation of Einstein’s theory, which 
includes Newton’s theory of gravitation, accounts for the observed facts up to 
a very high degree of accuracy. ‘This shows that the non-linearity can only give 
rise to secondary effects. 

Throughout this work we shall refer to the earlier paper (Gupta 1952) 
as I, and the equations of I will be referred to as (I, 1), (I, 2) etc. 


§2. COORDINATE CONDITIONS FOR THE GRAVITATIONAL FIELD 


EQUATIONS 
As mentioned above, we shall use the coordinate conditions of De Donder 
(1921), which are given by Mee Oy (1) 


Using (1), Papapetrou (1948) has been able to express the field equations for the 
gravitational field in an extremely simple form. In this section we shall give a 
brief account of his work and discuss its physical implications. 

The Lagrangian density for Einstein’s gravitational field may be taken as 


ee HE GHG) 


Regarding & as a function of q”’ and 9% =dq”"/éx", we get for the canonical energy 
momentum pseudo-tensor density of the gravitational field 


on op le 
aga mie Ween Eo ae (3) 


je 
t= 


We can further obtain the symmetrical energy momentum density by 
Belinfante’s method, as we have done in $5 of I. We consider an infinitesimal 
linear transformation 


on = —e bw, x” with 6@;,,=—OW,,, | — sevnes (4) 
and define a quantity f“”°, given by 
fle — _ bg” OE 
bw,,, O97? 
oh Of apy av Bu = 
= 2a gles —e gr). mete (5) 
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Then the symmetrical energy momentum density 6” will be 


v v 1 0 v A v, 
alg = 4th — 5 ane (EY a fer E Fer), it nee eee (6) 
and the angular momentum density is given by 
HAE 2 = xe HPC — yp” Que } 
=e th ig” ete Le eee (7) 


. 
Adding the symmetrical energy momentum tensor density of the matter field — 
<“& to (5), we obtain for the total energy momentum density 0” of the system — 


v v v 1 0 v, Pa Y, 4 
OM" = 4 (T+ t}) — 5 aye (e eo base de (8) =@ 
In order to simplify (8), we note that (Tolman 1930) 
7] 
fi a eT a pecerse 
T= 2-0, ws (9) 
; on ox 
‘ith RM = He OF Cone eS ee eee 10 
Wi dar? g 3 One? q (10) 
Moreover, according to (5) and (10), we have 4 
fi ea ORY eM are ee (11) 
Using (9) and (11), we can write (8) as 
OY” = ea [eM (ARG” + RB”) + €”P(ARGY + BAG) — e*F(HRG” + HRM). ww. (12) 
Evaluating 3%” from (2) and (10), and substituting in (12), we find 
a 
K2Q" = ere (get Sgt ats gee aye) ee (13) 
If we now use the supplementary conditions (1), we get the simple result 
07 g"” 
a8 = 2 @H” 
axon Ores a le ee ae (14) 


Since the supplementary conditions (1) lead to the simple field equations (14), 
and also reduce to Einstein’s supplementary conditions (I, 12) in linear approxi- 
mation, it is reasonable to expect that they are the appropriate conditions. 

The field equations (14), obtained by using (1), are similar to those for an 
electromagnetic field, the only difference being that, instead of the four-vectors 
A,, and j,, we now have the symmetrical quantities g”” and ©” respectively. Thus 
we may take the field equations for a ‘free’ gravitational field as 

2 ly 
ob PAT cee Oicindin ne ae anne ee (15) 
Ox*dxP 
and regard the entire interaction as caused by ©’, which represents the energy 
momentum density of the whole'system, including the gravitational field itself. 
‘This shows that the choice of the linear field (15) as the free gravitational field 
simply means that we are putting the interaction of the gravitational field with 
its own energy momentum pseudo-tensor on the same footing as its interaction 
with the energy momentum tensor of the matter field. 
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It was pointed out in I that the supplementary conditions (I, 12) are only 
approximately compatible with the field equations (I, 11). But here, on account 
of the conservation equation 


Daren eK ee (16) 
(1) and (14) are exactly compatible. 


§3. FLAT SPACE EXPANSION FOR THE GRAVITATIONAL 
LAGRANGIAN DENSITY 
We shall now expand the Lagrangian density for Einstein’s gravitational 
_ field in powers of the coupling constant x, and show that this expansion can be 
regarded as an infinite series in flat space. For this, as already explained, we put 


i PART at act a Nog See anil, agli ar Pies (17) 
Then, on account of the relation 
a \h 
La” |=(\/—2)* |g? |= Ys - a ag in deat (18) 
2 
we get g=—l1tke,,y + = (Eno €va— Env Ere) ye + O(k?), wee ee (19) 
and ese 
=<? —Ky"” ae de” eypy i. ge yp yh yt 
+ frre en eg, yO ye + hire” cap eag yy + O(k3), oe ee. (20) 


where the tensor «,, is reciprocal to <“”. In order to obtain the expansion for 
Eu» We note that 


; us 
= a ) =—(—2) M2 minora’ ya 6 saat o (21) 
But we easily find that 
minor (g””) Ey + KEE gp — KE ng E vB ye 


A A 
7. Keay EA €Bo yey : i én €Bv Edo yey : 
A 
— dr Ey egpe roy Y? + dK € yy Ear €po Vy Ve + Ok). wees (22) 
Hence we get 


Suv aa Eu ae dE yy Egg Y + Key € yp Vo 


So € au € av € Bo ft he “x Heyy Epy © do ft aM 
+ BK ny Eas C19 VO Ve — fr E ny Egak soy y? + Ol)... 2. (23) 
Now, in order to expand the Lagrangian density (2), we use the relation 
9 =1_~ Bus fs Brn Bur ot ees (24) 
hy es Ont. POxe Se ax" 


and get, after some simplification, 


Ox? dx Ox ax? 


x Lea OL x0 O86 5 | . 


L=-— = g”” 2% Bo 7 Sua OS v0 y) uv O2 oa 
K 


ox! Ox” Ox! Ox” 
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Using the expansions (17), (20) and (23), we find 


1 dyhu Cyr, Oy? dy Bi Oy2e Ay? 

eS 3 | 2 <b Ox” Oak © “eSB xe Ope eR Ox! Ox? 
yoy’ dye oy” @ dy ap OY Au Oye” 
KE ga esp Y" Ox” Ox — BKC A9 fap Y" Ox" Ox? + 2xewat pe” Ox” Ox" 


10 OY? Oy" 
€apSoitneY ax" Ox” 


p00 
nat Pe +0(«2) |. ..(26) 


According to (1) and (17) the supplementary conditions will now be given by 


tKelve 


Oy Cx =U | tetas tea ener (27) 
which enables us to simplify (26). For, the first term in (26) may be written as 
Ory PH Oy 0 pve 0 (oy? 
Pail pia fat pee aa st eae 1 
2 op Ox” Oxt! ax” : (decor Alt | +4 9 ep ax Ox” ceewan (28) 


Since the first term on the right-hand side of (28) is a divergence and the second 
term vanishes on account of (27), we can drop the first term in (26). 

We can write the above results in a more compact form if we use the usual 
flat space notation, i.e. we allow the Greek indices to take the values 
1, 2, 3, 4, instead of 1, 2, 3, 0, where x,=2%). Then it becomes unnecessary to 
introduce the Minkowskian metrical tensors ¢”” and «,,, and all the indices may 
be written as lower indices. Thus, using the flat space notation and dropping 
the first term in (26), as explained above, we may write (26) and (27) as 


— LL OY nr OY 1 OY an Ov ap OYaup OY up 4 OYan OY 8 
EIA Oe Oey OD eee ee Ox, Ox, 2°" Ox, Ox, 
: OY ua. OY vp OY 68 OY ny OY up OY vp 2 
+ 2eyu, Ps Yur aie. oe Dy a Oh, + O(« | Sree (29) 
OY, OK =O. a rer (30) 


It may be noted that the terms independent of « in (29) are identical with the 
Lagrangian density (I, 16). 


§4. INTERACTION OF GRAVITATIONAL AND 
ELECTROMAGNETIC FIELDS 
In order to show how to deal with the interaction of the gravitational field 
with other fields, we shall here consider the interaction of the gravitational and 
the electromagnetic fields. 
In the presence of the gravitational field, the field equations for the 
electromagnetic field are given by (see Tolman 1934) 


. 045 moa a 
waa a Fe oh gO eee (31) 
Of” og" 
re ee I (32) 
where GP a/—gF, Bag A", J=,/—-gJ*. ...... (33) 


The total Lagrangian density for the gravitational and the electromagnetic 
fields may be taken as 


Lig Ech hegee Dee (34) 


eee eer en Pe er 
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where & is given by (2), and 


Za=—4V-2 FY FP, = —4/ go g Big Buy. evens (35) 
Now we have A,=(-g) 1" g,,a° 
=Eyg ee ate eS (od eee RE (36) 
oa’ oa’ o 
so that Fe = Em 5 Fae SHO Baee + KE go, € "os 2 (y? a’) 
— KE gy Eng a (vy? AV OK). anemic (37) 


Using (17), (20) and (37), we may expand (35) as 
Z,=- (eM ee?” —Ke™ > ase — Kee? yy ak dice <*"egy™) 
0a 0a fe 0a 
(Bag San Bigr Sau) ( Baer Se an 
oa ie ur oa" BAe uv 
“Ox! ana Wea» Ox & Ces) a Pes «| pela ada: 
The terms within the square brackets i in (38) may be written as 
Rl cs 
ax" \Oxe 
The first term in (39) vanishes on account of the Lorentz condition 0A“/dx* =0, 
while the second term may be absorbed in O(«?), because <”0? A%/dx" 0x” = O(k). 
Thus (38) may be reduced to the simpler form 
i= He «@” — xe y?” — Kee? y# ae deM ecg) 
oa oa oa oa 
( CaS Aue ) é € Bo — axe cn) ats O(«?). SOS (40) 


Using the flat space notation, as we have done in the last section, we may write- 


(40) as 


+ 


Sax) Rep (e — sas) 2 y" €g€aa + divergence terms. . (39) 


Ze= —4F,, Fy heyap (Pou pu— 808 Pus”) + Ole), . 2.0 (41) 
where, according to our flat space notation, F’,,=0A,/dx,—0@,/0x,. We notice 
that the first order interaction term in (41) is given by —}xy,,7,,, where T,, is 
the usual flat space energy-momentum tensor of the electromagnetic field. This 
result agrees with (I, 54) when we take into account the fact that the trace of 


T,,, vanishes. 


§5. QUANTIZATION OF THE COMPLETE GRAVITATIONAL FIELD 


As already explained we may regard the linear part of the gravitational field 
as the ‘ free’ gravitational field, while the non-linear terms may be treated as 
a direct interaction between the gravitons. ‘Thus, by putting «=0 in (29) we 
obtain the Lagrangian density for the free gravitational field, which is precisely 
_the same as that given by (I, 16). It follows that the entire treatment of 
quantization described in I still holds for the free gravitational field. Among 
other results we can still use the vacuum conditions (I, 40), the commutation 
relations (I, 23) and (I, 24), and the vacuum expectation values (I, 58). 

For the treatment of the complete gravitational field we have to use the 
complete Lagrangian density (29). ‘hen we can pass over to the interaction 
representation in the usual way by treating the linear part as the uncoupled 
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gravitational field and the non-linear part as the interaction. In the interaction 
representation we can again make use of the commutation relations (I, 23) and 
(I, 24) and the vacuum expectation values (1, 58). Similarly, we can deal with 
the interaction of the gravitational and the electromagnetic fields by adding the 
Lagrangian’ density (41) to (29) and treating all the non-linear terms as 
interaction. 

It may also be noted that the supplementary conditions for the complete 
gravitational field in the Heisenberg representation are still given by (1, 56). 
For, according to the field equations (14), dy,,/dx, always satisfies the 
wave equation for plane waves 

rp tus =e Ou =(). 
Ox, Ox, 
Therefore dy,,/0x, can always be split into the positive and the negative 
frequency parts, and the supplementary conditions (I, 56) are maintained in 
course of time. In the interaction representation the supplementary conditions 
can be obtained in the usual way. 

We conclude this investigation with some applications of the present method 
of quantizing the gravitational field. We have seen that according to the above 
method it is possible to treat the complete gravitational field up to any desired 
order of approximation. However, for simplicity we shall use only the linear 
approximation of the gravitational field in the applications to be discussed in 
the following two sections. A more complete account of the applications of the 
present treatment will be given in a separate paper, where the influence of the 
non-linear terms will also be taken into consideration. 


§6. SELF-ENERGY OF THE PHOTON 


We shall first consider the problem of the gravitational self-energy of the 
photon, because the corresponding problem in quantum electrodynamics has 
attracted much attention. As already remarked, for simplicity we shall here use 
only the linear approximation of the gravitational field. 

According to (I, 54) or (41), the Lagrangian density for an electromagnetic 
field interacting with the linear gravitational field is given by 


L= —iF one 3(Fe2) 1 (Fe OV iy 1 oy Oy 


foNire ct | | dt hee ame meant Bef a) 
uv 2\ On, 4\ 0x, Ox, 20x, =) DY op Pal pu — Pup Ll us")- 


Denoting the canonical conjugates of y,,,, y and A, by z,,, 7 and z,, respectively, 
we have 
1 OYy» 1 oy 
Pe) at ) OO a eee 3) @. 8 alee (43) 


a, = — 1B 4+ kY4j Fy, = KY; 5 Fy Be 21K( Yao, =z 2y44) Pia, 
a, =10A,/ OX, 


Obtaining the Hamiltonian density from (42) in the usual way, and expressing 
it in terms of the field variables, their space derivatives and the canonical 
conjugates, we get H=H,+ Hp+ Hi, where Hg and Ay are the Hamiltonian 
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densities for the gravitational and the radiation fields respectively, and Hin, 
is the interaction energy density given by 


Aine on. BKV aol 7? + KY Py Fy — IKY 457%; F;; Z 

+ $4 ( Yee, — 2Y 44) M2? — BY THT e+ BEY, — 2V 44)?" 

— BV 5 (Youn — 24a) MM; + BY HiT 4M] — BUOY (Yo — 244) HP; 

Fas eT EE — HOV aan el mt Ol) wee eee (45) 
The above equation gives the interaction energy density in the Heisenberg 
representation. In the interaction representation the interaction energy density 
will have the same form as (45), but with the difference that now the field variables 
and their canonical conjugates will refer to free fields. ‘Thus, passing over to 


the interaction representation and expressing the canonical conjugates in terms 
of the field variables, we get* 


Fin = 2k Vapl F, ast Cian 9x6 F, -) 
ll ain Bae az 3K 2 V Bye tol ag 13 FV Vp sal ay 
— Fy 4V pal, ieee Fv AV BE nol ua tH Van pal sal up + O(K? A ey eet (46) 


Now according to Dyson (1949) the S-matrix for the interaction of the 
gravitons and the photons will be given by 


S=1+(— [de Hy(s)+(—dP [df de’ PEE), Hyi')} + « 


The terms in the S-matrix contributing to the second-order self-energy of the 
photon are 


Sa | ds Hx)-} | ds i di PEA), Ey) vee (48) 


where H, and H, denote the terms involving the first and the second powers 
of « in (46). In order to evaluate the P-bracket in (48) we have to find the 
vacuum expectation values of P-brackets containing a pair of y,,’s or a pair of 
F,,,’s. It follows from (I, 23) and (I, 58) that 


{Plu %)s Viol X’)) do= = 3(0 Oo + S90 va)Dy AC x Je oS) (enle Sthe (49) 


where D,’(x—x’) is Feynman’s singular function for the Bravia ona! field. 
Also, we have for the radiation field 


p (2x2) 24se)\\ * 
( On; ° tax, Jif oe ia 


which gives 


(A, (x), A(x’) »o or 18 190 4 240(% =e :) 


Th 18 0 149 440(X aut: ds 


07D, (x—x 0? Dy (x — x’ 
(PR (2), Fag )o= Ba pager + Be 


Q 
_45 B°Dy(x— x’) _45 02D, (x — x’) 
fas ON 0X as Ox 0x) 
+18 9.49 240(% — X") +15 ,,39 45 940(X — x”) 
10970 pho pO — He) 80,00 40440 —H ). ae sees (50) 


*Jt should be remembered that for simplicity we are using the same symbols for the 
corresponding quantities in the Heisenberg and the interaction representations. 
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It can now be easily verified that the terms involving the 5-function in (50) give 
rise to a contribution (48) which exactly cancels the contribution of H. In fact 
it is well known that such a cancellation always occurs. Thus, instead of (48)% 
we may take 


See | dx { df Pi ah Hae be (51) 
provided we also take effectively 
02D y(x« — x’) 0? Dy(x« — x’) 
a) hee pirat 2 PR eet MB 
(PCF #), F(x )) »o= 2 (>, 0x,0x) Oia Ox, Ox, 
02 Dy (x — x’) 0? Dy (x —x 2). 
en Ox,,0x?, O55 Ox, Ox Ox ceccee (52) 


Substituting the value of H, from (46) in (51), we get 
Sproton= =e? | die | de! PE (yas) Fy) Fus®) — 4a #) Fol #) Ful) 
(Yaq(%”) Fag ') Fag’) — a7 pp(%') Fag(#’) Fag(#'))}. ++ +++ (53) 


Extracting the photon self-energy contribution from (53), we get, after some 
simplification, 


4 = ' 02D.(x— x’ Bh Cpe 
Spen=— ie [dv | de’ Fy) (SHE? + 1,, 2) 
=e oe age ava 


“DW — x LF ele «© 4 puma (54) 
We now pass over to the momentum-energy representation by putting 


F(x’) = F,,(p) exp (tpx’) 
Dy(x— x") = On aa dk(k? + A?) exp {ik(x — x’)} 


Dy(x—*x’) = a dk'(k’? + p?)-1 exp {tk'(x — x’)} 


where, following Feynman (1949), the photons and the gravitons have been 
assumed to have small masses A and p respectively, which have to be put equal 
to zero ultimately. Then, we get 


af? hike +BBuok? 
Shaion ral cael dx F(X) F,,(p) e foe [ dk (k2 +A2)[(R—p)? + p?] 
=—if de Fale) Pil) ° iP), < (| a een (55) 
k Rs ot 58 fe 
where J (Pp) = ny yi a | +k —pP ate] ee aOSG0 (56) 


Since the photon self-energy density Wyyoton is connected with Sypoton by the 
relation 


co 
S photon — — U a. dx WE acto, ’ 


eqn. (55) shows that W>noton involves only the electromagnetic field strengths. 
It follows that our result is gauge-invariant, and therefore the photon self-energy 
vanishes unambiguously. ‘This situation is much simpler than the corresponding 


situation in quantum electrodynamics, where the photon self-energy is found 
to be ambiguous. 
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Though we have seen that the self-energy of a free photon vanishes, it will 
be interesting to see what divergencies are involved in an internal photon 
self-energy ‘graph’. For this we have to evaluate the integral (56), which may 
be easily carried out by regularization and using familiar devices of integration. 


Thus we find 
2 fl 
Jul?) = — 7s { dx [38,,{2x€2 log 2—A log (wé2) + A log (2A) — A} 
0 
+X*( PP» + 20x”) {log (xé) —log (2A)}], «+ (57) 
where A=.x(p?+p7)+(1—«x)A?—px?, and the parameter €, which has the 


dimensions of mass, tends to Pepe We may write (57) as 


J (Pp) = — = Pings, — == cd Hye to") log finite terms, ~~ ..-.. (58) 


The quadratically divergent term in (58) may be absorbed by a renormalization 
of the electromagnetic field strengths (Gupta 1951), but the remaining 
logarithmically divergent term remains a source of difficulty. 


§7. SELF-ENERGY OF THE ELECTRON 


According to (I, 54), the Lagrangian density for the electron field interacting 
with the linear gravitational field is given by 


eat A Op _ Of gE 1 (O¥u» OYuy sar st] 
i= 5 (Prax aa gach id eye F raiags ar re 


— €(Yp»— B8urY ”4 (Frise ie aa se yt) + + (hy os _ se vt) f ‘ 


As in the case of the electromagnetic field, we may take the interaction energy 
density in the interaction representation as 


Hs) = tel Pon) { (Tee wee Lah) + (hese mee set) by +++ (60) 


provided we also take = iad 


et) P(SE, ie’) ) = Lett) POMS), He’) 


Oyp(se) Of(x’)\\ a 
ca Ca © ata) Bier, 
in addition to the nae result 
c(t, t')( P(A), Wx)) y= —FSA#—®). ees (61) 


It will also be convenient to put 


[e(t, 1°) PCM), 4(x")) dol, 


Vuy— SousY = Muy Bg hed (62) 
and we easily find that 
{PUT %)s 4%) Yo = 4(8,,29 po + Supra — SnvOae)Da(¥— 2%"). neces (63) 


Using (62), and remembering that h,,,=h,,,, we may write (60) in the simpler form 


H(x) =4«h,, (Frise = st vat) Ar ae SEe (64) 
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Now the second-order contribution of the S-matrix for the interaction of 


electrons and gravitatons is given by — 


seme | pa Af, | a dx’ P{H(«), H(x')} 
= = aes dx’ P { ols) (x) go outs) Yui (x ), 
Inde!) (Hots Bo He) ix yy ae (65) 


Extracting the electron self-energy a from Hie we get 


eisctron ar = ie * dx le a dx’ (9770, + CI oa SO Ap) D p(x s: x") 


= hi 0 ; 
x | —Hedynge Solvay, GO? + TED y, Sole — x yah) 
= one YpSp(*—*')y a ane) 


Weng ae Soe —w' rahe’) | 


Putting Wx!) =p) exp (ipx’) 

Dye — x)= ae { dk’ (k'? + p?)-! exp {ik'(x—x’)} 
= CE pty ie 
(27) yea | h2 + m2 
we find in the momentum energy representation 


Secouon= i {dv Ua )1(pWAp) explips) sees (67) 


Sy(x — x’) = ~~ exp {tk(x —x’)}, 


Ke 
where ip= “tet | dh: (3 SE Se Speen) 
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Using familiar devices of integration with regularization, we find 


1(P)= aga | | (Bim +-x)2p*—w(1 +9)*04( py) 


2 2 
x {log <5 +log “e — log at + {4um —(1+3x)(py)} 


x {300g (4)e— \e2 A (log 55 gee —log A+1) } ], 


where the parameter £, which is of the dimensions of mass, tends to infinity, 
and A =(p?+ p?)x + m?(1 — x) — xp?. 
Since for the free electron 


pPP+m=0, (“pytmyb(p)=0, > aaanes (70) 


we obtain from (69) 


Kem (3 4\ & Siaee BS 
I(p) = 128% Tamm (3!08 (3) 3 ~ 4 low ges - >), Se ras Wg a 


which is simply the gravitational self-mass of the electron. 


€ | 
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In order to find all the diverge ncies in an internal electron self-energy graph 
~ we have to evaluate (69) without using (70). Then it is found that J(p) is of the form 


I(p) =A + Blipy + m) + C(ipy +m)? + D(ipy +m)? +finite terms, ...... (72) 


where A and B are quadratically divergent and C and D are logarithmically 
divergent. ‘his situation is much worse than the corresponding situation in 
quantum electrodynamics, where the electron self-energy graph involves only 
two divergencies. 

According to our investigation of the ‘graphs’ for the gravitational 
self-energies of the photon and the electron, the gravitational field involves many 
more divergencies besides those occurring in quantum electrodynamics. This 
suggests that it is not possible to remove all divergencies from the gravitational 
field by mass and charge renormalizations. However, it must be noted that in 
our calculations we have not taken into account the contributions of the 
non-linear terms to the self-energy graphs. In order to see whether the 
renormalization theory really breaks down for the gravitational field, it would 
be necessary to carry out more involved calculations by including the non-linear 
terms. ‘This problem will be fully discussed in a subsequent paper. 
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ABSTRACT. A model of liquid helium is described, which incorporates the main 
features of the Landau (1941) and London (1939) theories, making use also of Mayer’s 
concept of clusters. The ordinary process of liquefaction is associated with the formation 
of clusters in coordinate space, while the A-transition is probably associated with the 
formation of small clusters in momentum space. ‘The latter process is probably prevented, 
or greatly modified, in helium 3. 

Possible reasons for the appearance of superfluid effects are given, and it is suggested 
that the unique properties of helium II may be traced to the fact that the interactions 
between two atoms is of about the strength needed to form a molecule. 
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§1. INTRODUCTION 

N a previous paper (1952), referred to as I, a comparison was made betweelll 
] the two most successful theories of boars helium II, usually known as the 
Landau (1941) and London (1939) theories. The following model | 
incorporates some of the most successful features of the two theories, and was 
put forward in I as a natural synthesis of them. It was supposed that one might | 
think of liquid helium II as a Debye lattice of helium atoms in statistical 
equilibrium with helium atoms in states of zero or small energy, and it was shown | 
that such states would be appreciably occupied only below a certain transition 
temperature whereas, above this transition temperature, any such nearly free 
atoms would be immediately incorporated into the Debye lattice. The condition 
determining this transition temperature was found to be 


(OS/ON)z~=0 or g=O ——— weweee (1) 


which, it was pointed out, is equivalent to the accepted criterion for Bose-Einstein © 
condensation in a ‘one-particle’ type of assembly, but, in this form, could — 
equally well be applied to a many-body assembly. 

A discussion of this model showed that it was, in principle, capable of 
accounting for many of the observed features of liquid helium II, such as the 
occurrence of a sharp transition temperature, the appearance of a thick film, 
the fact that the velocity of second sound at low temperatures is given correctly 
by Landau’s (1941) expression, the peculiar behaviour of the ‘viscosity of the 
normal fluid’ as a function of temperature and the successes of the phenomeno- 
logical ‘two-fluid’ theory in describing second sound and pseudo-heat conduction. 
There is, however, one experimental fact that shows clearly that the above model 
is still too specialized. Equation (1) requires that the chemical potential of 
liquid helium II should be zero or positive, whereas the experimental evidence _ 
indicates quite clearly that g is negative and of the order of 10-15 calories per mole. 
Although this figure is small compared with either the attractive or zero-point 
energies, it is large compared with NRT, and Ease ore cannot possibly be 
neglected. 

In this paper we shall show that a generalization of the model developed 
in I is possible, and we now make use of the concept of ‘clusters’ due to Mayer 
and Mayer (1940). We shall suggest that the behaviour of liquid helium II can 
be interpreted in terms of statistical equilibrium between ‘large’ and ‘small’ 
clusters in momentum space, while the ordinary equilibrium between liquid and 
vapour is interpreted on the lines of Mayer and Mayer’s original theory (1940), i.e. 
in terms of statistical equilibrium between ‘large’ and ‘small’ clusters in coordinate 
space. It will be suggested that the processes of clustering in momentum and 
coordinate space are of comparable importance in liquid helium II because the 
interaction energy between two helium atoms is approximately that required to 
forma bound state. If it were much larger than this, the assembly could probably 
be described in classical terms, the ‘clustering’ taking place in coordinate space, 
i.e. a net gain in energy results if atoms approach one another whereas, if it were 
much smaller than this, we should have a_ model resembling the perfect 
Bose-Einstein gas, the “clustering” then taking place in momentum space. 
A ‘cluster’ in momentum space is to be thought of as a number of atoms having 
nearly the same momentum, in which case (Mott 1949) a gain in energy can 
result in certain circumstances. Price (1951) has deduced a very similar result, 
that a weak interaction can result in a ‘ potential well’ in momentum space. 


———— 
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§2. MAYER’S CONCEPT OF A CLUSTER 
If we consider, for the present, a purely classical assembly, in which the 
interactions are supposed to be between pairs of atoms only and to be independent 
of their velocities, it is well known that the process of averaging over momentum 
space can be performed at once, and the evaluation of the equilibrium properties 
of the assembly could then be carried out completely, if we could evaluate the 
“configuration integral’ 


OT)=|...-fexm(-gpBe)due. ve (2) 


the integrations being taken over all possible values of the coordinates of all the 
atoms. Mayer and Mayer (1940) write this in the form 


Q(T)=J.... MCh)... dr {> Se eee: (3) 


where f,;=exp (—<«,/RT)—1 which is appreciable only if ¢,; is significantly 
different from zero, which, in turn, implies that atoms 7 and 7 must be fairly 
near one another. Mayer proceeds by expanding the product occurring in (3) 
and attempting to estimate the relative importance of terms of various types, 
such as 1, f,;, f;; f;, etc. after the integrations over the coordinates have been 
carried out. This product is expressed in terms of the ‘cluster integrals’ 5, 
defined as ; 


b= 77 - on [2 Wished 


in which =’ implies that the summation is to be restricted in the following ways 

(ESS 1 ea 

(b) Each of the terms f;, is to have at least one suffix in common with some 
other term. 

(c) Each of the / atoms is linked directly or indirectly with every other. 

It is then a simply matter to show that 

O(T) (Vb,)™ is 
NT eg RC AER Wolty as tise ‘aalreabe (4) 

where possible values of m, are restricted by the condition &/m,=N. 

By the familiar technique of searching for the value of m, which, for each /, 
gives the maximum term in the sum (4), we may replace (4) by the equivalent 
result 


log (om ?) gale DZ IN log Ze GaP Oe (5) 


where Z is determined by V%,/b,Z’=N. In these expressions, the method of 
constructing the cluster integral 5, implies that it is negligibly small unless every 
one of / molecules is near enough to at least one of the others for the interactions 
to be appreciable. Because of the nature of partition functions, the higher 
terms in (3) are a measure of the proportion of time a typical molecule spends 
in the neighbourhood of others, and the various terms in expression (4) measure 
the relative probabilities of the presence in the assembly of clusters of various 
sizes—a gas consisting mainly of single molecules (/=1)—while liquefaction 
occurs when the main contribution to (5) is given by /= N, which implies that 
nearly every molecule is physically near to several others. ‘Thus, all the 
equilibrium properties of the assembly are known if 6, is known as a function 
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of 7 and 7, though the relationship between ‘the cluster integrals and other 
accepted means of describing an assembly, e.g. by means of a molecular 
distribution function, may be very complicated. The ‘cluster’ concept is a 
very valuable one, easily understood intuitively, but hard to define in precise 
terms. (It could be made precise for an interaction energy that vanished beyond — 
a certain distance.) | 

In a quantum mechanical assembly, it is also possible to think of cluster 
formation in momentum space. Consider, for example, the well-known — 
expression for the partition function of a perfect gas of atoms obeying Bose 
statistics 


RT\3/2 2 8 
log f= V (=) leet hee ae (6) 
: 27mMkRT\3!2, 2 a8 ow 
with N= v(: Re ) oe Eas 


The resemblance of this expression to expressions such as (5) has already been © 
remarked upon by Kahn and Uhlenbeck (1938) and by Mayer and Mayer (1940). 
Examination of the method of deriving (6) shows that the sth term in the first 
series is related to configurations in which there are s atoms in the same quantum 
state in just the same way that the term 5,Z’' in (5) is related to clusters of size /. 
We find that, as long as the 5, are positive, the probable number of atoms in 
clusters of size J is proportional to /b,Z' while the probable number of groups 
of s atoms in the same quantum state is proportional to a°s~*?, 

The phenomenon of Bose-Einstein condensation of a perfect gas may be 
regarded as the simplest case of the formation of a cluster in momentum space, 
the occupation of the very lowest state, measured by the term «/(1—«) in (6), 
becoming appreciable below the transition temperature. Mott (1949) has called 
attention to another process that might accentuate such cluster formation. He — 
considers the gas of Bose atoms with small interactions, and, following a 
suggestion of Bijl, de Boer and Michels (1941), shows that the application of 
perturbation theory suggests that the formation of clusters in momentum space, 
i.e. the multiple occupation of one quantum state, is energetically favourable 
provided that the interaction between two atoms is not great enough to form a 
molecule. A similar point has very recently been made by Price (1951) but his — 
method fails to give any information if the interactions are greater than this. 
Very recently an attempt has been made (Matsubara 1951) to express the 
Slater ‘sum over states’ (for a quantum assembly with interactions between the 
atoms) in a form similar to (5). However, this can only be done at present by 
neglecting the lack of commutation between the momentum and the potential 
energy, which is equivalent to calculating the first term of a series in ascending 
powers of h, so that the validity of the results is doubtful. We shall not attempt 
a rigorous discussion on these lines, but we first assume the possibility of a 
development such as (5) or (6) in a quantum assembly, a result made plausible 
by Kahn and Uhlenbeck’s work (1938), and then make certain assumptions about 
the behaviour of the 6,’s in liquid helium. Some of our considerations will be 
based on the fact that the interaction between two helium atoms is certainly of 
the order of magnitude required to form a molecule and we need not expect 
them to apply to other liquids. In particular, we shall introduce a distinction 
between ‘small’ and ‘large’ clusters which can be made only if the interaction 
is near this critical value. 
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lf the interaction is large compared with the critical value, cluster formation 
takes place mainly in coordinate space. In this case, the atoms are already 
prevented by their repulsive fields from approaching one another too closely 
and the influence of the type of statistics obeyed by the atoms should be quite 
small (Wergeland 1946). If the interaction is smaller than the critical value, 
two atoms do not gain any energy by approaching one another (as they would 
if they behaved classically) because the gain in potential energy is more than 
offset by the increase in ‘zero-point’ kinetic energy. On the other hand, we have 
already noted Mott’s (1949) result that, in just this situation, it is energetically 
profitable to place a number of Bose atoms in one quantum state, and this 
‘clustering in momentum space’ would be prevented, or would take place in an 
entirely different way, if the atoms obeyed Fermi—Dirac statistics. 


§3. THE CONCEPT OF ‘LARGE’ AND ‘SMALL’ CLUSTERS 

Let us assume provisionally that the binding between two He, atoms is 
sufficient to form a stable molecule. We shall then have a system very like the 
deuteron, namely a molecule with small binding energy and only a single bound 
level. Similarly, if the interaction is just less than the critical value we can gain 
only a small amount of energy by assigning two atoms to the same momentum 
state. As we add more and more atoms to a cluster, there will evidently come 
a stage where the cluster can have excited levels, and a very large cluster can 
evidently go into a large number of excited levels without breaking up. We shall 
describe a cluster as ‘small’ if it has no excited levels and ‘large’ otherwise. 
We shall suppose further that the spectrum of levels of a large cluster is, to a 
first approximation, the same as that of a Debye lattice, and we shall consider the 
conditions for equilibrium between the two types of cluster. We shall then 
show how this first approximation can be improved. ‘The energy of the assembly 
is made up of various contributions: (a) the translational and rotational energies 
of each cluster; (6) the energy gained when one forms free atoms into clusters; 
(c) the energy associated with the internal degrees of freedom of the clusters; 
(d) the intraction energy between clusters. 

For a ‘small’ cluster, contribution (c) is absent because there are no excited 
levels. Contribution (d) can probably always be neglected. In a classical 
assembly Mayer and Mayer’s theory (1940) shows that PV =kTX,m,, where 
m,is the most probable number of clusters of size /, that is the assembly behaves 
formally as if it were made up of &,m, non-interacting particles, so that the 
- interactions between atoms are completely taken into account by the introduction 
of clusters without introducing an inter-cluster interaction. It does not follow 
that a similar result holds strictly for quantum assemblies (for which PV is 
not equal to NRT even in the absence of interactions, except in the limit of very 
high temperatures), but it is probably a reasonable first approximation. It holds 
for the perfect Bose-Einstein gas as well as for a ‘classical’ gas. 


§4, EQUILIBRIUM BETWEEN ‘SMALL’ AND ‘LARGE’ CLUSTERS 

We consider first the condition of statistical equilibrium between ‘large’ 
and ‘small’ clusters when both have zero or small velocities of translation, 
secondly, the conditions in which they can transfer momentum to one another. 
For this purpose, it is probably sufficiently accurate to think of p ‘small’ clusters 
each containing / atoms, where / is the most probable size of a small cluster, in 
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equilibrium with a large cluster containing N—N, atoms, where N,=/. 
(We measure interaction energies reckoning that of atoms a great distance apart _ 
as zero.) The total free energy of the small clusters is simply N.£, where —&, 
is the mean binding energy per atom in such clusters. The entropy is negligible — 
since such clusters have no excited states. The free energy of the large cluster is 


(N—N,) | Yo+ 2R0— © (eT) (5) - ‘2 “i pa (7) 


where Vy is the mean potential energy per atom of the interaction forces, 9RO/8 
is the zero-point energy per atom of the large cluster, treated as a Debye lattice, 
and the form in 7“ takes account of the entropy and extra energy associated with 
thermal agitation. The condition for minimum free energy of the assembly is 


aF/aN,=0 or 
dE, 9 ee ines ag Ws 
eT ae sho+ TRT(Z) 36. oe (8) 


if we neglect any changes in V,) or © due to adding another atom to the large 
cluster. (If the addition took place with the large cluster keeping constant 
volume, the terms involving 0@/0N, could not be neglected, but we are keeping 
the volume of the whole assembly constant, so the neglect is probably justified.) 

The behaviour of this model evidently depends on that of £,as a function of J, 
the mean size of the ‘small’ clusters, and / in turn depends on N,, the total 
number of atoms in the ‘small’ clusters. As an illustration, consider the result 
of assuming the total energy of the ‘small’ clusters to be given by 

NE NEG ANG ee eer (9) 

Assumptions of this kind are often used in discussing co-operative phenomena. 
Moreover, it is known (see for example Price 1951) that quadratic terms of this 
nature appear when one treats the quantum mechanical many-body problem by 
perturbation theory. 

Equation (8) now gives 

; ot T\3 
2AN,=F,)-£,- Fk (5) vo Wee od Peer eto (10) 
5 1) 

where we have written E, for V)+ 9R@/8, the energy per atom in a large cluster 
at absolute zero. The condition for stability requires A to be positive. If 
E, —E, is positive, this leads to N, taking a finite value, not necessarily as great 
as N, at absolute zero, then gradually diminishing to zero at some transition 
temperature. Behaviour similar to this was expected in I for N,, the number - 
of atoms in ‘3-states’. Substitution of expression (10) for N, into (9) leads 
to an extra term in the specific heat proportional to J’, which disappears abruptly 
at the transition temperature at which N, vanishes. Thus, this first approximation 
gives a result of the right form, though it is not possible to choose A and E,—E, 
to fit both the observed transition temperature and the observed height of the 
specific heat anomaly (the latter coming out about 20 times too small if we 
choose £, — Ey to give a transition temperature of 2:2°K). It has, however, been 
pointed out in I that, according to Blackman’s work (1935), strict concurrence 
with the Debye law is not to be expected even for a perfect crystal lattice above 
temperatures of the order of 0-010 (~0-3° for liquid helium). If we admit a 
deviation from the Debye law of the type discussed in I, as suggested by 
Blackman’s work (1935), together with some co-operative effect such as that 
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involving the smaller clusters in (9) and (10) (due to #, being a function of N,) 
it is easily possible to explain the observed specific heat anomaly quantitatively, 
but we have too many adjustable constants for this fact to be significant by itself. 
The increase in the number of clusters in the liquid phase as the temperature 
falls could, by analogy with Mayer and Mayer’s theory (1940), imply an increase 
in pressure at constant volume, that is a negative expansion coefficient, as is 
observed. For a classical assembly Mayer found that the pressure was strictly 
proportional to the mean number of ‘clusters’ present. 

So far, we have left open the question whether the small clusters are to be 
thought of as occurring in coordinate or in momentum space. Either would 
be consistent with £, for the small clusters being a fairly rapidly varying function 
of /, as the data seem to require, and either would be consistent with the 
interaction between two atoms being near the critical value. However, the 
following considerations seem to point to the clustering being in momentum 
space, though, on present knowledge, a third possibility cannot be excluded, 
viz. clustering partly in coordinate, partly in momentum space, that is in phase 
space as a whole. 

(a) The ordinary process of condensation of the vapour involves the sudden 
appearance of large clusters in coordinate space. It seems unreasonable that 
the large clusters should begin to break up into smaller ones again as the 
temperature drops still lower, which is what would have to happen if our 
“small clusters’ were in coordinate space. 

(b) The ‘small’ and ‘large’ clusters are supposed to be simultaneously 
present in the liquid. If the clustering were in coordinate space, opalescence 
would be expected. (In the ordinary liquid—vapour equilibrium, the large 
clusters form the liquid, the small ones the vapour.) 

(c) Clustering certainly occurs in momentum space for the perfect Bose gas. 
It is reasonable, therefore, to expect it in a liquid with weak interactions. 


§5. COLLISIONS OF “SMALL” AND “LARGE” CLUSTERS 

We are thus led to think of helium II as consisting of one or more ‘large’ 
clusters, resembling a Debye lattice, in equilibrium with ‘small’ clusters. On 
this basis, we can define three distinct ‘masses per unit volume’ and three distinct 
“stream velocities’, as follows: (i) p, and v, the mass and velocity of propagation 
of the centre of gravity of the Debye lattice or ‘large’ clusters. (1) p,, and vy 
the mass and stream velocity associated with the ‘phonons’ which describe 
the excited states of the Debye lattice. 

In the Appendix to I, it is shown that, given any assembly of phonons, one 
can determine, first the total momentum, p,), v), per unit volume, that they 
could transfer to a material body, secondly the apparent stream velocity 7, 
this being the velocity of a material body in statistical equilibrium with these 
phonons, hence the effective phonon density p,,, which is found to be identical 
with Landau’s expression (1941). (iii) p, and v, the mass per unit volume and 
mean stream velocity of the small clusters. In I it is pointed out that free 
particles of a mass comparable with that of helium atoms do not behave 
classically at very low temperatures, their lowest wave functions being extremely 
sensitive to fields such as that of gravity. ‘This could probably be allowed for by 
the usual device of introducing effective masses for these clusters, differing from 
their gravitational masses, but we shall not introduce this refinement in the 
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quantitative treatment. Effective mass in this context is defined not by Landau’s — 
method but in the normal quantum mechanical way as (total force)/(rate of | 
change of flux) (Seitz 1940). | 

We can find a number of relations between these six quantities and the four 
quantities p,, Up) ps V, used in the phenomenological two-fluid theory. (@) ppn=Pn— 
because the observed velocity of second sound at very low temperatures requires 
such an assumption, and also because the phonons are the means by which a 
Debye lattice exchanges energy with its environment, for example, a moving 
boundary such as an oscillating disc. (6) po, +p.%.=PrUn +ps¥; because sound 
waves, in both classical and quantum mechanics, never result in any transfer 
of matter as long as anharmonic terms are negligible, but v, and y% both — 
correspond to actual movement of matter. (c) p,+p;=pn+ps=p because either 
of these sums represents the total mass of helium atoms present per unit volume. 
(d) v, =p, because it can be shown (Dingle 1951) that the entropy of a Debye 
lattice is transported with the velocity v,,. These equations are sufficient to 
determine all four quantities p,, Up, ps, Vs. Penrose (1951) and Dingle (1952) 
write Up,=V,—v; instead of the last of the above equations, so we shall also 
examine the physical consequences of this alternative assumption. 

On the assumption that v,=v,, we derive without difficulty the equation 


yj CRO Con) Fela laa eee (11) 
P ~Pph 

(The absence of ‘superfluid’ effects means v,=v, not v,=0.) The two main 
facts that have to be explained by any theory of liquid helium are: (a) that 
relative flow of superfluid and normal fluid can be initiated and maintained by 
extremely smal] pressure or temperature differences; (5) that, once initiated, 
such relative flow is damped out only very slowly (some recent unpublished work 
by D. V. Osborne suggests that the characteristic time is of the order of seconds). 
For large values of the relative velocity, stronger dissipative effects do occur. 

We are suggesting that both of these facts can be accounted for if we 
postulate the presence of ‘small’ clusters in the liquid. Fact (a) may be a 
consequence of their having small effective masses. Equation (11) shows that, 
in an ordinary liquid, the absence of any ‘superfluid’ effects seems to require 
that v,=v,, even if p,,>p. ‘This is actually the condition for thermodynamic 
equilibrium in a moving Debye lattice (see Appendix to I). The mechanism 
that maintains this equilibrium is presumably the Umklapp process of Peierls. 
(1929), which plays an important part in determining the heat conductivity of a 
crystal, and which may be described roughly as a Bragg reflection of phonons 
by the lattice, the net result of which is a transfer of momentum between the 
phonons and the centre of gravity of the lattice. We shall assume that this. 
process is effective in maintaining the equality of v, and v,, even in helium II. 
We now have to consider what mechanisms are available for maintaining 
equilibrium between the ‘small’ clusters and the Debye lattice that is, equality 
of v, and v. The Appendix to I shows that a given assembly of phonons cannot 
be in true statistical equilibrium with a material body moving at any velocity 
other than wv, but it does not follow a priori that any effective mechanism for: 
maintaining equality between v, and v,, actually exists. We consider a collision 
between a ‘small’ cluster and the Debye lattice at rest. (Collisions between 
the small clusters themselves lead to no change in p,v, and can be ignored.) 


Sf 
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A collision between a small cluster and the Debye lattice may result either in the 
transfer of momentum to the centre of gravity (which, by eqn. (11), does not 
affect v, — v, since it leaves p,v, + p,v, constant) or in the creation or disappearance 
of a phonon or phonons. The latter are the only processes that can result in a 
change of v,, —v,, and we may consider two possibilities : 

(a) The lattice is knocked into a higher excited state by the collision. Since 
the lattice is so much heavier than the cluster, this process becomes possible. 
if the energy of a phonon is less than the kinetic energy of the small cluster, the 
situation being analogous to that of an electron colliding with an atom, losing 
some or all of its kinetic energy, and leaving the atom in an excited state. Now, 
according to the Debye theory, the sound quantum of lowest energy would have 
a wavelength of about lcm in an ordinary vessel, and the corresponding 
frequency would thus be v=2:5 x 10sec"! so that the critical velocity of a ‘small” 
cluster of / atoms would be given by equating its kinetic energy to hy, the energy 
of a phonon: $/mv?=2-5 x 10*h or v=7/I'2 cm/sec for m equal to the mass of a 
helium atom. ‘This calculated velocity is somewhat smaller than the observed 

critical velocities, which again suggests that the effective mass of the clusters 
may be less than their actual mass. 

(6) A process involving a large number of phonons but a resultant very 
small change of energy occurs. As an example, one may suppose that the small 
cluster initially collides with a few atoms in the lattice from which an elastic 
wave spreads out. In I it is pointed out that a process involving the initial 
displacement of only a few atoms in the lattice necessarily involves a large 
number of phonons, and there is therefore no lower limit on the energy required 
to bring about such a process. 

While we certainly cannot say a priori that processes such as (6) are 
ineffective, it does seem suggestive that the assumption of the predominance of 
processes (a) does predict critical velocities of the correct order of magnitude. 
If we accept this idea, it is also possible, because of the variation of the expected. 
critical velocity with size of cluster, to understand why the observed critical 
velocity is rather indefinite, at least in the bulk liquid, ranging, apparently, from 
values of the order of 1cm/sec in an oscillating disc type of experiment, to 
100 cm/sec for ‘second sound’, so that many of the facts can be accounted for by 
assuming instead a frictional term that increases rapidly with the velocity 
v,—v; (Gorter and Mellink 1949). 

We now discuss briefly the consequences of the alternative assumption that 
Uph =Un—Vz- If we assume this, we have some difficulty in accounting for the 
non-existence of superfluid effects in other liquids. v,, zero means that on an 
average every phonon is balanced by another one of equal and opposite 
momentum, which implies some fairly strong coupling between such phonons. 
It is hard to see why such coupling should be absent in the one liquid helium IT, 
yet should appear suddenly when the temperature is raised. 


§6. CONCLUSIONS 

We have arrived at the following tentative picture of liquid helium: 

(a) The condensation of vapour to form liquid helium I is probably quite 
normal, and should be describable by the Mayer and Mayer theory (1940) except 
that we must take account of the zero-point energy, which, in both ‘small’ and 
‘large’ clusters, is comparable with the attractive energy between atoms. 

41-2 
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(b) The A-transition is interpreted as heralding the appearance of ‘small’ 


clusters in the liquid itself (probably in momentum space). It is then necessary 


to assume that the binding energy, per atom in the average ‘small’ cluster, is_ 
numerically greater than that in a ‘large’ one, and that it varies fairly rapidly — 
with the size of cluster. All this is quite a plausible consequence of the known _ 


fact that the interaction between two helium atoms is of about the strength 


needed to form a molecule. ; 

(c) The appearance of ‘superfluid’ effects is interpreted as a consequence 
of the fact that the ‘small’ clusters cannot easily lose their kinetic energy by 
collision with a ‘large’ one. It is also possible that their quantum mechanical 


effective mass is less than their actual mass. 

(d) The possibility of the formation of clusters in momentum space as well 
as in coordinate space is probably confined to helium atoms because of the 
particular strength of the atomic interaction. Clustering in coordinate space 
is not much affected by the symmetry properties but clustering in momentum 
space will be quite different for the two helium isotopes (if indeed it occurs for 
He, at all). This seems in agreement with the facts known about He. 
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ABSTRACT. Gaseous GaF and InF are formed when the metals are heated with 
aluminium trifluoride, and sufficient pressures of the monofluorides for observation of their 
absorption spectra are obtained at about 1000° c for GaF and 1400° c for InF. The three 
ultra-violet band-systems, 4, *IU,, ?11):<-1X+ expected by analogy with GaCl and InCl have 
been photographed on spectrographs of moderate dispersion, and vibrational analyses have 
been carried out. The 1II states converge rapidly to give rather sharp dissociation limits : 
the products are probably Ga or In, (?P3;2), + F(#P1/2). The bearing of the present results 
upon the interpretation of the electronic levels of TIF and on the dissociation energy of AIF 
is briefly considered. 


$1. INTRODUCTION 


UR interest in the spectra of GaF and InF arose in the first place because 
C) of the possibility of determining their energies of dissociation. It was 
thought that values for these molecules would probably throw some 
light upon the very large discrepancy which existed between the results of the 
spectroscopic and thermochemical determinations of the energy of dissociation 
of AIF (see §5). Moreover, if sufficiently precise values were derived, they could 
be combined with results to be obtained from future equilibrium measurements 
to provide values for the dissociation energy of fluorine. GaF and InF were also 
of some interest as, apart from BI, they were the only two molecules of the group 
stretching from BF to TII about which there was no spectroscopic information. 
The possibility of studying these molecules by a simple absorption technique 
was suggested by the work of Gross, Campbell, Kent and Levi (1948) on 
heterogenous equilibria of the type 2Al+ AIC],;=3AICl. ‘Their results showed 
that it would be reasonable to attempt to produce small pressures of the gaseous 
monofluorides of gallium and indium by heating the metals with aluminium 
trifluoride. ‘This expectation was fulfilled, and sufficient GaF and InF for the 
examination of their absorption spectra were found to be readily produced 
in this way. 

Three ultra-violet band-systems have been found for each molecule, which 
correspond closely with those already known for the other gallium and indium 
halides (Miescher and Wehrli 1934, Wehrli and Miescher 1934). ‘The forbidden 
systems, A®II)+<-x1X* and B°II,<-x'2*, are situated not far from 3 000A, and 
appear with only low intensity: the other systems, c1I[<-x!X*, lie at shorter 
wavelengths, and seem to be much stronger. ‘The vibrational levels in the states 
c1II converge rapidly, and fairly precise dissociation limits have been obtained. 
Their significance is discussed in a later section. 
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§2. EXPERIMENTAL 


Mixtures of Ga or In metal with AIF, were placed in Pythagoras boats in _ 


aluminium oxide tubes. The tubes carried water-cooled brass end-pieces with 


silica windows and were heated in a Globar furnace: the effective column length ~ 


was about 20cm. The GaF plates were taken at about 1000°c, the InF at 
about 1400°c. 

A hydrogen lamp was used as source of continuum. The absorption spectra 
were photographed with Hilger Small and Medium quartz prism instruments 
on Ilford Q1 and Process plates. Cu and Fe arc comparison spectra were used : 
the wavelengths of the lines were taken from the M.1.T. list (1939). 


§3. ANALYSES OF THE SPECTRA 
(i) Appearance 


‘ The reproductions in figs. 1 and 2 give a general impression of the absorption 
spectra. Both molecules show two distinct regions of banded absorption. In 
GaF, there are weak violet-degraded bands of two kinds, single- and double- 
headed, in the region 2 900-3 1004, and a much stronger system of predominantly 
red-degraded double-headed bands lying between 2050 and 22004. ‘The systems 
of InF are similar, but shifted slightly to longer wavelengths, at 3050-3 300A 
and 2275-2475 A respectively. 

Apart from some GaCl bands, which appeared on one or two of the GaF 
plates, no other molecular absorption was observed; in particular the plates 
were quite free from AIF. Lines of Ini appeared strongly, and on the exposures 
at the highest temperatures, there was strong development of the series 77S, )j, 
Dojo, 5/2 — 9° P19, 32 Of Int with limits at 2145 and 2252A (Garton 1951). 


(ii) Vibrational Analyses 


The vibrational analyses presented no especial difficulties. It is most probable 
that the far ultra-violet transitions are !I] —1X, as has been demonstrated for 
the corresponding system of InCl by the rotational analysis of Froslie and 
Winans (1947). ‘The two heads are thus R and Q, or in the few violet-degraded 
bands, P and Q. In all the systems, B’=B”, and the 0,0 bands are much the 
strongest. ‘The 0,0 band of the c—x system of InF is of unusual type, with two 
Q heads, an R head and a rather diffuse P head. This kind of band has been 
encountered in other systems in this group of molecules, for example in BBr 
(Rosenthaler 1940): in the present instance its formation arises from the fact 
that whereas F’'(J)>F’(J) for small J, F’(J)<F"(J) for larger J, i.e. By’ > By’, 
but Dy >D,’. 

Much less information was obtained about the longer wavelength systems, 
but a sufhcient number of bands was observed to confirm the ground-state 
vibrational intervals deduced from the analyses of the c—x systems. ‘The weakness 
of the systems A, B—x and the observation that the a—x bands possess only one head, 
whereas the B-x bands have two, support their identification as a?1I9+, B31], —1+ 
as in other molecules of this group (see, for example, Wehrli and Miescher 1934). 

The vibrational analyses are given in tables 1-6: the measurements given 
for the B, C-x systems are for the Q heads, where resolved; those for the 
A~X systems are for the P heads. 
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Table 1. Deslandres Scheme for the System a*II)+—x!2+ of GaF 


P heads. 
4 33858-6 (0) 
2952-60 
3 338311 (1) 
2955-00 
653°8 
2 33800°8 (2)  % 331953 (1) 
2957-65 8 3011-60 
642°7 
1] 33767-7(2) ‘& — 33158-1 (3) 
2960-55 8 3014-98 
651°7 


0} 33116-0 (5) 
‘ 3018-81 


Table 2. Deslandres Scheme for the System B?II, —x12X+ of GaF 


Q heads, except where the measurements are in parentheses, indicating unresolved 


heads, Q+P. 
5 ; (34220-2) (1) 
2921-40 
4 (34190-5) (2) 
2923-94 
3 (34160-2) (3) 
2926-53 
(636) 
2 34137-2 (2) %- 33528-4(2) & (32915-8) (0) 
2928-50 9 2981-68 8 3037-17 
647°2 647°2 
1| 34107-1(3) = 33490-0 (5) cs 328812 (0) 
2931-09 ~ 2985:10 2 3040-37 
\o 
658-2 655°2 
0] 33448-9(10) 5 32834-8 (0) 
| 2988-77 3044-67 
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Table 4. Deslandres Scheme for the System a®II)+—x!X+ of InF 


P heads. 
4 31106:8 (0) 
3213-81 
3 31081°6 (1) 
3216-41 
550°4 
2 31054-1 (1) ot 30531-2 (1) 
3219-26 a 3274-40 
558°4 
1| 31023-5 (1) oe 30495-7 (1) 
3222-44 ¢ 3278-21 
563°9 


0| 30459-6 (2) 
+ 3282-10 


Table 5. Deslandres Scheme for the System B*II,—x1X+ of InF 


Q heads, except where the measurements are in parentheses, indicating unresolved 


heads, P+Q. 
5 (31940-4) (0). 
3129-93 
4 31925-6 (0) 
3131-38 
3 31900-0 (2) 
3133-89 
552°8 
2 31871-2(3) 2 31347:2(1) >  (30816-0) (1) 
3136-72, S 3189-16 th 3244413 
558°5 (560) 
1] 31841-9(2) % 31312-7(5) © (30783-9) 2(1) 
3139-61 @ 3192-67 i 3247-52 
565°5 (571) 
0] 31276-4(6) % (30738-6) (0) 
4| 3796-38 = $3 3252-30 
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(iti) Rotational Constants 


An attempt may be made to interpret the R-Q and P-Q head separations 
in these systems as follows. We first estimate the value of 7,”. Let us 
consider AIF. Now the increase in internuclear distance on going from the 
diatomic molecule BX to AIX is almost constant: we find (from Herzberg 1950) 
AIH-BH, 0-413A; AICI-BCl, 0-422A4; AIBr—BBr, 0-4084; AlO-BO, 0-4134. 
The mean difference is 0-414A, which, with ,7,"(BF)=1-262A, gives 
7,’(AIF)=1-684. For InF, InH-BH, 0-6054; InCI-BCl, 0-6024, giving 
7. (InF)=1-87A. No internuclear distances have been measured for diatomic 
gallium compounds: we should guess, however, that r,"(GaF) lies between the 
values for AIF and for InF, and that it probably lies nearer to the value for AIF 
(see, for example, Sidgwick 1950, p. xxix). We take r,"(GaF)=1-75A. Then 
B. (GaF) =0-369; B,"(InF) =0-295, cm-1, 


Table 7. Band-head Separations, Q—P, in the System B—x of GaF 
Band (v’, v”) 1,0 Pi 0,0 (li DD) 0,1 i 
Av(cm~)obs. 15:2 10-7 11-0 10-2 8:1 11-1 10-2 

calc. 14-9 11:1 14-0 10-6 8°5 10-2 8-2 

Table 8. Band-head Separations, R—Q or Q—P, in the System c—x of GaF 

Band (v’, v”) 2,0 10 Dal Bee 0,0 1,1 BD Re) 
Degradation R R R R R R R R 
Av(cm~!)obs. 6°8 15:8 8-4 3-6 49-1 28-6 14-7 6:2 
calc. 6°8 12-4 8-4 6:3 53-6 18-5 10-9 7-6 
Band (v’, uv”) Ot {2 2,3 3,4 Or?) 1e3 2,4 355) 
Degradation Vv R R R V V Re R 

Av(cm—!)obs. large 40-7 22°3 9-1 50:8 large 31-4 Sei! 

calc. large 3:78 P57 9-7 pS.) large 28-6 (S325) 
‘Table 9. Band-head Separations, Q— P, in the System B—x of InF 
Band (v’, v”) 1,0 pest 32 4,3 0,0 Lat PP) 
Av(cm~)obs. 8:5 7:2 Tee 8-5 5-9 6:5 6-9 
calc. Pee 7:8 8-1 8:3 6:2 6:3 6:4 


With these assumed values of B,”, values of «” were derived from the 
head—head separations, Av, in the v’-constant progressions of the C—x systems, 
using Av=(B’ + B”)?/4(B’—B”). Values calculated from B,” = B,” —«"(v" + $) 
were then used to obtain B’,, and hence B,’ and «’ for the systems c—x and B-x. 
The values of the rotational constants derived in this way are given in table 11, 

The observed and calculated band-head separations are compared in 
tables 7-10. It will be seen that the calculated values reproduce quite well the 
general trends observed in the values of the separations with v’, v". 

It should be emphasized that this agreement means no more than that a 
self-consistant set of rotational constants has been obtained on the basis of the 
assumed values of B,’. The values of B,” can be varied in fact within quite 
wide limits to give other sets of rotational constants which in turn give reasonably 
satisfactory separations. If there were a convenient criterion of best fit, it would 
of course be possible to determine B,” by a trial-and-error method. ‘This can hardly 
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be expected to be a fruitful procedure, however, since for the small separations, 
the error in the experimental values is proportionately large, while for the large 
separations, of which there are good experimental values, the calculated values— 
because of the importance of D-terms—are no longer well represented by 


Av=(B’ + B")?/4(B’— B”). In the case of GaF, the effect is still further 
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complicated by small isotope displacements. 


Table 10. Band-head Separations, R—Q or Q— P, inthe System of c — x of InF 


, Band (v’, v”) Degrad. obs. 


Si 


GaF 


= 


24-7, 10-*4g. 


InF 


f= 


27-0,  10-%4g. 


Notes : 


R 


ARAM ADAAAAKS DAD AAR 


Av(cm~) 
calc. 
4-0 2 
QE 8:6 
6:9 5-9 
4-2 4-4 
24-7 26:0 
YD ven 
5-0 Sal 
~35 28:2 
SHED large 
20°5 16:8 
10°5 9-0 
6°5 6:0 
4-0 4-4 
Ds 365 
18-2 15:8 
large large 
29-8 Sy 
14-2 AEP 


Av(cm~?) 
Band (v’, v”) Degrad. obs. calc. 
4,5 R 7:2 73 
5,6 R 4-9 5 
6,7 R 25 5-9) 
02 V 10:8 10-9 
13 Vv ~20 21-4 
2,4 R ?42°7 large 
3.5 R 20-8 19-2 
4,6 R 10-4 9-4 
sys R 5:8 6-1 
6,8 R 2°9 44. 7 
DS Vv 16:3 34-9 
3,6 R 30-2 45°8 
4,7 R 14-0 13-7 
5,8 R 6°3 TOS 
6,9 R 5-4 5-2 
Bil ? large large 
4,8 R 18-2 22-7 
5,9 R 9-7 9-9 


Table 11. Summary of Constants for GaF and InF 


State 


ell 
Bal 


x tat 


col 
Bay 


Ve 


47362:3 


We 


542°3 


33430°5: 661-, 
A ®IT)+33096-, 


0 


663°. 
O23e2 


42807:5 463-9 


31256°8 


A IIo, 30441-, 


ee Sain 


0 


573-0 
S72" 
534-7 


Kee 


7-75 
3°8 


VeWe Be a r(A) Do(ev) 


—0°9, [0:368,] [0-008,] [1:75] 0:4, 
—  [0:377,] [0-0006] [1-73] = 


— [0-369] [0:003.] [1-75] eis 


—0-5* [0-300¢] [0-006,] [1:85] 0-4, 
—  [0-310,] [0-002,] [1:83] — 


— [0:2955)) f0002;) [STIRS 45 


on 


So RE Apr ee eit a ~ leew 


Ke(dyn/ene 
x 10-8) 
2:58 
3-85 
3°87 
3-41 


2°06: 
3°15) 
3:14 — 
2°74 


The rotational constants are of relative, rather than absolute significance : see text. The 
least certain value is that of a for B*II, of GaF : most probably the value given here is too 


small. 


* There is some evidence that levels with v’>7 in cI of InF begin to converge 
rather more rapidly than is given by this expression. 
The constants have their usual significance : thus for InF, c—x, 


where u=v-+3. 


v(Q head) =42807°5 + [463 -9u’ —7-3 ,u’®—0-5u’3] — [534-7u” —2-56u”?], 
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In spite of these limitations on the absolute accuracy of the rotational 
constants, it is more than probable that the values are correct relatively, i.e. that 
they show the correct trends from one electronic state to another of the same 
molecule. In particular, it may be concluded that (i) the equilibrium internuclear 
distances vary much less among the states of these molecules than do the force 
constants, (ii) the states c are anharmonic to a high degree, as shown both by the 
large values of x,w,, Vow, and of «,. 


1 §4. ENERGIES OF DISSOCIATION 


The vibrational levels in the states c'II converge rapidly, both for Gal’ and 
for InF. Numerical extrapolation, using the vibrational constants given in 
Bable 511, gives. Dy (c: GaF)=3545 cm; D,'(c: InF)=3 465 em. ~The 
dissociation limits above v”=0 are 50870 and 46240cm™ for GaF and InF 
respectively. 

It remains to identify the atomic products of dissociation of the states c. 
When the molecules of this group are considered as a whole, it proves difficult 
to draw a consistent picture of their dissociation energies except on the 
assumption that the process is MX(c14I)=M(?P3.)+X(?P,)). The arguments 
may be best illustrated by reference to the molecules InCl and InBr. 

For InCl, an upper limit for Dy” was found by Wenk: (1941) from a 
study of atomic fluorescence. ‘The effective wavelength for the process 
InCl—In(?5,/.)+ Cl(?Ps.) was found to be 16404, ie. 60980cm™+4. The 
interval In (?S,).,—?Pj)s) is 24373 cm“ so that Do” <36610cm™. Now v9 9 for 
the c-x system of InCl is at 37415cm™!: a very short extrapolation of the 
vibrational levels in c gives a limit at 37950cm™. The excitation energy of 
Cl(??P,.) is 881cm™ so that if the products at this dissociation limits were 
In(?P/2) + Cl(??Pj)2), Do” =37 950 — 880 = 37 070 cm, which is higher than Wenk’s 
upper limit. Thus one of the products of dissociation of state C must be 
In(??P3).), with an excitation energy of 2210cm™, and D,” is at least as low as 
37950—2210=35740cm™. It will appear that, with great probability, 
InCl(C)—+ In(?P 3/2) + Cl?P;;2), so that Do” 35 740 — 880 = 34 860 cm?. 

In the case of InBr, there is no indication that the state C possesses a potential 
minimum, the c—x system appearing as a set of fluctuation bands. Analysis of 
these absorption maxima leads to a dissociation limit at 33 940 cm above v” =0. 
The atomic fluorescence of InBr has been studied by Sedov and Filippov (1934) 
and by Wenk (1941). ‘The former give 1960+ 304 as the long-wavelength limit 
for production of the lines In(?S,/.2—?Psj 4j2): the latter finds that the most 
effective wavelength is 1830A. With the last figure, Dy” <30270cm™!. Now it 
is reasonable to assume that the products of dissociation of the states c in InCl 
and InBr are in the same electronic states. It has been shown that the 
In atom is éxcited, *P3,, from state c in InCl. For InBr, we have 
Dy” =33 940 —2 210 =31730cm™, if InBr(c)—In(?P3.) + Br(?P3j2). This figure 
is, however, compatible neither with the upper limit derived from Wenk’s 
work, nor with D’~26600cm™! from Sedov and Filippov’s observations. 
It must therefore be concluded that Dy” =33 940 — (2210 +3 685) =28 040 cm“ 
corresponding to the process InBr(c)—In(?P3).) + Br?P}j). 

It is likely that the other c'II states in this group of molecules behave 
similarly, and indeed there is a good deal of evidence for this view, which 
however need not be considered here. For GaF and Inf, assuming that the 
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states c dissociate into M(?P3.)+F(?Pij2), we find Do’ =6:15 and 5-41 ev | 
respectively. Ai 

It appears then that the two major problems 1 in determining the dissociation _ 
energies of these two molecules, (i) the precise measurement of a dissociation limit, © 
and (ii) the identification of the atomic products at that limit, have been solved” 
reasonably satisfactorily. It should however be noted that there is a possibility 
that the convergence limits correspond, not to dissociation limits at infinite” 
separation of the atoms, but rather to potential maxima. According to our 
present information such a situation appears to arise in the case of AIF (see 
below, §5). ‘Thus it is possible that the true values of Dy” for GaF and InF are 
somewhat lower than those given in table 11. It is hoped to elucidate this point, — 
and to check the identification of the atomic products of dissociation of the 1II 
states, by thermochemical work on these molecules now in progress. 

Finally, it is of interest to compare the values of Dy” with those obtained by _ 
direct linear Birge-Sponer extrapolations for the ground-state vibrational 
levels. It is found that the latter are only some 57% and 63% of the above ~ 
values for GaF and InF respectively. This supports Gaydon’s view (1947) 
that ionic forces contribute appreciably to the binding in this group of molecules, 
and it must be concluded that linear extrapolations of vibrational intervals are 
quite without value in these cases. 


§5. DISCUSSION: TIF AND AIF 

The values of the constants derived in the present work are collected together 
in table 11. A comparison of these figures with those for the other gallium 
and indium halides (see, for example, Herzberg 1950) shows that GaF and InF 
have very much the properties that would have been expected. ‘The binding 
in the states of the fluorides, as measured by dissociation energy or by the force 
constant, is appreciably stronger than for the corresponding states of the other 
halides. ‘This effect is particularly marked for the c1II states, where, for example, 
the force constants for GaF, GaCl, GaBr and Gal are 2:58 x 10°, 0-20 x 10°~0 
and 0 dyn/cm respectively. ‘The indium halides behave similarly. On the other 
hand, the multiplet splittings, $11, —°II)+, vary little throughout the series of 
gallium halides (the values, —F to —I, are 334, 332, 370, 330cm~) or indium 
haliles (816, 796, 786, 649 cm~). 

Finally, the present results throw some light on two anomalies existing in 
th : literature of the molecules of this group. 


(1) The Electronic States of TLF 

The *If—*X absorption systems of Ga’, InF and 'TIF appear to be closely 
analogous: for each molecule, the system at longest wavelength consists of 
single-headed bands, while the shorter-wavelength system shows two heads 
in each case. ‘hus, contrary to Howell (1937), the minimum of the potential — 
curve for the °IIy+ state of TIF lies, as for the other molecules, below that of 
the *II, state. ‘Thus Howell’s assignments of the components for TIF must be 
inverted (see also Herzberg 1950). 


(ii) The Dissociation Energy of ALF 
The present results suggest that the value Dy” for AIF should lie somewhere 


between the value for GaF, 6-15 ev, and the ee less certain value, 8-4 ev, 
which can be derived by extrapolation of the vibrational levels in the state a1II 
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of BF (Chrétien 1950). ‘Two values have been put forward for D,’(AIF). 
Gaydon (1947), from an extrapolation of the ground-state vibrational levels, 
following Rochester’s analysis (1939) of the absorption bands, gives 
D* = ?1-8+0-3ev, but notes that ionic forces may be important, i.e. that the 
true value may be higher. A second value, D’ =6:3;ev, may be derived from 
the measurements by Gross, Campbell, Kent and Levi (1948) of the standard 
heat of formation of AIF gas. A revised vibrational analysis of the 1II —1X system 
of AIF, on the basis of emission spectrograms (Rowlinson and Barrow, 
unpublished work) leads in particular to a smaller value of x,”w,” than was 
derived by Rochester, and a somewhat higher value of D’, 4:4 ev, is obtained by 
the ground-state extrapolation. We should expect, by analogy with GaF and InF, 
that this figure represents about two-thirds of the true value of D’, i.e. D’~7ev. 
A better spectroscopic estimate of D” may however be obtained by extrapolating 
the vibrational levels in the 'II state; this gives DY =7-2ev. It is thus established 
that the dissociation energy of AIF is somewhere around 6 to 7ev, but there 
is still a small discrepancy, 0-8; ev, between the thermochemical value and what 
is considered to be the most reliable spectroscopic value. Levels in the 1II state 
are known high enough to give a lower limit for D” at 6-65 ev, 0:3 ev above the 
thermochemical limit. This fact suggests that there is a small maximum in the 
potential curve for the !II state of AIF. The potential curves for the c1II states 
of GaF and InF may be similar, but work on the thermochemistry of these 
molecules is needed before any definite conclusions about the existence or 
otherwise of potential maxima in these states can be drawn. 
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ABSTRACT. 'The determination of the product pv for particles leaving tracks in photo- 
graphic emulsions is discussed on the basis of measurements of the scattering in the tracks 


ee ae 


~~ 


met 


made by the ‘coordinate’ method (Fowler 1950), allowing for errors of measurement due ~ 


to grain size. The use of Fowler’s method with the comparatively large interval lengths 
used in practice is shown to be justifiable in spite of its neglect of much information 
contained in the track. (This neglect increases the root mean square error in the estimation 
of the scattering over the minimum attainable error by a factor less than approximately 
2-6.) A possible method of estimating pv, which approximates to the ‘best possible’, 
based on the use of mean squares, is considered, and the root mean square relative error in 
the estimation found, the necessity for cutting off the scattering distribution at large angles 
being taken into account. A possible method of estimation based on the use of mean moduli 


is also discussed. 


§1. INTRODUCTION AND SUMMARY 

HEN a charged particle passes through matter it undergoes multiple 
coulomb scattering and, over a region in which the loss of energy by 
the particle is negligible, the probability distribution of any quantity 
associated with this multiple scattering depends essentially on the value of the 
product pv, where p is the (relativistic) momentum of the particle and v is its 
velocity. For an ionizing particle leaving a track in a photographic emulsion 
a measurement of the multiple scattering in the track therefore affords a means 
of estimating pu for the particle: in conjunction with a measurement of its 
velocity (specific ionization) or energy (residual range), this estimate can be 
used to give complete information about the mass and energy of the particle. 
Scattering measurements are particularly useful for fast particles which do 
not come to rest in the emulsion or which produce a nuclear event before coming 
to rest, since for these the range-energy relation cannot be used. It is thus of 
importance to discuss ways of measuring the scattering and of estimating pv from 

the measurements, and to find the probable error in these estimates. 
In discussing this problem, specific allowance must be made for the 
unavoidable errors caused by the finite size of the grains forming the track; 


a measurement of the position of the centre of a grain will give information | 


about the position of a point on the actual path of the particle only within a 
certain accuracy, since developable grains off its direct path may be formed by 
the ionizing particle. It will be assumed here that these errors have a gaussian 
probability distribution, of width approximately equal to that of a grain diameter, 
and that this distribution is independent of the multiple scattering. No other 
sources of error will be considered, and it will be further assumed that the loss 
of energy by the particle is negligible along the portion of its track that is measured. 


* Now at Department of Mathematics, University of Southampton. 
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The theory of multiple scattering has been discussed by several authors 
(Williams 1939, Moliére 1948, Snyder and Scott 1949, 1950). The results of 
Snyder and Scott have been found to be in satisfactory agreement with 
experiment (Corson 1950) and will be used below when numerical estimates of 
quantities associated with the multiple scattering distribution are required. 

The maximum amount of information about a track would be obtained if 
the position of every grain init were measured. Apart from the practical difficulty 
of making such a measurement, the analysis of the results would be complicated 
(Moyal 1950). ‘This paper will therefore be concerned with the treatment of the 
information obtained from the ‘coordinate’ method of measuring up tracks 
developed by Fowler (1950), which is the most accurate so far described and 
which can in principle be extended to yield almost all the information contained 
in the track. 

In Fowler’s method the projected angles between successive equal chords 
of the track are measured. ‘These measurements give the true multiple scattering 
angles combined with errors due to the finite grain size, and since these effects 
are statistically independent, the obvious method of obtaining information about 
the multiple scattering is to use mean square values of the scattering angles. 
This introduces a difficulty, however, since, because of the Rutherford scattering 
law, the mean square value of the chord angle diverges logarithmically in small 
angle approximation, and to give it a well-defined meaning some method of 
cutting off the observations at large angles must be used. ‘This difficulty would 
not arise directly if mean moduli were used instead of mean squares, but one would 
then be faced with great analytical difficulties in removing the effects of the 
errors of measurement, and some sort of cut-off would probably still be necessary 
in order to minimize the statistical errors. It has therefore seemed most practical 
to discuss the problem with the use of mean squares and to introduce an explicit 
method of cutting off the observations. 

In §§2 and 3 of this paper, the problem is discussed on the (incorrect) 
assumption that the mean square quantities introduced are well defined. 
§2 is concerned with the ‘minimum possible’ error in the estimation of the 
scattering, while in $3 the corresponding ‘best possible’ equations are discussed. 
Their form leads to the consideration of a particular method of estimation which 
is shown to approximate to the ‘best possible’ when the interval length used in 
the measurements is not too small. ‘The discussion also leads to a justification 
of the use of Fowler’s method with comparatively large interval lengths, in 
spite of the consequent neglect of possible information about the scattering 
which might be obtained from the track. 

In the remainder of the paper the effect of the introduction of a cut-off is 
considered. In §4 a cut-off criterion, depending on a cut-off parameter, is laid 
down and its effect on the results of §3 is discussed, resulting in the postulation 
of an equation for determining the cut-off mean square chord angle and an 
expression for the mean square error of this estimation. § 5 is concerned with the 
estimation of the product pu and the determination of the best conditions for 
measurement, including the best value to choose for the cut-off parameter. 
Finally, in § 6, a systematic error in the estimation of the scattering, due to the 
method of cut-off used, is discussed and shown to be negligible, and in §7 a 
possible method of using mean moduli is considered. 
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Notation. Apart from symbols defined in the text, the following general 
notation has been used. If £,, is an angle derived from the mth of the measurements 
of angle made on the track, its mean square expectation value is denoted by do, 


while the expectation value of the product £,8,,, is denoted by 6, The — | 


corresponding mean values taken over the actual measurements are indicated — 
by enclosing the symbol in triangular brackets, thus ¢b,). Finally, the relative — 
root mean square error in the estimation of some quantity x from the measure- 
ments is denoted by F(x), so that the actual root mean square error is xF(x). 


§2. THE ‘MINIMUM POSSIBLE’ ERROR* 

(i) Correlations between the Chord Angles and Statistically Independent Variables _ 

In Fowler’s method the track, which is nearly straight for a fast particle, 

is aligned along the x axis of a microscope stage and its y coordinates, y,, are 

read off at equal intervals 7 along the x axis. Since the scattering angles are 
small, the quantities 


1 
Ba= = (nt ~ 2Yn + Int) Senate (1) 


are the projected angles between successive equal chords of the track. In 
accordance with the programme outlined above, one writes 


Bu=%a+nr..>0) bon acute (2) 


where «, 1s the true chord angle due to multiple scattering, and «, is the error 
in its measurement. ‘Then, since the «,, and «, are statistically independent 
(using the notation described in the introduction), 


by = dye, “OLS GETS 10) 26 eee (3) 
on the assumption that the mean squares are well defined. With this assumption 
it can be shown (Appendix A) that, independently of the exact form of the « and e 
distributions, 


by =a) /4—2eq/3 3 ba eg/6> De 2) ee (3 a) 


and that a) is proportional to r and eé) to 7 *, so that the ‘signal/noise’ ratio 
v = d)/éq is related to the interval length by 


vs (¢/7)2o' 1h teem aa) Ge eee (4) 
where 7, is the value of 7 for which v= 1, and depends on the characteristics of the 
track under examination. 

The calculation of probable errors is greatly facilitated by the introduction 
of certain linear combinations of the 8,, which are almost statistically independent. 
If one adds to the set of measured angles £,, »=1, 2,...., N the fictitious 


angle fy =0, one may regard this whole set of angles as coefficients in a Fourier 
expansion of a set of variables 


2 \a . 2anp 
= == Sy —— ,p= SIRE So TOLER 
by (ws i) 2 Pa eageee Ouse Se (5 a) 
s Pa Se . 2anp 
then Bo=(xa3) Rett sin 7? 7 Oy le dacs ay Nee cen (5 d) 


and these equations allow one to express results either in terms of the ¢, or the 
measured angles f.,. 


* Some of the results for the gaussian case given in this section were derived in unpublished 
calculations by Dr. M. Krook. 


On the Utilization of Multiple Scattering Measurements 6.43 


From eqns. (5) and (3) it is readily shown that the correlations between 
the ¢, are given by 


4 . Qap oe 27g 
$rbg= 0,” Sngt N+ Wai 2 Sorta Ey ce Noe pO We aa? kav enerens (6) 
where 
2np 4np 
——4 # 2 
®)=0; 0, by + 2b, cos TF +2b, WS lesa) p> Pal BN ciate tte oes (7) 


If the second term on the right-hand side of (6), which is of order 1/N, is 
neglected, the ¢, may be regarded as statistically independent, with a mean 
square value given by (7) which becomes, on substitution from (3 a), 


2rp 2rp \2 
Bes, 1s : 2 = 
® 7 =a, (2+ cos 5) + 2e ( cos Pi] 


= : “| (1-¢ os er) + Ai 7 (2+ cos wa) |: DEG fe (8) 


This approximation materially simplifies the calculation of probable errors, 
and is quite justifiable when N is large, since one does not require great accuracy 
in the estimation of probable errors. 


(ii) The ‘Minimum Possible’ Error 
In order to judge any particular method of estimating a) one must know 
what the inherent limitations on the accuracy of estimation are, as determined 
by the probability distribution of the measurements. Using the ¢, variables, let 


P(d1s «+ + +» Py [40s 0) 41... - dy 
=constant . exp{—42(¢,, ....,¢y |, &9) }ddy....-. Gig stirs (9) 


be the probability of ‘observing’ ¢,,...., dy for given a and é. ‘Then, from 
the principle of inverse probability, the se pes of a) and ey having specified 
_values for a given set of observations ¢,,...., dy 1S 


P(ao, £9 lpr, +++) hy) dagdey=A exp, — 38(b1,--- shy aq; €) }day deg, 


where A contains the prior probabilities of aj and é). : 

The most probable values of a) and é) for given ¢,,...., dy will be those 
which make (10) a maximum. ‘That is, assuming that the prior probabilities 
of a) and éy are constant or only slowly varying functions, the most probable 
values of a, and @) are the solutions of 


Opi0d ee Oce0n de, =). Bae) ot 0 tenn (11) 
These are the ‘best possible’ equations for estimating a) and é), making the 
maximum use of the observations. 
Suppose the solutions of (11) are aj™, e9™, and put x=a)—a ™, y =ey— e™. 
Then, in virtue of (11), 


(a0, &) =g(ag™, e6™) + 5| (x2 +y~) | 4+ 
E40) &0) =8\ 40> &o 3 Oa, de, aes 


=f, + (ax®+2bxy+cy")+.... 
where A=40"2/0a_?; b=40"g/Oaydey; C=40%g/Oey?. =. ss ae (12) 
Thus P(dp, y|¢1, .. - -> Py) day dey = Py exp{ — 4(ax* + 2bxy + cy")} dx dy. 


42-2 
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Using this expression for P, it is clear that the mean square error associated 
with the ‘best possible’ determination of ap is given by 


[apFy(a))P=x2=[a—-Bfef* ae (13) 


BP |! 
The definition of a, b and c by (12) depends explicitly on the ‘observed’ values ~ 
of the ¢,. It will, however, give rise to only a small error if these ‘observed’ ; 
values are replaced by their expectation values after the differentiations have — 
been performed, and one may assume that this has been done in (13). 

To proceed further one needs to know the function g(¢,,...., dv], o) _ 
determined by the combined probability distribution of ¢,,...., dy. The © 
exact expression for the distribution cannot be put into manageable form, but © 
one may obtain an approximate evaluation by assuming that the «, have a 
gaussian distribution. On this assumption, the ¢, will also have a gaussian 
distribution since each ¢, is a linear combination of thea, ande,. Then, because — 
the ¢, are taken as statistically independent, 


N 
Plrs + + +» P| Gos 0) = (27) ™ oo exp {— 365 /®,°} 


N 
which*gives.  g(¢i,.. << «, dy |p, 9) = 2 (loge OF oO eee (14) g 
e p= c 


where ®,*, given by (8), is a function of a and e. The quantities a, b and c 
can be found from (14) and, on replacing the observations ¢,? by their expectation 
values ®,” after differentiation and transforming sums over p to integrals, one 
obtains from (13) 


ues 1 
Fi"(%0) = 9 (N+ Typ 20 — Fir"/Loe)™ bon ok (15) 
where 
— 1p = (2 cos kl cos ky 
Le) a oe [( — cos RP + 34(2 +cos AYP Ai sae (16) 


To discuss the variation of the ‘minimum possible’ error as given by (15) 
with the interval length one must express the number of observations in terms 
of 7 or vy. Suppose that the total length of track available for measurement is T. 


Lace Nels Boa v8Crim)) 2 eee (17) 
and (15) can be written as 


(T/1)Fy2(4y) =32v-*[Iog—Jy32/Ipa 2/9. wn en ws (18) 


For a given track, i.e. given T and 7, (18) exhibits the variation of F.y(ap) 
with v=(z/79)?; the quantity (7/r,)!?Fy(a)) is plotted as a function of r/z) in 
fig. 1 (the broken curve I). 

It may happen that for a number of similar tracks in an emulsion an accurate 
estimation of the ‘noise’ ey is possible. In this case the ‘ best possible’ equation 
for dy) is simply dg/da)=0 and the corresponding ‘minimum possible’ error 
is given by (18) with the term in J,,”/Z), (which arises from uncertainties in the 
estimation of é) omitted. ‘This ‘minimum possible’ error is plotted in fig. 1 
as the broken curve II. The two errors coincide for small values of T/T) as one 
would expect, and both show a very slow decrease in the ‘minimum possible’ 
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error with interval length. Evaluation of the integrals in (18) for 7/7y)<1 shows 
in fact that for small interval lengths 

(Pag) eR (Gg) 2s ir) ee Rn eens (19) 
This slow decrease of the minimum possible error with small interval lengths 
will be shown, at the end of the next section, to justify the use of Fowler’s method 
with the comparatively large interval lengths that are convenient in practice. 


§3. THE ESTIMATION OF a, 
(i) The ‘Best Possible’ Equations 


The explicit expression (14) for the function g(¢,,...., dy|a@, &) enables 
one to write down the ‘best possible’ equations (11) for a gaussian multiple 
scattering distribution. These equations are 
Oe ee odes SAB! pe SN OR fa as 
04 Ga, [¢,2-9,7]=0; -—=z-5 Ge, \Po =D 2A Oe oe a (20) 


? = a 
dey D ®, 


and, using the expression (8) for ®,?, may be written in the form 


(aoa rey SHE nee (21) 


with a similar expression for (@),.,. ‘The form of (21) and the explicit expression 
for the k,(v) obtained from (20) both result from the assumption of a gaussian 
scattering distribution. The ‘best possible’ equation for a general non-gaussian 
distribution will differ from (21), but one may assume that, provided the 
distribution has a well-defined mean square, an equation of the form (21) with 
suitable choice of the coefficients k,(v) will give results lying close to the ‘best 
possible’ with a statistical error of the order of magnitude of the ‘minimum 
possible error’. The fact that (21) must become an identity when (a)),,, and the 
‘measurements’ ¢,” are replaced by their expectation values a) and ®,” imposes 
general conditions on the k,, namely that when in the right-hand side of (21) 
¢," is replaced by ®,” the coefficient of aj must be unity and that of é) must 
vanish. These conditions are, from (8), 


1 2rp Eo np a5 
N+ 1)» Pol) | 2+ Soe | =. ae) [1- cos | =={), 


and imply that (expanding k,(v) in a Fourier series) the ‘best possible’ 
equation (21) is necessarily of the form 


1 ‘4 2rp 4np 
(4o)ess= TT a3 tte cos Wz + (44-6) cos 5 
4 2rpn : 
+ ae Cos wo |e A soda (23) 


where the actual values of the coefficients p(v) and A,(v) will depend on the 
scattering distribution. Application of the transformation (5) to (23) gives the 
general form of the “best possible’ equation for the estimation of a) in terms 
of the actual measurements f,,: 


(age Ge (bo) + wkbr) +44 6){Ba)+ BD ABe dbs vse (24) 
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(Since the derivation of (23) neglected the small correlations that do exist 
between the ¢,, the application of (5) to (23) will produce some spurious terms ; 
these are easily recognizable and have been dropped from (24).) 


(ii) A Possible Method of Estimation 


One would expect the terms in (24) involving (b,) with r>2 to have a small — 
effect on the estimation, since they have zero expectation value, and thus that a 
good approximation to the ‘best possible’ equation (24) would be* 


4 
(4o)est = rer {{bp) + Xb) + (44 —6){bo)}, we eee (25) 
with the parameter p suitably chosen. A particular form of (25) corresponding 
to w=8/5 is 
(gett i= Col G8 aimee (25 a) 


where b,” refers to the angle between adjacent equal chords of the track of 
length 27. This form of the equation can be directly derived by considering 
the dependence of a, and eé) on the interval length, and is particularly simple to 
use in practice. 

The root mean square error associated with the estimation of a) by (25) can ~ 
be evaluated by writing the equation in terms of the ¢,,; neglecting small terms 
of order 1/N due to end effects, the equation is of the form (21) with 


ee aap 4ap 
ky= ze, | 1+ 1 008 yy + du —6) cos 5 |, ip a. (26) 


and consequently the mean square error in the estimation of a, by (25) is 


(Ne D 
since the ¢, are taken to be statistically independent. 


Now ¢, is the sum of two statistically independent parts depending on the 
a, and e, respectively; 6,=¢,(«)+¢,(<). Hence 


by'—O,4=$,'(a) — Da) + $Me) —OpNe) +40, 7(a)P 26), -.--. (28) 


and, since the «,, are assumed to have a gaussian distribution, ¢,4(<) =30,(e): 
Thus the part of (27) depending on the last three terms in (28) can be evaluated 
by using (26) and the appropriate expressions for ®,>(x) and ®,(«) obtained 
from (8). Since the gaussian assumption is mot made for the «,, however, no 
such simple calculation is possible for the remaining terms of (27). The best 
way to find these is to rewrite them in terms of the «,; the various expectation 
values involved (of the type «,4, «,3«,41, etc.) can be found by the methods of 
Appendix A, and expressed in terms of a) and the quantity 


[4,F(a,)*2[uoak= == Dhige tom eae (27) 
p 


3 =at/c%9%, 
and the summation can then be carried out. 
The final expression obtained for the relative mean square error in the 
estimation of a) by (25) is 
1 8 
F'(a)) = —— ——5 (4? - 
(a) Wei G@tp) (4 — Bat Ch> 5 Seen (30) 


* This equation was suggested by Dr. F. J. Dyson. 
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Bose lun 361 
where A= SE tea fail ap okeo pt 
354 1208 
Bh LE rele mg caine gene ho Cal, 5 (31) 
one ial 7 2304 L012 


Pe a foes 
The actual shape of the scattering distribution thus affects the error only through 
the occurrence of the quantity « in A, B and C. 


(iti) Comparison with the ‘ Best Possible’ Method 
Expressing the number of observations N in terms of T, 7) and », as in (17), 
one may write (30) in a form analogous to (18): 


(T}z,) F(a) = a PRA Gonr Ce ety eee (32) 


Comparison of (32) with (18) will then show how near, in fact, the estimation 
of a by (25) approaches the ‘best possible’ estimation. This comparison is 
made in fig. 1, taking « =3 in (31), the value appropriate to a gaussian scattering 
distribution. In the figure the solid curve A shows the envelope of (T/79)"?F(ap) 
calculated from (32) when the best value of 1 is chosen for each value of v, while 
the curve B is obtained from (32) with » =8/5, giving the error in the estimation 
of ad by (25a). 


(T/1,)" Flay) 


T/Ty= VP 


Fig. 1. I: The ‘minimum possible’ relative r.m.s. error in the estimation of a) when ép is 
unknown ; II: the same quantity when é is known from other measurements. 

A: The relative r.m.s. error in the estimation of a) from (25) when the parameter p is 
chosen appropriate to each value of v ; B: the same quantity when in (25) w=8/5. (All 
these curves assume a gaussian scattering distribution.) 

Figure 1 shows that (25) gives results close to the ‘best possible’ provided 
7/T)=0°7, while for smaller values of 7/7, it becomes very bad. ‘The general 
form (24) of the ‘best possible’ equation indicates that the estimation can be 
improved by adding terms in ¢b,) with r>2 to the right-hand side of (25), and 
an analysis similar to that given above shows that if a term Ab; ) is added to (25) 
the resulting estimation, with proper choice of ~ and A, remains good down to 
7/7) 0-45, after which the error again increases rapidly. Because of the slow 
decrease of the ‘minimum possible’ error with the interval length, this improved 
method of estimation can only increase the accuracy of an estimation based on 
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the best use of (25) by 8%, and it is evident that to obtain any marked 


improvement on (25) very small interval lengths would have to be used, with 


a correspondingly great increase in the labour of estimation. Now the finite 
spacing of the grains forming the track imposes a lower limit t,;, on possible — 
values of 7 and consequently on the minimum possible attainable error. If one — 
uses (19) for this minimum error, taking tj, to be the average spacing of the 
grains in the track and 7, the interval length corresponding to unit ‘signal/noise’ 
ratio, to lie between 100-500 , it appears that the greatest possible improvement — 
that can be made on estimates based on measurements with 7/T)=1 is a reduction 
of the root mean square error by a factor lying roughly between 1-8 and 2:6. — 
This comparatively small decrease in the error does not justify the elaborate 
technique that it would require for the estimation of a), and for practical purposes 
one may therefore conclude that Fowler’s method with the comparatively large 
interval lengths implied by 77, is quite justifiable, and that for these intervals 
an estimation based on (25) will give good results. 

This result has been obtained on the assumption of a gaussian scattering 
distribution, which is incorrect. Its qualitative aspect, however, should be 
largely independent of this assumption and may reasonably be applied to the ~ 
cut-off scattering distribution. 

In order to find the value of v corresponding to a given interval length, it is 
necessary to estimate é. This can be done in a way analogous to (25); the 
equation for é) corresponding to (25 a) is 


(Co)est = 7{2€B0)— <8" J. nan ee (33) 


§4. THE CUT-OFF CRITERION 
(i) The Effect of the Cut-off on Equation (25) 

The results obtained above were based on the assumption that the quantities 
@ and é are well defined. As pointed out in the introduction, however, this 
assumption is quite unjustified, and in order to give a well-defined meaning to 
the mean square values of « and 8 some sort of cut-off must be applied to the 
observations. When a cut-off is imposed, the equations (3) on which the results 
of the previous paragraphs were based are no longer exact and it is necessary 
to investigate the effect of the cut-off on them. 

If W(a«)d« is the distribution of the multiple scattering chord angle, 
Snyder and Scott* have shown that for « measured in units of the quantity 


p= Ka? pee La Es Pca (34) 
(where K is a constant depending only on the properties of the emulsion and 
the charge of the ionizing particle) so that «=px’, one may write 
Wa) da =w(a") da’, 
where the distribution w(«’) is very nearly independent of the interval length. 


Suppose that the f distribution is cut off at B= +6. 'Then denoting angles 
measured in units of p by a prime, the cut-off f’ distribution is 


P(g’) dea 1} UCB’ —e') dB; IB<O, 3 
ae | 0 ; IB I>e, x: 


* See Appendix B for a discussion of the Snyder—Scott results. 
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and the mean square value of f’ over this cut-off distribution is 
(oe) 0’ co 0’ i 
bo =| de’ aprsrra(e'yw(p’—e)/[” ae’ a'gleeu6'—c'), ... (86) 


where q(e’) is the probability distribution for the error e’.. Now q(e’) has a sharp 
maximum for e’ =0 and falls off very rapidly for large values of <’. Thus if the 
cut-off angle 6 is chosen so that q(0’) is small and if the ‘signal/noise’ ratio is 
not too small, one may assume that any value of the measured angle that is 
comparable with @ is due mainly to an exceptionally large multiple scattering, 
the error of measurement contributing only slightly to its size. Then a reasonable 
approximation to (36) is obtained by expanding the integrands in a Taylor series. 
about «’ =0 and retaining only the first few terms. The result of this expansion 


ives 
E Oye fee, PG = vive ee (37) 


where a,’ is the mean square value of «’ when the «’ distribution is cut off at 
angle 6’, v,=d),/e) and y is a quantity depending on 6’ and »,, which is small for 
large &’ and v, not too small. (For the values of 6’ and », finally chosen as the 
best for measurements, y ~0-1.) 

This expression for by,’ is to be compared with the first of equations (3) 
which, on writing @)=a,)/v, becomes by)’ =a) [1+1/v]. The only difference in 
the form of the two equations is the presence of the small number y in (37), 
so that it appears that the results obtained on the basis of (3) will be approximately 
correct for the cut-off distribution, if a) and v are replaced by their cut-off values 
dp, and v,, the cut-offs for the « and f distributions being identical. 

In particular, since w(«’) depends only slightly on the interval length, it is. 
reasonable to neglect the change in it when the interval length is doubled. 
With this approximation, y is a constant, a),=p"d», is proportional to 7 and v, 
to 7®, so that the form (25 a) of the equation developed in the previous paragraph 
for estimating a) may be retained for the cut-off distribution with confidence: 


(4oc)est = *[4 Does > aa (boc pls A BIGOG (3 8) 


where b,,” refers to measurements made with interval length 27. ‘This equation 
will give a good approximation to a),, and on the basis of the arguments developed 
above, the value of a, estimated from it should have a mean square error close 
to the minimum possible, provided that the value of », for the measurements 
is correctly chosen by suitable choice of the interval length. (From the way 
in which (25) was derived, it is evident that the procedure for using (38) is to 
cut off the original sample of values of £,,, and to form the ,,” from those that 
are retained.) 

The expression (30) for the error associated with an estimation by means of 
(25) may also be retained if v is replaced by »,, N by N, the number of 
observations kept in the cut-off set of measurements, and « by «,, its value for 
the cut-off « distribution, while » is put equal to 8/5, the value appropriate 
to (25a). The approximations involved here are the postulation of equations (3) 
with a, replaced by ap,, and the neglect of the local destruction of correlations 
caused by the rejection of a particular f,,. ‘These approximations are reasonable 
since no very accurate estimate of the error is required. 
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(ii) Postulation of a Cut-off Criterion 


The method of cut-off used is rather arbitrary. The criterion discussed 
below is simple and its consequences easily evaluated. Define 


0 0 j | 
t= [_2°V(e) dx | We) de with @2=Aaj,  «.- se (39) 


where A is the cut-off parameter. (39) defines ay, uniquely for each value of A, 
and the estimation of a), from the measurements proceeds simply as follows, 
since one assumes that any large value of £,, arises from a large multiple scattering © 
deviation, with a comparatively small contribution from the error of measurement: — 
find (4p)e; from (38) for some arbitrary, small value of the cut-off and determine — 
67.4=A(Aoc)estr Then reject all the 8, for which 8,?>6,,, and recalculate 
(2 c)esp and 6?,., with this new cut-off. Carrying on this process will eventually 
result in a steady state, and the final values of (dp,)og, and 67 4 are taken as 
estimates of a), and 6. 
Suppose that da), is the error in this estimate of a,. ‘Then 


Sdao= (Saya), + (day, aoeyee#, = ane (40) 


where 86? is the error in the estimate of 6, and (8a),), is the error due to statistical ~ 
fluctuations in (38). Now 67,.=A(dpc)est) therefore 66?=Ada,, and substituting 
for 60? in (40) gives 

day, =[1— \(dag,/d6”) |" (Sage) s- 


‘Thus the relative root mean square error, F,(d),), in the estimation of a, defined 

by (39), from (38) is 
F.(a6,).= | De A dagidik FG aes eee (41) 

where F(a.) is given by (30) with the changes described in (i) above. This 


equation will be used in the following section as the basis for deciding on the 
best value to choose for the cut-off parameter A. 


§5. THE ESTIMATION OF pv AND THE BEST CONDITIONS FOR 
MEASUREMENT 


Expressing the equations (39) in terms of angles measured in units of p one has 


Ay, =p” G(8'), 
0’ 6’ 
where G6’) = | a en( cr’) de’ / | w(a") da’, O2=AG(O). —....e (43) 
—6 —i97 


Now the distribution w(«’) depends on the interval length by being a slowly 
varying function of m, the average number of scatterings suffered by the ionizing 
particle in the traversal of one interval length. This number can be determined 
if the velocity of the particle is known, and so can be estimated independently 
of the scattering measurements. Assuming it to be known, eqns. (43) can be 
solved to give the value of 0’ and G(6’) corresponding to the value of the cut-off 
parameter A chosen. Hence, by substituting the estimated value of ap, in (42), 
p and so pw can be calculated. 

The error dp in the calculation of p from (42) is related to the error Sap, in 
the estimation of ag, by dp/p = $5a9,/49,._ Hence F(pv) the relative root mean square 
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error in the determination of pv when ag, is estimated in the manner indicated 


in §4 is 
F(pv) 7 F(p) a BF (Apo), Tenner (44) 


where F,(a),) is the expression (41). Using (43), the term in (41) depending on 
day,/d6? can be expressed in terms of 6’, w(6’) and X and from (30) and (41) 
one finds 


F(pe) =41—0-1)8'«(@)/00 PF (sq) Fla Bu cys, (45) 


6’ 
where (6) = | ee eo ee (46) 
= 6" 


Here A, B and C are given by (31) with v replaced by », and «x by x,, and 
je =8/5 the value appropriate to an equation of the form (38). From the definition 
of p(@’), the number of £,, retained after the cut-off has been applied will be 
N,=p(@)N on the average. Then, writing N in terms of 7, and the total 
length T of measured track, and substituting the numerical values of A, B and C, 
(45) can be written explicitly as 


| 54/2 {p(0’)}12 1-28  0-27\ pe 
1/2 at eee SS aes ee Ty ee teak 
(The Oe) = ea AL (het S744 SE + SE) | 


This expression for the relative error in the estimation of pv depends on two 
parameters, viz. v, the ‘signal/noise’ ratio for the track and interval length 
considered and A the cut-off parameter. For any particular track, the 
‘signal/noise’ ratio can be made to take any value (within limits) by suitable 
choice of the interval length. The parameter 4, on the other hand, should be 
fixed once and for all for every track, to enable comparisons of different tracks 
to be made without the necessity of calculating pu from ap, each time. 

To find A,,, the best value for the cut-off parameter A, one first, for given A, 
finds the minimum of the right-hand side of (47) with respect to v,; Ay 1s then 
the value of A for which this minimum is least. Carrying out this procedure, 
one finds (taking as a reasonable value of the average number of scatterings per 
interval length n=1-5 x 10°) A, =18, and that A can vary somewhat about this 
value without affecting the error appreciably. 

Using A=18, the quantity (T/7))"?F(pv) is plotted as a function of (7/T») 
in fig. 2. From this plot one can find the probable error in the estimate of pu 
when the value of v, corresponding to the interval length used for the 
measurements is known. It is evident that one should always try to work with 
interval lengths such that 7/7)20-75. If this is done, the error, in terms of the 
total number of measurements JN is 


F(pv) = 1-73-12, a8) 


For purposes of comparison, the ‘minimum possible’ relative error in the 
estimation of pv, calculated from (18), is plotted in fig. 2 as well (curve B). Since 
the gaussian assumption made in deriving (18) probably does not affect the 
numerical values appreciably, the difference between the two curves (as compared 
with the similar curves B and I of fig. 1) shows the effect of the extra uncertainties 
resulting from the necessary introduction of a cut-off. 
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§6. SYSTEMATIC ERROR IN THE ESTIMATION OF 4), 


The method developed above for the estimation of a, is valid in the sense 
that the estimate of ag, will tend to its correct value as the length of the track, 
and hence the number of observations, increases. One would expect, however, 
that if the estimates of a), obtained from a large number of similar finite tracks 
were averaged, the result would differ systematically from the true value of a, 
for the tracks, since the cut-off used for any one track depends on the scattering 
observed in that track. 


(7/1)? F(pv) 


T/T, =p 
Fig. 2. A: The relative r.m.s. error in the Fig. 3. Illustrating the notation of 
estimation of pu by (38) and (42) when Appendix A. 


the large angle scattering is cut off. (Based 
on the Snyder—Scott distribution.) 

B: The ‘minimum possible’ relative 
r.m.s. error in the estimation of pv for 
gaussian scattering with no cut-off. (For 
the use of these curves note that 
T|7)=N(1/70):) 


Suppose that one has a very large number of similar tracks of the same 
energy and that a,, 0, are the estimates of a), and the cut-off angle @ for the 7th track. 
Then, since for a large number of tracks the statistical fluctuations in the 
estimation of a, from (38) after the cut-off has been applied will — on the 
average, one may write 


a= (6) Boe a (49) 
me f(0)= fea) da [ 20); 7(6) = [We de. wae (50) 
Then aaa + FZ o2-0 P44 aap O + Petit (51) 


and averaging over all the tracks, as the number of tracks tends to infinity, the 
systematic error in the estimation of 400 is 


Ado. = 2 ay —~, (8622, pt er, ee ee ee Oa ta (52) 
since 66," = 6,2 — 0” will vanish on the average. Now from (49) 
[56,7] am [Sa,]2 i r?[ docF (0c) |? 


ee qd? 
Hence AG oo/ 400 = $r2 Ta oF 02) ae ter (53) 
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and since for large 0W(0)~1/68 
— m\ r / 12 dw 

AG w(0’) +6 WW 


writing @ Me W(@) in units of p, and putting p(?)=1, which is a reasonable 
approximation for large 6. Under the best conditions for measurement, the 
results previously obtained give F,(ao,)=3-46N-1?, A0’w(0’) 0-23, so that 
Daye) ao — 1-4/N. 

Taking N=50 as a reasonable value for the average number of observations 
made on a track, Adp,/a),=—3%. It therefore appears that while there is a 
systematic error in the estimation of a), using the cut-off procedure developed 
above, this error is small and is negligible in comparison with the inherent 
statistical error of the estimation. 


Fo'(Goc) = — BAG’) Fo*(Aoc)» ve vee (54) 


$7. A POSSIBLE METHOD OF ESTIMATION FROM MEAN MODULI 


The method of estimating a), developed above demands the calculation of 
mean squares from the observations. In practice, however, it is much simpler 
to evaluate mean moduli and it would therefore be advantageous to have a 
method for estimating the mean modulus of the scattering chord angle from a 
knowledge of that of the measured angles. The general treatment of this problem 
is difficult but a possible method can be postulated on the basis of the results 
obtained above. 

Denoting the square of the mean modulus of the true chord angle taken 
over a cut-off distribution by |a,|, with a similar notation for the measured 
angle f, eqn. (38) suggests that one might assume an equation of the form 


1 uM 
Tl Se a 0G =. Oe 9 FE aie eter eis 
[a lo= + [41200” I~ LBoe (55) 
for the estimation of [ape|- 
1 F 
Now | 2c] — al4 Poe V— [Poel lap. PS See. (56) 
Boel 4 | Boe’ | 
where r= 547+ | aed kde ele \, Pde 57 
FW, Naval Goad sie 


and the postulation of (55) will therefore involve a systematic error in the 
estimation of | a ,| of amount given by 


Bice iialets te (et: 9 GSR. (58) 


The calculation of [ requires the evaluation of |b9,| as a function of the 
cut-off angle and the ‘signal/noise’ ratio v,. If the contribution of the error of 
measurement to large values of 6 near the cut-off angle is neglected, one may 
write approximately 


(Soe)? = | a | ie de] a +€|W(a)q(c) / iy a | = aan 


ss w | te W(x) {2 | “a de+| ‘ <q(<) ae| Age Hake (59) 


where p(6) is defined by (50). (This approximation is equivalent to the neglect 
of the small number y in (37) in the analogous calculation of bo,._ It is reasonable 
for the purpose of estimating the order of magnitude of I.) The integrations 
can be carried out if a definite form for q(c) is assumed, the parameters occurring 
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therein being expressible in terms of ao,, and hence 6 and »,. The integration — 


in (59) was carried out for v,=1 using the Snyder—Scott results for w(x) (with 
n=1:5>10°) and assuming (i) a gaussian distribution and (ii) a rectangular 


distribution for the error «; in both cases the resulting value of I’ was found — 


to be less than 1%. 


>| 
f 


Thus it appears that (55) gives a reasonable basis for the estimation of 


| @.|=(|«|.)? provided that the ‘signal/noise’ ratio is of the order of unity and 
that the cut-off angle is large; the error of the estimation will be of the order of 
that found above for the estimation by mean squares. (If the cut-off 
criterion (39) is written in terms of mean moduli, the value A=18 for the cut-off 


parameter corresponds to a cut-off of approximately 5-5 times the mean modulus.) 
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ALP PEN DD 1 xe 
CALCULATION OF CORRELATIONS 
(i) Evaluation of a, 


Consider a typical portion of the track as shown in fig. 3. Using the angles 
as marked there, one has 
on =Pna—PnitYn NGOS (Al) 


and it is obvious that the angles ¥,, ¢; are all statistically independent except for 
those with the same sufix. Thus 


a= $2 +2226; a,=¢b—-W; a,=0, r>1. ...... (A2) 


On the assumption that the individual_coulomb scattering angle has a 
well-defined mean square, one can evaluate ¢?, ¥?, dx as follows. Suppose that 
during the traversal of one segment of the track the ionizing particle is scattered 
nm times, and that at the 7th scattering it is deviated through an angle 6, at a 
distance rx, from the end of the segment. Then, using the small angle 
approximation a A 
x=y[T= et Da Re ee Ne eee (A3) 

= 


i=1 
In (A3) the 6, are statistically independent, and since x, varies from zero to 
unity, the mean value of x?*1 is 1/(p6+1). Thus 


40 GXiXG = gn, 

ij=1 

= n ona 

gp =. LOO x= ene, ae a eee eee (A4) 
ij=1 

—_ nN =: — 

ae oy a 
= 6,0, =n0?, | 
ij=1 
and from (A 2) Qg=Fnl®; .a,= 3nl*=4a,, 3.) eee (A5) 


showing that a, is proportional to m and hence to 7. 
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(ii) Evaluation of e, 
Let z, be the error made in measuring the actual coordinate y,, of the track. 


Then 


ey = acer rex) Bo tely (A6) 

and since the 2,, are statistically independent with the same mean square 
ey= 5 (6% wader (A7) 
and Se CRIN ee ON ro Aer soe ne eM op (A 8) 


while (A7) shows that é is inversely proportional to the square of the intervah 
length. 


ARP BINDS 
ON THE SNYDER-SCOTT DISTRIBUTION 
Let the mean free path for coulomb scattering of the particle in the emulsion 
be A, and let p(y)dy be the probability that in an individual scattering the 
particle is scattered through a projected angle 7. Then the probability that 
after travelling a distance z into the emulsion the particle is moving at a projected 


angle 7 to its original direction and has suffered a lateral displacement x is 
(Snyder and Scott 1949) 


Wa, x)2) dy dx = — | a. | ” “ieee Wes | - = (hs+ ts) = 19} |. 
ioe (B1) 
where h{s)= de [1—q(u)]du; q(u)= le Pipe da ee ae (B2) 


(If this is integrated over x, the resultant expression for W(y| 2) is identical with 
that obtained by Moliére (1948) written in terms of the projected angle.) From 
(B1) one obtains the probability distribution for the chord angle « (Scott and 
Snyder 1950) as 

W(« |r) da = =| ds exp (ixs) exp { ee “I eas (B 3) 


RY Ss 


(Since small angles are assumed throughout, the chord length is equated with 
the length of the segment of the track cut off by the chord.) In order to evaluate 
this distribution, Scott and Snyder assume 
bea, No" Pa gS eal Bay 
p 7) AZ (72 + n92)32 ) r a2 pee ee! noe’ oc ep ete ( ) 
where 7 is a ‘screening angle’, Z the atomic number and N the density of the 
scattering atoms and Z’e the charge of the scattered particle. (‘These expressions. 
arise in Born approximation from a shielded coulomb potential of the form 
V =Z(e?/r) exp(—r/a).) ‘Then expressing angles in units of p=(2mn,2)"2, where 
n=7/A is the average number of scatterings suffered by the particle in traversing 
one interval length, Scott and Snyder give numerical results from which the 
% distribution can be found for given 2. 
Moliére (1947) has investigated in how far the actual state of affairs can be 
represented by the approximation (B4) and, basing his results on an analytic 
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representation of the Thomas—Fermi potential, concludes that if y is chosen — 
so that 
No = Nol 133" 0a") ae ee eee (B5) 


where yy =A/27a, with A the de Broglie wavelength of the scattered particle, — 


a the Thomas-Fermi unit of length and «=ZZ’/137f, then the mean square 
value of 7 over the distribution (B4) will agree with its mean square value over 
the actual distribution when both are cut off in a get where the Rutherford 
scattering law applies. 

With Moliére’s choice of the screening angle, the Snyder—Scott results are 
probably a good approximation since the unit of angle employed does not 
depend on 7, and the distribution is very insensitive to the value of »~1/n)? 
if n is large. 
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ABSTRACT. Examination of the relatively long-lived isotopes produced by pile irradiation 
of samarium using the proportional tube spectrometer leads to the elucidation of the decay 
schemes of 4Sm and !°Eu. ‘The conclusions are supported by experiments involving 
B-e~ and f-y coincidences, utilizing Geiger—Miller and scintillation counters. 


$1. INTRODUCTION : 
HE object of the work described in this paper was to obtain information 
| upon the long-period activities resulting from the bombardment of 
samarium with neutrons. ‘The radiations from these isotopes appeared 
to be well suited to a detailed investigation by proportional counter and related 
techniques, and it was hoped that results might be obtained which would clarify 
the decay schemes for the two isotopes, ©'Sm and ! Eu. 


Neutron capture by pure samarium results directly in the formation of the 


radioactive isotopes Sm, 'Sm, 85m and Sm (Inghram, Hayden and 
Hess 1947). '°Sm is a pure K-capture source of half-life 410 days (Butement 
1951). ‘Sm is a f-emitter of many years half-life (Inghram, Hayden and 
Hess 1947, 1950, Marinsky 1949). 85m is short-lived, decaying into stable 
13Bu. 1°°Sm is short-lived, decaying into B-emitting }°Eu of half-life 2 years 
(Hayden, Reynolds and Inghram 1949, Winsberg 1950). Since the capture 
cross sections of *Sm and its ultimate product '°Eu are high, a relatively very 
small amount of B-emitting 4Eu can also be expected. 
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Samarium in the form of oxide, spectroscopically pure, was irradiated at 
Harwell for eight weeks at a thermal neutron flux of 3-5 x 101°/cm?/sec. ‘The 
source was allowed to decay for several months. Preliminary absorption 
studies on the f-radiation showed that most of the B-rays were of energy less 
than 250kev, approximately the limiting energy of the f-spectra of 1°Eu. 
_ A further well-defined end point showed at about 75 kev, agreeing well with the 
end point of 44Sm as given by Agnew (1950) and Marinsky (1949), viz. 75:5 
and 74kev respectively. Evidence was obtained too of the presence of 4°Sm, 
which, being a pure K-capture source, is not an appreciable obstacle in the 
investigation of the radiations of !4Sm and Eu. The sources of radiation 
151Sm and 15Eu were present with relative intensities of roughly 1 to 4. It had 
been hoped that one or other of these radioactive isotopes would be predominant, 
but in view of the involved nature of the chemistry, the two activities have been 
studied without chemical separation. 


§2. APPLICATION OF THE PROPORTIONAL TUBE SPECTROMETER 
(1) Gamma Radiation 


A proportional tube, diameter 5-5in., active length 8in., filled with methane 
and argon (partial pressures 15 and 60 cm Hg respectively), was used to examine 
the energies of the soft y-rays from the source. The histogram of fig. 1 indicates 
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Fig. 1. y-ray peaks at 15-2 and 19-3 kev as measured in proportional tube spectrometer. 
Inset : spectrum below 10 kev using counter with aluminium cathode. 


the presence of radiations of energies 19-3 kev and 15-2kev. ‘The peak at about 
8 kev is due to copper K x-rays excited at the wall of the counter by the y-rays 
of the source, and was eliminated when a glass counter with an aluminium 
cathode was used to examine the region below 10 kev (inset to fig. 1). Concerning 
the two y-radiations of 19-3 and 15:2kev, Scharff-Goldhaber et al. (1950) have 
detected y-radiation of energy 21 kev emitted by Sm. ‘The radiation here 
estimated with energy 19-3 kev is presumably of this origin. Since they did not 
observe any radiation corresponding to the peak at 15-2kev in fig. 1, we are 
naturally led to the conclusion that this softer radiation is emitted in the decay 
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; 


of 155Eu (their source was a fission product, presumably frée of europium isotopes). _ 


Marinsky and others have found two y-radiations of °Eu of energies 85 and 


99kev. These radiations were measured by means of internal conversion — 


electrons in a magnetic spectrometer, and the very low energy radiation 


(15-2 kev) would presumably escape observation. It is natural to relate the | 


radiation of energy 15:2kev to the difference between 85 and 99 key, in which 
case the 85 and 15-2 kev quanta would be emitted in cascade. 

The inset to fig. 1, obtained with the second counter, shows a broad peak 
around 5kev which can be tentatively ascribed to the L x-rays emitted by the 
source (in K-capture and internal conversion). "Taking account of the absorption 
in the aluminium window, it must be concluded that the absolute intensity of 
this radiation is high. 

An investigation of the region 20-70 kev showed the K x-radiations centred 
about 40 kev. The source peak was compared with a calibration peak obtained 
by exciting the K« fluorescence x-radiations from a thin layer of samarium oxide 
with an X-ray equipment. The position and width of the source peak suggested 
that x-rays of higher and lower energy than Ka of samarium were emitted by 
the radioactive material. It is thought that the peak is mainly due to K-capture 
in 1°Sm and internal conversion of !°Eu y-rays (x-rays of Pm and Gd 


respectively). - 


(ii) Beta Spectra 


A thin source was mounted on a thin nylon film, and was inserted into a 
proportional counter similar to the first one described above. The gas mixture 


was methane to a pressure of 10 cm and argon to a total pressure of 34 atmospheres. — 


A magnetic field of 1400 gauss was applied (Curran 1950) and the counter was 
end corrected (Cockroft and Curran 1951). 

The spectrum analysis achieved with the help of a single-channel kicksorter 
showed discontinuities in the spectrum at 70-80 kev, ~150kev and ~250 kev. 
These various end points are associated with ©4Sm (70-80kev) and Eu 
(150 and 250kev). Marinsky’s work indicated end points for ®Eu at 154 and 
243 key. A very low intensity spectrum stretched beyond 250 kev, showing end 
points at 400 and 700kev. These were associated with ™4Eu. Previous work 
by absorption technique on 4Eu indicated that it emitted f-spectra of energies 
0:3 and 0-7 Mev (Hayden et al. 1949), 


§3. BETA-GAMMA COINCIDENCES 


The arrangement for measuring B-y coincidences consisted of a thin window 
B-counter (1-7 mg/cm?) and a lead cathode G-M tube. The latter was later 
replaced by a scintillation detector to improve the efficiency of y-detection. 
This detector was shielded by 5/64 in. of copper to stop electrons and soft x-rays. 
In each case the resolving time was 0:7 usec. The results obtained were similar, 
and a set using the scintillation detector is shown in fig. 2.- Here By coincidences 
per f-particle are plotted in terms of the thickness of aluminium between the 
source and the f-counter. 

This coincidence absorption curve receives a ready interpretation in terms 
of the decay schemes proposed in fig. 3. he first steep rise to approximately 
80 kev on the curve is due to the elimination of the f-rays of 41Sm, for the 
19:3kev y-radiation in cascade with the soft f-rays is not detected by the 
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y-counter. The rather less steep rise to approximately 150kev is probably 
associated with the absorption of the f-rays of higher energy (~250 kev) of Eu 
since they do not contribute to the B-y coincidences. The elimination of 
coincidences arising from the 150 kev spectrum of Eu results in the break 
in the curve at approximately 150 kev, assuming the proposed decay scheme is 
correct. The gradual falling off to lower values is mainly due to the 
contamination of the source with Eu, which gives rise to f-rays of energy as 


high as 0-7 Mev. 
The results of work on B-e~ coincidences were consistent with the above 


interpretation. 
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Fig. 2. Ratio of number of B-y coineidences to number of B-counts as a. function of ° 
aluminium absorber thickness. 
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Fig. 3. Proposed decay schemes of Eu and 1°tSm. 


Marinsky found that the lower energy f-rays of !°°Eu were coincident with 
y-radiation, but he concluded that there were no y—y coincidences. We believe, 
as the result of the above work, that the 99 kev y-ray corresponds to the difference 
in energy between the two f-spectrum end points and that the soft y-ray detected 
here (15 kev) corresponds to the difference between the energies of the harder 
quanta. one 
It is easy to explain the difficulty of finding y—-y coincidences with 5Eu. 
The 15kev quanta are not detected efficiently by most normal G—M counters, 
and, moreover, the presence of intense L-radiation and theoretical considerations 


suggest strongly that the quanta are very highly converted. However, a direct 
43-2 
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test of the scheme of fig. 3(a) was arranged. This test sought to establish that 
the electrons of the softer B-spectrum of Eu were followed by 15 kev 
y-radiation. A B-counter was covered with sufficient material to stop the f-rays _ 
of %1$m (75kev). Opposite this was placed a special G-M counter designed 
to detect soft quanta efficiently, containing argon to a pressure of one atmosphere 
together with ethyl alcohol (1-5cmHg). This second counter was covered with 
polythene of thickness 270 mg/cm2, adequate to stop all B-particles from the source 
while transmitting quanta of 15kev fairly efficiently. B~y coincidences were 
recorded as a function of the thickness of aluminium placed between the source 
and the y-ray counter. The results shown in fig. 4 indicate the predominance 
of two radiations. The lesser slope corresponds to the absorption coefficient 
in aluminium of 40 kev quanta, which is of the order of the K x-ray energy of 
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Fig. 4. Logarithm of the number of (8—soft y) coincidences as a function of y-absorber thickness. 


the source. The greater slope, after corrections for the harder radiation, 
corresponds to quanta of about 6kev, which is about the L x-ray energy. The 
test does not therefore directly confirm fig. 3 (a) but, on the other hand, this is 
the result which would be obtained in the presence of strongly converted 15 kev 
quanta. Removal of the absorber in front of the B-ray counter (total window 
thickness 1-7 mg/cm?) showed that the number of By coincidences per 8 remained 
almost constant (1:18 x10-? and 1:14 10-*), although the f-count increased 
considerably. From this we deduce that the %'Sm f-particles are likewise in 
coincidence with soft y-quanta (19-3 kev) or with the L x-rays arising in their 
internal conversion. ‘T’his affords direct support to the scheme of fig. 3 (b). 
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ABSTRACT. "The energies of the lowest states of ®*Li have been calculated from the 
standpoint of the shell model, using recent data for the proton—neutron forces. It is con- 
cluded that it is necessary to take into account the tensor force in discussing the lowest states 
of this and other odd—odd nuclei. 


consisting of an alpha-particle core, and an ‘orbital’ neutron and 

proton in pg states. ‘These can combine to give four possible ground 
states of spins 0, 1, 2 and 3. Since the nucleons are in the same shell, each of 
these configurations has the same energy in the potential due to the core (even 
including the coulomb force), and separation between them can only occur as 
a result of interaction between the odd nucleons themselves. ‘This separation 
must be considerable, for ®Li has but one excited state at 2-4 Mev and no others 
are known below 7 Mev. Moreover, the ground state of °He is about 3 Mev above 
that of ®Li, and so those states of ®Li corresponding to levels in °He, in particular 
those of spin 0 and 2 referred to above, must be at a similar or higher excitation 
above the ground state of °Li. It is the purpose of this note to show that the 
tensor interaction between nucleons must be taken into account to explain the 
observed spin 1 of the ®Li ground state and the large separations between the four 
states of this multiplet. 

The interaction 


eae to the j-coupling shell model *Li must be regarded as 


—kr 1) 2) 
V0) =Vo— {1 + Jel. 0—1) +p (AE ts) _ (10.9) 


Nie 


was assumed between two nucleons, and two cases were considered. 
A. V,=81:04x10-%mevcm, k=0-8580x10%cm™4, g=0-1602, y=0. 


These are given as best values for proton—proton scattering by Jackson and 
Blatt (1950). These authors do not consider the tensor interaction, but as this 
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has no effect in the singlet state it is reasonable to suppose that the value of , 
which is deduced from these data, is correct in any case. 


B. V )=46-31 x10-%mevem, k=0:5755x10%cm, g=0-:107, y=0-73. 


These are given as permissible values from the deuteron data by Hu and 


Massey (1949). Qualitatively these authors expect that to assume the larger 

valued k given above requires that V, should be smaller and y much larger, 

but no exact figures are given. After fitting the core potential to be consistent 

with the data in first approximation the additional binding energy, due to the 

interaction, for each of the four states was calculated to first order approximation. 
The results are given in the table below: 


Case A 
State of spin 3 = 2 1 0 
Contribution due to V) (Mev) 15-760 8-833 8-152 9-707 
Contribution due to Vog (Mev) 0 —1-120 +0-318 +0-606 
Total binding (Mev) 15-760 7713 9-707 10-313 
Case B 
State of spin 3 2; 1 0 
Contribution due to V9 (Mev) 3-875 3-102 2°636 2:247 
Contribution due to V 9g (Mev) 0 —0-276 +0-161 +0-484 
Contribution due to Vyy (Mev) —0:779 —0-472 +1-123 —(Q-185 
Total binding (Mev) 3-096 2-354 3-920 2°546 


It will be observed that the only effect tending to cause the state of spin 1 to 
be the lowest by a considerable margin is the tensor force (contribution of Voy), 
and that this effect is large in proportion to the size of y. 

For the reasons previously stated it is almost certain that an even larger 
separation may be obtained by choosing a value of k more in line with the 
proton—proton scattering data. It should be remembered, however, that Inglis 
(1937) obtained a ground state of spin 1 without using tensor forces, by employing, 
essentially, an L—S coupling scheme. This is inconsistent with the large 
departure from L—S coupling implied by modern shell model theories. The 
purpose of this note is fulfilled, however, by pointing out that, if the shell model 
is valid, it is necessary to take into account the tensor force in any consideration 
of the spins of odd—odd nuclei. 

My thanks are due to Mr. J. M. C. Scott for suggestions and discussion. 


REFERENCES 
Hu, Ts1-MInG, and Massry, H. S. W., 1949, Proc. Roy. Soc., 196, 135. 
INcuIs, D. R., 1937, Phys. Rev., 51, 531. 
Jackson, J. D., and Biart, J. M., 1950, Rev. Mod. Phys., 22, 77. 


663 


The Photo-Ionization Cross Section of Methane 


By A. DALGARNO 
Department of Applied Mathematics, Queen’s University, Belfast 


Communicated by D. R. Bates; MS. received 14th March 1952 


ABSTRACT. The photo-ionization cross section of methane is calculated using 
self-consistent wave functions. ‘The value at the spectral head (105 850 cm™) is found 
to be 9-4 10-17 cm?. An approximate general formula is employed to obtain results in 
a wave number range up to 200 000 cm~ and the associated oscillator strength is computed. 
It is pointed out that exchange effects probably reduce the cross section. With the aid 
of the general formula an estimate is also made of the photo-ionization cross section of argon. 


§1. INTRODUCTION 

HE method of self-consistent field, in which each electron is regarded as 
occupying a stationary state in the field of the nucleus and of the other 
electrons, has achieved considerable successes as far as atoms are 
concerned, but its extension to molecules is a matter of great difficulty because 
of the lack of spherical symmetry. Buckingham, Massey and ‘Tibbs (1941) 
have shown, however, that in certain cases it is legitimate to average the nuclear 
distribution over all orientations about a central nucleus and so obtain a 
spherically symmetric potential field; thus they point out that for systems which 
possess a high degree of symmetry, such as the tetrahedral molecules CH,, CCl, 
and SiH,, the distributions of the electrons obtained using the averaged nuclear 
distribution should not differ greatly from the actual distribution. The 
self-consistent wave functions (without exchange) of methane, derived after 
making this approximation, were used by them to evaluate the diamagnetic 
susceptibility, the polarizability and the van der Waals force between two 
molecules. The elastic scattering of slow electrons by the calculated potential 
was also studied. From the results of these investigations, Buckingham, Massey 
and Tibbs concluded that the treatment for methane is only slightly less 

satisfactory than that for atoms. 

Optical transition probabilities are sensitive to the detailed nature of the 
wave functions and hence provide a useful criterion for their accuracy. In this 
paper, the photo-ionization cross section of methane is evaluated using the 
self-consistent wave functions of Buckingham, Massey and ‘Tibbs. It is 
understood that measurements of this cross section are being made by 
Professor R. W. Ditchburn of the University of Reading and by Professor G. L. 
Weissler of the University of Southern California, so that data should soon be 
available for comparison with the theoretical result. As Professor Ditchburn 
is also measuring the photo-ionization cross section of argon a preliminary 
estimation of this is made. 

§2. THEORY 
If y,; and yp represent respectively the initial and final wave functions, 
suitably normalized, the photo-ionization cross section is given by (cf. Bates 1946) 
4m2e2 2 
Semrnes Tee | x On)xedr meer) 
A] 


a= 
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where e, m, c and h have the customary meanings, @; is the weight of the initial 
state, v is the frequency of the incident radiation and r, the position vector of the 4 
jth electron; ©, denotes a summation over the initial states and 2, a summation 
over the final states. The averaging of the proton distribution effectively replaces — 
the four protons of methane by a uniform shell about the carbon. nucleus. | 
Consequently the ground state may be represented by the configuration 
(1s)2, (2s)?, (2p)® and is closely similar to neon, except that the potential field 
of the nuclei is —~4—(6/r) or-—(10/r) according as 7, the distance from the 
centre of symmetry, is less than or greater than the radius of the proton shell. 
This similarity enables (1) to be reduced as in the atomic case (cf. Bates 1946b) 
giving the formula 


d,=1‘71(1+¢) @p(@ (6)? F208, (ey 10 em (2) 
in which J is the ionization potential and « the energy of the ejected electron in 


rydbergs; @p is a factor, close to unity, allowing for the distortion of the passive 
electron wave functions by the transition; and 


®ald=|f rP(n,1\7)P(e, 1+ 1)|7) dr|, 


where + !P(n, /|r) is the radial part of the initial bound wave function, normalized 
to unity (m and / denoting the principal and azimuthal quantum numbers) and 
r1P(e,/+1]|r) are the radial parts of the final free wave functions, normalized 
so that P(c,/+1|r) has an asymptotic amplitude of 1/e*. ‘The transitions 
which must be considered are of course 2p—>es and 2p—ed. 


2.1. Electron Wave Functions 


The radial equations for the ejected electron are 


P'+ {v= — te} P=0 


where in the usual self-consistent field notation (Hartree 1946-47) 
2 
Vq= = {6-+4r—5Y, (2p. 2p | r) —2% Y, (ns, asir)h, PL 
n=1 


2 
V(r)= 5 10-5Y, (2p. 2p Ir) —2 a y, (ws,ns |r), yA 


the radius of the proton shell being 2 atomic units. The Yo’s depend upon the 
wave functions of singly ionized methane. ‘These were obtained from the 
self-consistent wave functions of neutral methane and neutral and singly ionized 
neon by the method described in the following paragraph. 

A scaling factor was applied to the 2p wave function of methane, chosen so 
that the magnitude of the maximum of this scaled function equalled the magnitude 
of the maximum of the 2p neon wave function. A unique correspondence could 
then be established between the radial distance at which the scaled function 
attained any particular value and the radial distance at which the neon wave 
function attained the same value. ‘This correspondence was used to derive, 
from the 2p wave function of ionized neon, a wave function which, after 
application of the inverse of the scaling factor, should be a good approximation 
to the 2p wave function of ionized methane. ‘The 1s and 2s wave functions were 
found in a similar way. It is unlikely that the procedure introduces any 
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significant error; for the wave functions are required mainly for the calculation 
of the potential field in which the ejected electron moves and as will be seen later 
the cross section of methane is insensitive to the details of this potential field. 
The radial equations of the ejected electron were solved by numerical 
integration for «=0, that is at the spectral head. Normalization was effected 


Wave functions for methane. Curve a: rP(2p|7r); curve 6: P(0,s|7); curve c: P(0,d|7r). 


using the procedure suggested by B. Stromgren and described by Bates and 
Seaton (1949). The figure illustrates the form of P(0,s|7) and P(0,d|r) and 
includes for reference the function rP(2p|r) of Buckingham, Massey and Tibbs. 
It may be remarked that owing to the dominance of the centrifugal term 6/7? 
P(0,d|7) approximates closely to the corresponding hydrogenic function. 


§3. RESULTS AND DISCUSSION 
3.1. The Spectral Head Cross Section 


The wave functions described in the previous section were used to compute 
@,and@®,,,. ‘Taking the ionization potential of methane as 0-97 rydberg (Lossing, 
Tickner and Bryce 1951) and substituting in (2), the value of a, at the spectral 
head was thus obtained. The results are: @p=0-98, ®_1(0)? =4-6 atomic units, 
®,,,(0)? =26-°8 atomic units, a, (spectral head) = 9-4 x 10-1? cm?. 

Each of the radial integrals ®,, , may be divided into two portions of opposing 
sign. Bates (1947) has shown that a convenient measure of the degree of 
cancellation is provided by the quantity 1— D, where D is the ratio of the smaller 
to the larger of the two portions. For methane, 1 —D for ®,_, is 0-65 and for ®,,, 
it is 0-98. The degree of cancellation in the greater of the two integrals is 
therefore slight and there appears no possibility of gross error in the value derived 
from the cross section (cf. Bates and Seaton 1949). However, it is apparent 
from the figure that a considerable part of the integrals ®,,, comes from large 
radial distances where the discrete function is ill-determined.* As a 
consequence high accuracy cannot be expected. 


*In calculations which involve a single central potential, no advantage is gained by 
Chandrasekhar’s device (1945) of reducing the contribution® from large radial distances by 
transforming the transition matrix element from the dipole length to the dipole velocity form, since 
the two forms are exactly equivalent. In the calculations carried out here the difference between 


the potentials is trivial. 
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An obvious refinement would be to make allowance for exchange. Since | 
the contribution of the transition 2p->es to the cross section is small compared e| 
to that of the transition 2p->ed, we need discuss only the latter. Seaton (1951) = | 
has found that the Fock terms scarcely alter the «d wave functions of neon owing 
to the control exerted by the large centrifugal potential, and what little effect 
they have is repulsive so that the exchange solution is actually rather closer to 
the corresponding hydrogenic function than is the non-exchange solution, 
In contrast the 2p-wave function is made appreciably more compact. On referring 
to the figure, it may hence be seen that the value derived for the methane 
photo-ionization cross section is probably an overestimate. 


3.2. Variation of the Cross Section with Frequency 


Bates (1946b) has given a general formula for photo-ionization cross sections 
based on the assumption that the free electron wave functions may be represented 
by the corresponding hydrogenic wave functions. It has been indicated above — 
that this assumption is satisfactory in the case of P(0,d|7r); and, indeed, the 
value obtained for ®,,,(0) using the general formula agrees exactly with the 
value obtained directly. For small values of ¢, the approximation should still 
be satisfactory. As ®, ,(c) is relatively unimportant it was judged sufficient to 
apply the general formula multiplied by an external correcting factor chosen so 
as to make the value of the cross section at the spectral head equal to that given 
by the detailed computations. ‘The results finally obtained are displayed in the 
table. For convenience the wave number of the incident radiation, rather than e, 
is used as the independent variable. 


“ 


Wave number of incident 


radiation v/¢ (cm7). 105850* 120000 140000 160000 180000 200000 
Photo-ionization cross 
section a, (10-'’cm?). 9-4 7:7 5-7 4:2 3-3 2-4 


* spectral head 


For an N-electron system, the fsum rule of Kuhn and Thomas states that 
2 ,,f(t,k) =N where f is the oscillator strength of the transition i>k and the 
summation includes an integration over the continua for which the oscillator 
strength and the cross section are related by the formula df=1-24 x 10!’ a,de 
with a, in cm? and « in rydbergs. From the results given above the total 
oscillator strength associated with the 2p+es and 2p—ed continua was found 
to be 6-3. ‘The fact that this exceeds the number of 2p electrons is in no way 
anomalous since the discrete transitions 2p+2s and 2p-+1s, whose oscillator 
strengths must be included in the summation (in spite of the fact that the inner 
shells are filled), give negative contributions. 


3.3. Photo-lonization Cross Section of Argon 


‘The photo-ionization cross section of argon at the spectral head (127 110 cem~) 
was evaluated using the approximate formula mentioned earlier and the 
self-consistent wave functions with inclusion of exchange given by Hartree and 
Hartree (1938). ‘The value a,,=3 x 10-17 cm? was obtained. 
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LETTERS-TO THE EDITOR 


The Cupric Ion in a Trigonal Crystalline Electric Field 


The paramagnetic resonance spectra of single crystals of cupric salts have been 
extensively studied. The results are expressed as parameters in the usual spin Hamiltonian, 
with S=4: 

EA = Bg, HS, + Bg (H, Sy => H,Sy) ae Sele a5 BS Ly mis Sly). 


Here we have assumed axial symmetry, and have omitted the terms due to the nuclear 
electric quadrupole moment of copper, and the term due to the direct effect of the external 
magnetic field on the nucleus; these have no effect on the point which we wish to make 
in this letter. 

The spectra are usually anisotropic, with g varying from 2-4 to 2-1. But in the 
case of copper fluosilicate hexahydrate CuSiF,.6H,O, diluted with the isomorphous zinc 
salt, the spectrum is found at 90° K to be nearly isotropic (Bleaney and Ingram 1950). 
The experimentally determined parameters are gy=g,=2-24, [A|=2x10° cm”, 
|B| =3=10-* cm™. Similar results have been obtained at 90° kK with copper bismuth 
double nitrate 3Cu(NO3),.2Bi(NO3),.24H,O, diluted with the isomorphous magnesium 
salt, and with copper bromate Cu(BrO;),.6H,O, diluted with the isomorphous zinc salt. 

The crystal structures of the fluosilicate, double nitrate and bromate are hexagonal, 
trigonal and cubic respectively. The structure of the bromate is similar to that of the 
alums, with four ions in the unit cell, each subject to a field of trigonal symmetry. ‘Thus 
in each of these three salts one would expect that the static crystalline electric field has 
no components of symmetry less than trigonal. ‘The consequence of this has been 
considered by Abragam and Pryce (1950). A field of trigonal symmetry does not split 
the degeneracy of the lower orbital doublet left by the cubic field; however, a general 
theorem of Jahn and Teller (1937) shows that the energy of the system will be lowered 
if the degeneracy is lifted. ‘This requires a distortion of the octahedron of water molecules 
surrounding the cupric ion, and there is a continuous set of configurations all giving the 
same reduction in energy. ‘The system resonates between these, and the experimentally 
determined parameters represent the average of the rapidly fluctuating values. 

Recent experiments have revealed that in each of these three cases there is a radical 
change in the spectrum at low temperatures, the spectrum becoming complex and 
anisotropic. In the case of the fluosilicate and the double nitrate the change occurs between 
60° K and 20°, and in the bromate between 20° and 12° x. The change from one 


668 Letters to the Editor 


spectrum to the other occurs over a rather wide temperature range, and it is possible to 
see the two different spectra simultaneously. ‘The change is reversible with temperature. 

The spectrum of the fluosilicate at 20° has been fully analysed, using 0:8 cm7 
radiation to obtain better resolution. ‘The fluosilicate has normally only one ion in the 
unit cell (Pauling 1930), but at 20° k the spectrum is that of three ions in the unit cell. 
If the trigonal axis of the crystal is considered as the space diagonal of a cube, then the 


a eS 


three ions have each axial symmetry about one of the three mutually perpendicular cube ~ 


edges. The following values are obtained for the parameters in the spin Hamiltonians : 
&y= 2°46, gy = 2-11, | A] =11 x10 cm, | B| <3x10-* cm™. These parameters are 
very similar to those obtained in the Tutton salts of copper (Bleaney, Bowers and Ingram 
1951), and a nuclear electric quadrupole effect is also visible. Thus it seems probable 
that in this case we have a field of tetragonal symmetry, corresponding to the fourfold 
symmetry associated with the edges of a cube. 


These results show that at low temperatures the rapidly fluctuating distortions — 


associated with the Jahn-Teller effect are replaced by static distortions. The nature of 
this transition is obscure, but further information may be obtained from a careful study 
of the resonance spectrum at intermediate temperatures. At present no satisfactory theory 
of the mechanism involved has been suggested. 

The authors wish to thank Professor M. H. L. Pryce for helpful discussion. The 
work on the double nitrate was carried out by Mr. R. S. Trenam, who very kindly informed 
us of his results prior to publication. 


The Clarendon Laboratory, B. BLEANEY. 


Oxford. K. D. Bowers. 


17th June 1952. 
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The M, Absorption Edges of Tm and Er and the Origin of the 
Er M, Absorption Lines 


We have recently obtained a series of photographs of the absorption spectrum of ®°T'm 
in the M, region by means of a bent crystal focusing spectrometer, using an absorbing film. 
of thulium oxide approximately 0:1 mg/cm? thick. The photographs show no sign of the 
Myy absorption edge, while in the My region a narrow absorption line only is apparent. 
This result may be compared directly with Lindberg’s result for the two neighbouring’ 
elements ®*Er and “Yb. Lindberg (1931) found no Myy edges for either of these elements, 
a single My absorption line for Yb, and a close pair of absorption lines in the My region for 
Er. Density curves of the Tm absorption line taken with a recording microphotometer, 


however, show that the line is definitely asymmetrical, having a steeper profile on the long’ 


wavelength side, which suggests that this is actually the My edge. A tracing of a typical 
density curve is shown in curve (a), where A is the point of minimum density (i.e. maximum 
absorption). Closer examination of these curves reveals that the long-wave side changes. 
its gradient rather abruptly by a small but noticeable amount at the point B. This effect 
appears so consistently in approximately the same position on the microphotometer record- 
ings of all the photographs taken that it must be considered to be real, and indeed suggests 
that we are dealing with two close unresolved absorption lines of unequal intensity of the 
type indicated by the dotted curves in (6). Our measurements show that these two lines 
would be separated by approximately 10 x.u. which would fall between the measured’ 
separation of 17 x.u. found by Lindberg for the two Er lines, and the zero separation for the 
Yb line. Assuming this analysis, the longer wavelength would be ascribed to a transition 


from the My level to an unoccupied 4f level, and the long wavelength profile of the shorter 


wavelength line would represent the My edge modified by the presence of the absorption 
line. ‘The conventional wavelength of the My edge would then be that of the point K in 
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(5), the point of inflection of the reconstructed My edge. Using Ka lines of Zn, Mn, Al, 
Ti and Ni in high order spectra as reference lines we estimate the wavelength of this point 
as 8469-5 + 3 x.u. and of the point A as 8461:-4+2x.u. The value of v/R for the edge is then 
found to be 107:-59+0-07. A computed value of 107-65 is obtained from the wavelengths 
of the Ly; edge found by Coster (1923) and of La, by Coster, Nishina and Werner (1923) 
(these are the wavelengths given in the tables of Cauchois and Hulubei 1947), so that the 
discrepancy is effectively zero, in contrast to the small but definitely negative values found 
in the case of other rare earth elements (Stewardson and Zandy 1950, Lee and Stewardson 
1951). A modified Moseley plot of all relevant frequencies (Ly, La, and My) for the rare 
earths shows that the accepted value of TmL«, is rather low compared with that of neigh- 
bouring elements and the low value of the L-My discrepancy might well arise from this fact. 
The wavelength of the assumed absorption line accompanying the edge cannot be arrived at 
with any great reliability but will obviously differ from that of the edge by a few x.u. only. 


(a) 
B 
A 
<j) v— 
(c) 
B K 
A 
ay Vv —> 


This result for Tm suggests that the results of Lee and Stewardson for the My edge of 
Er can be reconciled with those of Lindberg if it be assumed that Lindberg’s two lines have 
the significance suggested here for the supposed ‘I'm doublet, i.e. the longer wavelength 
corresponding to a transition to a 4f level, and the edge being given by some wavelength 
between the two lines at approximately 8813 x.u., instead of 8829 x.u. as found by Lee and 
Stewardson for a very thick absorption film. The hypothetical curves (c) (thin film) and 
(d) (thick film) show how this difference could arise, K’ being taken for the edge in the 
thick film, whereas the true edge is at K. Such progressive changes in appearance of an 
edge accompanied by a line absorption have been observed by Rule (1945) for °*Sm, and by 
Zandy, working in this laboratory, for various rare earths. ‘The modified value of the 
ErMy wavelength suggested here would reduce the L—My discrepancy from —15 ev as 
previously found by Lee and Stewardson to the still large but more reasonable value of 
—12 ev approximately. Further experiments on the My and Mry edges of these two 
elements are in progress. 


Physics Department, P. A. LEE. 
University College of Leicester. E. A. STEWARDSON. 
16th May 1952. , J. E. WILson. 
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The Nature of the Neutral ~-Meson 


It has been shown recently (Heckrotte, Henrich and Lepore 1952, Francis and 
Marshak 1952) that the 7°-meson would be expected to be produced much more readily ~ 
in elastic photoproduction from deuterium, y+-D=7°+-D, for asymmetrical meson theory 
than if the meson were coupled with the same sign to neutron and proton. It is the 
purpose of this note to consider the possibility of determining the relative sign of coupling, 
or indirectly of the parity, of the 7°-meson from differences to be expected in the | 
reaction n+p=7°+D. 

Table 1 shows the selection rules operative in this reaction for initial relative angular 
momenta 0, % or 2% and for 7°-meson produced in s, p and d states. The differences that 
occur for the Yukawa neutral meson, where gy=gp, and the Kemmer neutral meson, 
&y—=—gp are due to the antisymmetry properties of the two fermion wave functions 


y 
d 
= 
5 
7 


Table 1. Allowed Transitions for the Reaction p+n=7°+D 


Meson 


a Scalar Pseudoscalar 
P Y K Y K 


system 


aS Yes 
aS) Yes 


sie Yes* Yes 4 
sp Yes Yes* Yes. 


1D Yes* 
8—D Yes Yes* 


* Forbidden if the fourth component of P.S.(P.V.C.) is inoperative. 
NG: Ex Sp; See Ex aeep: 


Table 2. Excitation Values and Angular Distributions for Allowed Transitions 


Meson 


Scalar Pseudoscalar 
1S — Stee op) 
4S) — T3/2; cos? 0:(Y, p) 
ae —_ PLAIN CYSIS) i 
=P T3/2. cos? 6 (K, p) T22: 1 (K,-s)t 
1p Tei cos™ O(a p) 


sD T°? ; ~cos! 6 (Y, d) 722: cos* OY, p) 
* Small cos? 0 term. 
t Also T?/2: cos? 6 term from s—d interference. 


combined with the necessity of an isotopic singlet initial state for the Yukawa neutral 
meson, while since 73 @)|8(7))) =—(—1)‘|(z)9), the Kemmer neutral meson is only produced 
from isotopic triplet initial n—p states. Qualitative differences are evident from 
consideration of this table and are verified by first-order perturbation calculations in many 
particle single time formalism. ‘lable 2 gives the excitation function and expected angular 
distribution in the centre-of-mass system. Considering scalar mesons it is seen that the 
meson has to be created into d states if gy=gp while p states are available for &n= —2p- 
The excitation function varies as T°!, where T is the meson kinetic energy in the 
centre-of-mass system, for Yukawa mesons, while the predicted cross section has a sharp: 
forward peak, following an approximately cos*@ distribution. For Kemmer neutral 
mesons the cross section behaves as in the corresponding reaction p+p=a+-+-D; that is. 
to say an excitation function of T°/? and a cos? @ angular law (Gunn, Power and Touschek 
1951). For pseudoscalar mesons the selection rules are less effective. When accepting 
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the ad hoc assumption that the odd part of the pseudovector coupling is inoperative, the 
neutral and symmetric mesons must then come off in d or p states. This has been suggested 
by Francis and Marshak (1952) to explain the angular distribution for 7°-meson production 
in y—-p collisions and the forward cross section in 7+ production by protons on protons 
(see also Chew, Goldberger, Steinberger and Yang 1951). It is only through a dominating 
fourth component in the interaction, for p-state mesons or from s-d interference that a 
non-isotropic distribution is favoured. This is necessary in the symmetric case to explain 
the observed behaviour in charged meson production. 

‘The magnitudes of the cross sections are readily obtained. ‘That for the Kemmer meson 
can be obtained from the charged meson results from which it differs by a factor two, 
except for the difference in rest mass of the 7+-meson and 7°-meson. At 350 mev the 
integrated cross section is approximately 1:1 g?/he x 10-8 cm?, for neutral scalar mesons 
the cross section is approximately 3-6 g?/Rc x 10~ ®°cm?. 

If other evidence, such as from the reaction y+D=7°+D, gives information concerning 
’ the isotopic nature of the 7°-meson then the above interaction will distinguish scalar from 
pseudoscalar. On the other hand if the parity is decided by other means (Marshak 1951) 
it should be possible to: obtain a decision on the relative sign of coupling. Other reactions 
with more immediate possibility of verification, involving 7° production from nuclear 
collisions, such as from protons of energy greater than 200 Mev on deuterium, and also 
consideration of relative cross sections for neutral and charged production, will show 
differences between Yukawa and Kemmer neutral mesons. 


Mathematics Department, E. A. Power. 
University College, 
London. 
8th May 1952. 
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Multiple Resonances in Photodisintegration Cross Sections 


The reactions *O(y, 4x) and C(y, 3x) have been studied by irradiating Ilford nuclear 
emulsions (200 micron E1 and C2) with bremsstrahlung, using peak energies from 
20 to 70 Mev. Some 700 *O(y, 4x) stars and 1 800 C(y, 3x) stars have been identified 
and analysed by techniques already described (Goward and Wilkins 1950), and the 


= iS 
a 03 thee 
= E 
= _S) 
=. 2 
— So 
= O2 2 2 
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3S = 
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Gamma-Ray Energy (MeV) 


results reveal multiple resonances in the dependence of the cross sections on gamma-ray 
energy. ‘These resonances suggest that the reaction should be interpreted in terms of the 
excitation of levels in a compound nucleus. 

The figure illustrates the nature of the cross sections which have been found, and 
presents the “O(y,4«) cross section in the form of a corrected histogram. The 
uncorrected histogram recorded the gamma-quantum energy FE, appropriate to each 
1’O(y, 4x) star, and was then corrected for the bremsstrahlung spectra, for ‘ escape’ of 
stars from the emulsion layer and for other minor factors. Standard deviations are 
indicated by vertical lines and the experimental resolution (fairly closely known from 


672 Letters to the Editor 

= 
experiments on the 2Cly, 3) reaction using ’Li(p, y)*Be gamma-rays) by the dotted } 
gaussian curves. At least three maxima, or resonances, are established, at gamma-ray 
energies of about 22:6, 25:8 and 29-5 mev respectively. Comparison of the first two with — : 
the superposed dotted curves suggests real widths appreciably less than 1 Mev. i 

Making the assumption that the true cross section consists of a continuum of some 
form with relatively narrow resonances (<1 Mev) superposed thereon, it is clear that the — 
present techniques cannot decide the relative contributions of continuum and resonances 
with any accuracy, although the continuum must be small near E,=22-6 and 25-8 Mev. 
As an extreme illustration, the assumption of very narrow, discrete resonances at 22:6, 
23:5, 24:5, 25-8, 27:2, 28:2 and 29-5 mev, of relative magnitudes 1:00, 0°36, 0:29, 0-65, 
0:32, 0:53 and 0:94 respectively (the values being obtained by trial curve fitting) would 
yield excellent agreement with the experimental data between 21-5 and 30:0 Mey without 
requiring any continuum. 

The main trends of the C(y, 3x) cross section have already been described (Wilkins 
and Goward 1951). Detailed study of results now available (Goward and Wilkins 1952) ~ 
establishes at least four resonances at E,,=18-4, 21-8, 24:5 and 29-4 Mev respectively; 
the widths of the second and third must be appreciably less than 1 Mev. Several further — 
resonances are strongly suggested by the results, particularly below E,=18 Mev. Again 
adopting the extreme illustration used in discussing the “O(y, 4a) results, the assumption 
of narrow, discrete resonances at 15-0, 16:1, 16-8, 17-4, 18:4, 19-5, 20-7, 21-8, 23-2 and 
24:5 mev, of relative magnitudes 0-30, 0:34, 0-28, 0:57, 1:00, 0-32, 0:08, 0:32, 0-23 and 
0:70 respectively, gives excellent agreement with experiment between 14:5 and 25 Mey. 
Resonances at 1°C excitations of 16:1, 16-7, 19:4, 20:7 and 21-8 Mev have been reported 
in other reactions (Hornyak et al. 1950, Blaser et al. 1951). Such a series of discrete © 
resonances is too drastic a simplification however, since it does not account for the known 
cross section at E,=17-6 Mev (Nabholz et al. 1952, Goward and Wilkins 1952), or for 
certain forward/back asymmetries in the alpha-particle angular distributions (Wilkins 
and Goward, to be published) which indicate overlapping of !*C states. 

A significant feature of the results for both the ’C(y, 3x) and *O(y, 4x) reactions is 
that the general trends of the cross sections show minima where the corresponding 
(y, n) cross sections have maxima. This again suggests a compound nucleus interpretation. 


Atomic Energy Research Establishment, F. K. Gowarp. 
Harwell, Didcot, Berks. J. J. WILKINS. 
23rd May 1952. 
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The Scattering of Positrons by Nuclei 


The scattering cross sections for positrons by the coulomb fields of nuclei have been 
calculated for the energy range 1-10 Mev of the incident positrons. A procedure similar 
to that used by McKinley and Feshbach (1948) for electron scattering has been applied. 
This method starts from Mott’s formula (Mott and Massey 1949) for the differential 
scattering cross section o(@) : 


k*a (0) =q">(1 —B?) FF* cosec? 40+ GG* sec? 40, — ...... (1) 


where q=a/f, «=Ze*/he=Z/137, B=v/e, kh is the momentum of the incident positron 
and F and Gare obtained by replacing a with —a in the similar quantities for electron 
scattering defined by Mott and Massey. F and G are expanded in powers of a, 
keeping terms up to a* and the coefficients in the expansion calculated. The ratios R of 
the differential cross sections so obtained for positron scattering to those obtained using 
the Rutherford formula, 

o'p(0) = (Ze®/2ene*)*{(1— B2)/B*} cosec* 40, ra (2) 


were calculated. 
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The error in R due to this procedure is of the order of at and should be reasonably small 
for values of Z less than 50 (at=0-018). ‘To obtain values which should be accurate to 
within about 2°% up to Z=80, the exact positron scattering calculations of Massey (1942) 
for Z=80 were compared with those calculated using the above expansion and the size of 
the necessary correction term thus obtained. (The extrapolation of Massey’s values to 
higher energy was carried out by plotting his differential cross sections against # for each 
scattering angle.) It was then assumed that the correction to R should be proportional 
to Z4 for Z less than 80 and in this way values of R were obtained. 

These values are given in the table. 


Values of R=o(6)/o,(8) 


@=the angle of scattering in degrees, H=energy of the incident positron in Mev, 
Z=the atomic number of the scatterer. 


6 30 45 60 80 90 100 120 135 150 180 


\ 
by 


0-903 0°8255 0-731 0°591°~ 0-519 0-448 0-314 0-232 0-168 0-115 
0:897° 0-812. 0-711 + 0:558 0-480 0-403. 0-258 0-169 0-099 0-042 
0:894 0-807 0-702 0-544 0-463 0-383 0-233 0-140 0-068 0-009 
0-893 0-806 0-700 0-541 0-460 0-379 0-229 0-134 0-063  0-0024 


0:875 0-791 0-697 0-560 0-492 0-425 0-300 0-224 0-166 0-118 
O:367— 0-707 0-675) 0-526 0:4520°379F 0-243) 0:1592-0-096 07043 
0:863 0-770 0-665 0-511 0-434 0-358 0-217 0-130 0-068 0-009 
0-862 0-769 0-663 0:507 0-430 0-354 0-212 0:124 0-058 0-0024 


10 


— 
ounN eR 


20 


— 
ounNnr 


0-854 0:769 0-676 0-543 0-478 0-414 0:296 0-222 0-167 0-121 
0-845 (0-752) 10-651) 0-507 0-436 «0-365 0-236 0-155 0-095 0-044 
0-840 0-745 0-640 0-491 0-417 0:343 0-209 0-125 0-062 0-009 
0-839 0-744 0-638 0-487 0-413 0:339 0-204 0-119 0-056 0-0025 


0:836 0-751 0-659 0-529 0-465 0-403 0:289 0:220 0-170 0°127 
0-826 0-733 0-633 0-490 0:420 .0-352 0-227 0-156 0-095 0-047 
O82 1unO 724 0-020 -47206 02399107329) 510-1998 10-119 (0-061) 0:01 
0°820 0-723 0-618 0-469 0-396 0:325 0-193 0-113 60-055 0:005 


30 


= 
oOuUN re 


40 


b= 
oOuUnNrF 


O°8235 (07/4. 80,6545 0-522 50-460 10-4007 0-291 0.224 0-175" 0-134. 
0-012 0-7 2160-024 9 0-483) 0-413) 10°346" (0:2242 0-151) 0-098" (0-051 
0-300 0-7 11 0-61000°463.980°392)10-321- 0-194 2 0-119 0.0639 10-014 
0-805 0-709 0-607 0-460 0-388 0:316 0-188 0-113 0:056 0-007 


50 


_ 
oOuUnNre 


0-813) 0-732 0-646" 0-517 710-456 (0:398) 10-294 0-229) 0-172 (0-141 
0-800 10-7108 0-610980:476 0407 0-341 0:222) 0-156) 02102) 0057 
0-793 0-700 0-601 0-455 0-385 0:314 0:190 0:120 0-067 0-019 
0-791 50-698) 10597) 0:450, 10-379. 07309) 0-184) 0-113 202060 0:013 


60 


sy 
oun re 


0°303, 40-722 = 0-635— 0-508" 0-446; <0:394" (0:297 -0:237 (0:191 90-149 
0-790 0-697 — 0-598 0-466" 0-396 0:335 0-219" 0:162--0-111 ~—0:065 
0-781 0-686 0-586 0-443 0-373 0:306 0-186 0-124 0:076 0-027 
0-780 0-683 0-583 0-436- 0:367 0:298 0-179 0:116 0-067 0-020 


— 
ounNrR 


0-795.) 0-703") 0-612" 0-487 70-427) 0:385 0-300 0:245 0:201 90-157 
0-780 0°677 0-577 0-445 0:376 0°322 -0:215 0:167 0-123 0-077 
0-770 0-663 0:560 0-420 0°353 0-290 0-181 0:131 -0:090 0-038 
0:768 0°660 0-556 0-413 0-345 0-280 0-173 0-123 0-080 0-032 


80 


—s 
oOuNnNrF 
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suggestions during the course of the present work. I also express my sincere thanks to 
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22nd May 1952. 
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Neutron Energy Measurements with a Helium 3 Filled 
Proportional Counter 
Fast neutron energies are usually measured by methods based on the scattering of the 


neutrons by hydrogen. The incident neutron energy spectrum is obtained either by 
differentiation of the complete recoil proton pulse height distribution or by measurement 


of the energies of the protons recoiling in the forward direction. Scintillation counters _ 


with organic phosphors are quite sensitive detectors which can be used for the first method, 
and recently satisfactory results have been obtained (Poole 1952). One way to accomplish 
the second method is by use of the photographic plate as detector (Nereson and Reines 1950). 
However, the grain size and sensitivity set a low energy limit of about 500 kev for this 
technique. It has, in addition, the disadvantage of being laborious and time-consuming. 
The drawback of any method designed to measure neutron spectra which is based on 
scattering lies in the conversion of the definite energy of the incident neutron into a con- 
tinuum of energies of the recoiling particle. For some considerable time a suitable 
nuclear transformation induced by neutrons has been sought, in which one energy release 
corresponds unambiguously to one neutron energy. In order to be successful the cross 
section of such a reaction should be fairly large and its variation with neutron energy should 
be smooth and free from sharp resonances. Also the energy release for zero neutron 
energy should be small, and there should be no low-lying excited states of the residual 
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nucleus leading to the emission of particles of reduced energy. Of the various possibilities 
(GAN, 1°B, ®Li, 7He), some of which have been tested experimentally (Keeping and Roberts 
1950), the most satisfactory is based on the disintegration of *He : 


gHe-+ in iH 7 7Olkevy 4 bege Meeeperceeres (1) 


The thermal neutron capture cross section of this reaction has been measured by 
Batchelor et al. (1949), 3750 barns, King and Goldstein (1949), 5000 barns, and Coon and 
Nobles (1949), 5 040 barns. The variation of cross section with neutron energy has been 
measured in the thermal energy region (King and Goldstein 1949) and in the 1 Mev energy 
region (Coon 1950). ‘The value of the cross section for 1 Mev neutrons is about 1 barn 
which is comparable with the scattering cross section of neutrons of this energy by hydrogen. 
In reaction (1) the total energy of the proton and triton is (H,+0-77) Mev, where Ey is the 
kinetic energy of the bombarding neutron. 

The energy range of this detector is limited to a maximum neutron energy of 1 Mev because 
then the neutron—*He elastic collision produces a maximum *He recoil energy of 3 Mev, 
which event cannot be distinguished from a *He slow neutron disintegration. Above 
1 Mev neutron energy the results of this method are ambiguous. However, if the neutron 
spectrum to be investigated has a definite upper energy limit E,,, then the maximum He 
recoil energy will be #/,. It follows that the energy release in the interval from 3, to 


(E+ 0°77) Mev can be assigned directly to the primary neutron spectrum in the range 
(ZE, —0-77) to Ey. 2 
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To put this scheme into operation a proportional counter containing *He and xenon to: 
increase the stopping power has been constructed. Since the wall effect is a function of 
neutron energy, it is important to make it small. The presence of xenon unfortunately 
increases the y-ray sensitivity of the counter and so increases the low energy background... 
The counter is 11 in. long and 3 in. in diameter and is filled to a pressure of 110 cm xenon 
and 5:4cm of helium containing 10% *He. This helium was supplied by the Harwell. 
Thermal Diffusion Group. 

Figure 1 shows the pulse height distribution obtained with thermal neutrons for a gas: 
amplification of 4. The width at half height is 55 kev. ‘The counter was tested with. 
approximately monokinetic fast neutrons from the 7Li(p, n)’Be reaction. Protons were 
accelerated by the Harwell Van de Graaff machine on to a lithium fluoride target. Curve a 
of fig. 2 shows a complete pulse height distribution including y-ray background, thermal 
peak (due to thermal neutrons from wall scattering etc.) and a small fast neutron peak for 
neutrons of energy 560 kev. Curve 6 is drawn with a much increased ordinate scale to show 
the fast neutron peak more clearly. Part of the width of 150 kev is due to spread in the 
incident neutron energy introduced by target thickness and geometry (50 kev and 30 kev 
respectively). Similar response is obtained with neutrons of energies up to 800 kev, above 
which tests have not been made. With 350 kev neutrons, the peak is not completely 
resolved from the background. 

I wish to thank Mr. P. Tunnicliffe for help at the start of this work, Dr. W. D. Allen for 
use of the Van de Graaff generator, and Dr. E. Bretscher for helpful guidance and discussion. 
I also thank the Director, Atomic Energy Research Establishment, for permission to publish 
this work. 


Atomic Energy Research Establishment, R. BATCHELOR. 
Harwell, Didcot, Berks. 
19th May 1952. 
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The Decay of the 2°4 Mev Level in °Be 


The 2-4 mevy level in °Be has been studied through the reactions *Be(p, p’)*Be* (Davis 
and Hafner 1948, Rhoderick 1950, and Browne et al. 1951) and ™B(d, «)®Be* (Buechner e¢ al. 
1951). Neutron emission from this level is energetically possible as an alternative to gamma 
radiation, and we have found that neutrons are indeed emitted. 

A target of separated 4B was bombarded with 425 kev deuterons. Alpha-particles 
emitted at 90° to the beam were detected by a thin anthracene screen and E.M.I. photo- 
multiplier. A liquid phosphor (1-4 litres of a terphenyl solution in toluene) was used to 
detect the gamma-radiation or neutrons. It was mounted at a distance of about 50 cm 
from the target in a plane normal to the deuteron beam, and could be set at various angles to 
the first detector. The delay spectrum, i.e. the coincidence rate as a function of the delay 
between the pulses from the detectors (James and 'T'reacy 1951), was studied using a circuit 
of 8 mysec resolving time, similar to that described by Curling and Newton (1950) and 
modified later by Newton (1951). 

If the reaction “B(d, ~)*Be* is followed by neutron emission to the ground state of 
8Be, recoil effects produce a variation in the neutron time-of-flight between 56 musec per 
metre at 180° to the alpha-particle direction and 187 musec per metre at 0°. For the same 
reason the rate of alpha—neutron coincidences should rise rapidly in going to 180°. 

Such a sharp increase in the coincidence rate was indeed found in a preliminary study 
using a radiation detector placed close to the target. A delay spectrum with the radiation 
phosphor 40 cm away from the target and at an angle of 70° to the alpha detector gave only a 

single peak, at a delay of 8 msec in the radiation detector, and we must ascribe it to 
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particle-gamma coincidences. (The 8 mpsec shift is to be expected if the gamma-ray 
pulses are larger than those from the other detector.) 

The delay spectrum was next measured at 180° with the liquid scintillator 50 cm away 
from the target, giving curve 1 of the figure. Here a delayed group is clearly seen. The 
background (mainly from reaction *O(d, p)!”70*(y)!7O) was obtained with a foil inserted 


to cut off the alpha-particles from 1B(d, «)*Be*, and is shown in curve 2. The position of | 


the particle-gamma coincidence group served to confirm the zero point of the delay spectrum. 
‘The two spectra were obtained in one session and under the same conditions ; the average 
true-to-chance ratio for curve 1 was 1 : 1 and the true counting rate nearly one per second. 
The calculated delay of the alpha—neutron group is 28 mysec. The observed value is 
27+5mypsec. From these spectra we estimate the proportion of alpha~gamma coincidences 
at 180° to be 2-5+2-9% of the total coincidence rate. 


8 
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Arrow indicates expected delay of neutron group. 


Assuming the alpha-gamma contribution to be small we studied, in steps of 10°, the 
correlation of the alpha-particles and neutrons with a radiation detector close to the target. 
(The background was determined as before and subtracted.) The rise towards 180° was 
confirmed, the intensity at this angle being about ten times that at 70°. This anisotropy, 
however, almost completely disappeared when the correction for the centre-of-mass effect 
was applied (disregarding the change of sensitivity of detection with neutron energy). 

The small width of this level (< 3 kev, Browne et al. 1951) could be due to the angular 
momentum barrier for the low-energy neutrons emitted if a large spin change (2 or 3) is 


involved in the transition. Such small neutron widths have been observed (Bockelman 
et al. 1951). 


We wish to thank Sir John Cockcroft and the Atomic Energy Research Establishment 
for providing the separated 4B targets. We are much indebted to Professor W. E. Burcham 
for suggesting the experiment. One of us (G.A.D.) is grateful to the Ceylon Government 


for a scholarship, and the other (J.O.N.) to the Department of Scientific and Industrial 
Research for a grant. ; 


Cavendish Laboratory, 
Cambridge. 
26th May 1952. 
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Internal Conversion Electrons from the ”°U Decay 


In the artificially produced (4n-+-1)-series (Seaborg, Gofman and Stoughton 1947) 
previous investigations of y- and x-rays indicate the existence of at least two excited states. 
of the thorium resulting from the process 7°°U -+ 22°Th (Studier 1947, West 1951). Using 
8-sensitive nuclear emulsions, a preliminary study has been made of the internal conversion 
electrons produced in this transition. Ilford G5 emulsions, thickness 200 1, were impreg-. 
nated with uranium 233 in solution as the complex ammonium uranium citrate, adjusted to 
pH=7. ‘The exposure time was 120 hours. 

Out of 6 012 alpha-tracks, 548+ 60 (i.e. 9:1+ 1°) were associated with electron tracks. 
The uncertainty arises mainly from the difficulty of distinguishing the very short tracks. 
The grain count distribution of the electron tracks is shown in the diagram. ‘Two groups of 
tracks are clearly resolved, the ratio of their intensities being 3-6 to 1. 
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Grain count distribution for 548 electron tracks. 


Using the relation between total grain count and energy established by Zajac and Ross 
(1949), the first group corresponds to an energy of about 20 kev, the second to one of about 
33 kev. This is in agreement with the existence of an excited state of energy about 40 kev, 
assuming that the two groups mainly consist of L and (M-++N) internal conversion electrons 
respectively. The interpretation is consistent with the observations by West (1951) of a 
y-ray of 43+0-5 kev. There were 53 cases in which two electron tracks undoubtedly 
came from the origin of the same a-track. 44 of these could be interpreted as 43 kev conver- 
sion electrons associated with L Auger electrons. Eight events agree with the assumption 
that the weak 56+ 0-5 kev y-line reported by West (1951) occurs in cascade with the 43 kev 
transition. The remaining double event consists of a 10 kev and an 82 kev electron track,, 
which may be interpreted as due to an L Auger electron and an internal conversion electron 
arising from a direct transition with an energy of 43-56 kev. ‘This interpretation is sup- 
ported by the existence of three single tracks attributable to L or M conversion of 99 kev. 
An alternative interpretation of this double event in terms of K Auger effect is somewhat 
less plausible. The intensities observed in this work are consistent with those published by 
West (1951). If then the 43 kev and 56 kev y-rays are in cascade and if the 56 kev y-ray is 
weak compared with the 43 kev ray both in emission and in internal conversion, then the 
first excited state is at 43 kev and the second at 99 kev. A check on the density of the 
background electron tracks showed that not more than four of all «-tracks might be wrongly 
interpreted as connected with an electron track because of chance superposition of the two 
events. The probability of a double electron event being spurious is negligible. 

By using total grain counts instead of range measurements the time required for obser- 
vation is much reduced. The use of grain counts as a measure of energy was justified by 
measuring both the three-dimensional range and the total number of grains for a limited 
number of tracks. Consistent energy values were obtained from the two measurements. 
An extended investigation of the electron energy distribution is contemplated. 
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REVIEWS OF BOOKS 


Advanced Five-Figure Mathematical Tables, by C. AtTtwoop. Pp. v+69. 
(London: Macmillan, 1951.) 4s. 6d. 


The meaning of the word ‘advanced’ in the title of this collection of tables is somewhat 
obscure, and is not greatly clarified by the explanation in the Introduction that it is a 
continuation of another collection issued by the same author under the title Practical 
Five-Figure Mathematical Tables. But what’s in a name? It is, in fact, a handy set. 
On a good-sized page, we find the trigonometric functions to five significant (not decimal) 
figures, tabulated both in degrees and decimals, and in radians, and with a good deal of 
trouble taken to make linear interpolation possible. There are also tables of hyperbolic 
and exponential functions and of the gamma function, which here is very sensibly treated 
as the factorial function and written x!. Short tables of the complete elliptic integrals 
-of the first and second kind take the square of the modulus as argument, an innovation 
which experience may well show to be advantageous. ‘There are also short tables of the 
ordinate and of the area under the Gaussian curve; those interested in statistics would 
probably have wished for Student’s t-distribution at least as often as they wish for the 
normal error curve. 

The sources of the tables are not mentioned, but the author makes clear the great 
attention which he has given to the arrangement of them, with a view to reducing the 
occasions when linear interpolation is insufficient, and the care which he has given to 
avoiding numerical errors arising in the setting-up. J. Eee 


Ferromagnetism (Bell Laboratories Series), by R. M. Bozortu. Pp. xvii+968. 
(New York: Van Nostrand; London: Macmillan, 1951.) £6 10s. 


The subject of ferromagnetism nowadays pervades a wide field of pure and applied 
science. One of the oldest of the sciences, it is as important to the metallurgist, the chemist 
and the engineer, as to the physicist, whether an instrument designer or a student of the 
quantum mechanics of the solid state. Such an encyclopaedic subject needs an 
encyclopaedia to do it justice, and Dr. Bozorth’s book succeeds in some measure to satisfy 
that need. ‘The author is well qualified to carry out this monumental task. He is amember 
of the Bell Laboratories, whence so much of our recent knowledge of the nature and 
applicability of ferromagnetism has sprung, and his name is known wherever ferro- 
magnetism is studied. 

The book is written with emphasis on practice rather than on theory, although topics 
like domain theory and the molecular field concept also find their place. This accentuated 
practical approach may easily disguise the importance of phenomena which are of 
considerable interest in theoretical, and hence eventually, it is hoped, practical investigations. 
‘Thus, in discussing Ni-Co alloys we find the statement that ‘“‘ none of these alloys have 


~ 
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any important industrial uses. Even the magnetic properties have no interesting 
peculiarities . . .”’. To the theoretical student this very lack of interest, smooth curves 
rather than spectacular peaks, may be a substantial boon. Elsewhere in the book, 
especially where the author’s own particular specialities are described, the approach is 
much happier. In fact, some theoretical sections, such as that on paramagnetism, could 
well have been omitted in favour of even more detailed experimental results, such as the 
important magnetic measurements by Sucksmith on the temperature variation of the 
magnetization of Fe—Ni-Al alloys, basic to modern permanent magnet materials. Many 
of the experimental sections, such as those on magnetic annealing of the permalloys and 
on the AF effect, are excellent and should prove to be most stimulating to teachers and 
investigators. The sections on powder patterns include photographs of the beautiful 
“tree patterns’. ‘These and other sections are frequently supplemented by useful practical 
hints, although it may be questioned whether they would be equally useful to one who 
was not a member of the Bell Laboratories. 

The book opens with two introductory chapters on basic concepts and phenomena, 
followed by eight chapters on specific magnetic materials. Pure metals as well as many 
important (in the sense indicated above) binary and ternary alloys are described, those 
based on Fe-Si and Fe—Ni being naturally of particular interest as they provide most of 
the magnetically soft materials of technology. Permanent magnet alloys are discussed 
in a long chapter, where slight digressions are sometimes made : thus, an account of the 
systems Co—Pt and Co-Pd is followed by a brief description of the corresponding nickel 
alloys, which have never been found to have permanent magnet properties. Here and 
elsewhere there is apparent an endeavour to include everything. 

The next nine chapters are collected under the perhaps unfortunate title ‘ Magnetic 
phenomena and theories’. A first thought would lead one to think, quite erroneously, that 
these are entirely divorced from the previously considered properties of materials. 
Actually, many important and interesting topics are dealt with, such as magnetic anisotropy 
of single crystals, magnetostriction and temperature and time effects, all of which have a 
profound influence on the technologically important properties of the materials. The 
final chapter in this part is entitled ‘ Special problems in domain theory’. Many of these 
problems are, however, not ‘ special’ and might well have been included in the earlier 
parts of the book. Thus, a discussion of domain geometry could well have accompanied the 
section on powder patterns. The last chapter in the book deals with magnetic measurements, 
and in the main appendices are tables of the properties of soft and hard magnetic materials. 
The wide range of the characteristic properties of magnetic materials (for example, coercive 
force values ranging from 2 x 10~*, supermalloy, to 6 x 10%, silmanal) becomes particularly 
astonishing when these properties are presented tabularly. 

Many of the disadvantages of the occasionally unsatisfactory arrangement of the subject 
matter and of incomplete cross references are offset by a detailed subject index. Here 
alloy systems are well arranged in alphabetical order and, since among the 800 or more figures 
there are very many binary and ternary phase diagrams, parts of the book constitute a 
miniature Hansen. The diagrams are consistently clear, but the captions are sometimes 
not sufficiently detailed. ‘The author index refers the reader to the list of references 
containing some 1700 publications between 1842 and 1951. ‘The first of these is by Joule 
on ‘ A new class of magnetic forces’ and the last by Weil on ferrites, in which magnetic 
forces of a new kind are again important. ‘The wheel seems to have come full circle, 
but the reader of this book, if he can afford it at the exorbitant price, will have gained the 
impression that the wheel is capable of many more revolutions, leading through domains 
yet unexplored and fields yet untrodden. E, P. WOHLFARTH. 


The Measurement of Radio Isotopes, by D. 'Taytor. Pp. viti+-118. (London: 
Methuen, 1951.) 6s. 6d. 


This little volume will be welcome to the rapidly growing number of people (engineers, 
biologists, doctors, etc.) who have to employ nuclear counting techniques without being 
nuclear physicists. Their job has become much easier now that nearly all the equipment 
required is manufactured, and [Dr. Taylor’s book quite rightly refrains from describing 
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any of the techniques for making counters and electronic circuits which constituted much 
of the art of nuclear physics in earlier days. Instead he briefly describes the physical 
principles involved in some of the commercial instruments (counters, scalers, counting 
rate meters, etc.) and gives common sense instructions for their use and testing. The 
book also contains in some detail the simple mathematics involved in the interpretation 
and planning of radioactive measurements; the elements of differential calculus are assumed — 
and employed to derive formulae to allow for the effects of radioactive decay, self-absorption, 
solid angle, resolving time and statistical fluctuations. Most of the book refers to the use 
of Geiger-Miiller counters, but proportional counters and scintillation counters are 
briefly described, and there is a chapter on monitors used for detecting contaminations 
and radiations dangerous to health. 

The style of the book is pleasant and readable,-and the mastery of the subject by the 
author is obvious throughout, occasionally from his use of terms or concepts which have not 
previously been explained. For this reason, and because of the stress on mathematical 
treatment, this is no book for bedside reading. But for the serious student it should not be 
difficult, and there are many numerical examples for the lovers of brass tacks. 

O. R. FRISCH. 


ve 


Facts from Figures, by M. J. Moronry. Pp. 472. (Harmondsworth: Penguin, 
195.1) 5s: 


Neither in scope nor price can this book be classed among the usual run of Pelican 
primers. ‘To anyone who had not already some knowledge or interest in. statistics its size 
might be discouraging. But to those who realize the possible value of statistical methods, 
and want a general introductory survey of the subject, this book should be most valuable. 

Starting from the simpler ideas of probability, the earlier part of the book leads up to a 
descriptive account of the more important statistical distributions. Subsequent chapters 
show how these ideas lead to various applications, ranging from quality control to the 
analysis of variance ; a particularly interesting section describes the technique of sequential 
analysis. ‘The author assumes that his reader has only a very slight knowledge of mathe- 
matics. "This means that most of the formulae have to be given without proof. Instead, 
there is usually a careful explanation of the purpose of the technique under discussion and a 
fully worked out éxample showing its application. 

The treatment throughout is remarkably entertaining. There are a large number of 
diagrams and tables to assist in the solution of the many examples, which are often both 
original and amusing. Sometimes the attempt to be easily understood is carried too far—it 
hardly seems necessary to explain how to square a number—but the reader can easily skip 
those parts which are too elementary for him. The author does not attempt to hide his 
political and other prejudices; this cannot be regarded as a fault since it adds to the general 
liveliness of the text. The treatment is particularly good in stressing both the error 
inherent in any statistical calculation and the many pitfalls which may be encountered in 
understanding the meaning of the results. Another virtue is that the author frequently 
points out that it would be dangerous to use the formulae which he gives without consulting 
a more advanced textbook. Even so, this book is likely to be worth a place in many 
research laboratories. J:2OUBe 


The Challenge of Atomic Energy, by E. H. 8. Buruop. Pp. x+137. (London: 
Lawrence and Wishart, 1951.) 6s. 6d. 3 


This booklet reviews the various aspects of atomic energy from the point of view of : 
the general reader, and it fulfils its aim admirably. The physical introduction is written 
in a very simple and straightforward style, and it is no mean praise to say that the directness 
of the scientific approach is maintained in the later chapters discussing the economic, 
social and political aspects of the problem. 'The professional physicist concerned with 
this problem will find the first part useful if he has to prepare a popular lecture on the 
subject, and from the second part he will learn as much as anybody else. 


L. ROSENFELD, 
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Elements of Thermodynamics and Statistical Mechanics, by E. O. Hercus. 
Pp. ix+153. (Melbourne, Australia: University Press; Cambridge: 
University Press, 1950.) 21s. 


This book is divided into two distinct and separate halves, one on thermodynamics and 
the other on statistical mechanics, representing the lecture courses given to second and third 
year students respectively. 

In the first half a review of the deductions which can be made from thermodynamics 
about systems of different sorts—chemical, radiating, electrochemical or purely mechanical 
—is made, preceded by a short review of the two main laws of thermodynamics themselves. 
The reviewer feels that this development is excellently carried out, but that the basis for the 
laws on which it rests is treated very scantly. In this he differs from the author and, since 
in fact the belief in the truth of the first law is one which grows up with each student during 
his study of mechanics, electricity, chemistry and other subjects, it must be conceded that it 
is not very important whether, before starting on formal thermodynamics, he is or is not 
given a summary of the evidence for a law which by that time is almost axiomatic to him. 

The second part, on statistical mechanics, is very clear indeed, and can be recommended 
without any hesitation. 

For the student there remains only the question whether he can afford the price, since 
he will in any case have to get a book on heat which will itself treat both these subjects. The 
answer depends on the depth of his purse and the amount of time he can devote to reading. 
If both conditions permit, he will do well to buy the book. Tobey 


Heat and Thermodynamics, by Mark W. ZEMANSKY. Pp. xiv+465. 3rd Edition. 
(London: McGraw-Hill, 1951.) 42s. 6d. 


The author of this textbook is Professor of Physics in the City College of New York, 
and, therefore, has practical experience of teaching, and in fact says that the book represents 
an amplification of lecture notes. This is the third edition, which in itself shows that it 
must have been very popular in its country of origin. 

To a reader who was not acquainted with the earlier editions it appears at the outset 
to attack the subject in a formal and rather pedantic way, but this impression disappears 
after the first ten or twelve pages. It has, in fact, quite a realistic and modern outlook, 
especially on the side where the subject approaches the field of physical chemistry. The 
chapter on heat transfer is also written from an extremely practical point of view. 

The chapter on the physics of very low temperatures is quite masterly, and must have 
been re-written since the earlier editions, for its whole treatment of liquid helium IT, 
supraconductivity and related matters reflects the ideas now current in this field, which 
had hardly crystallized when the earlier editions appeared. 

The scope of the book is adequate for a degree course, and will often provide reference 
material for the research worker. It is well written and well produced, and if the price 
seems very high we may ask what commodity does not now seem very dear. 

J. H. A. 


Heat, by R.C. Brown. Pp. ix+270+ viii. (London: Longmans Green, 1952.) 
13s. 6d. 


This is one volume of a series which the author is writing to cover the branches of physics 
at the Intermediate Science level. At this stage the descriptive and experimental approach 
is called for, and the aim must be to give the student a broad familiarity with the terrain. 

Dr. Brown writes well, and has given clear descriptions of the major phenomena usually 
treated under the heading of ‘‘Heat’’ and also of one hypothesis, the kinetic theory of gases. 
He has naturally not taken up the deductive side of the subject exemplified in thermo- 
dynamics. 

His treatment of thermometry is particularly commendable for its accuracy in sorting ° 
out the concepts of heat and temperature, and in explaining why and how we formulate a 
scale of temperatures.. The sections on calorimetry and on radiation are also particularly 
good, and the book can be confidently recommended. It has some examples for the 
student’s practice. J. H. A. 
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Vektor- und Dyadenrechnung fiir Physiker und Techniker, by ERWIN LOnR. 
Pp. xv +488. 2nd Edition. (Berlin: de Gruyter, 1950.) DM. 24. 


Professor Lohr originally intended to carry out a complete revision of his book but, 
owing to the war, he has had to be satisfied with the addition of about 70 pages of 
supplementary material to an otherwise unchanged first edition. Being a physicist, the 
author writes primarily for physicists and engineers, but has provided at the same time 
a most adequate reference book for the applied mathematician working in a variety of fields. 

In the first part a vector is defined geometrically in terms of its cartesian components, 
the idea of scalar products, complex vectors and dyads following in a simple manner. 
The quadric surface associated with a dyadic receives a fair amount of attention, including 
discussion of invariants and the Cayley-Hamilton equation. 

The second part deals with differential operations together with a detailed treatment of 
vector and tensor fields, line, surface and volume integrals, sources, vortices and double 
layers. 

The third part, which covers about half the book, deals with a wide range of physical 
applications, including theoretical mechanics, deformable media, differential geometry of 
curves and surfaces, elasticity, hydrodynamics of perfect and viscous fluids, electro-magnetic 
theory, optics, matrix and quantum mechanics. These are treated in surprisingly great 
detail considering the space at the author’s disposal. 

The additional material in the Supplement includes such topics as potential theory, 
special relativity, four-dimensional formulation of Maxwell’s equations and Dirac’s theory 
of the electron. 

Although the book covers familiar ground and contains nothing substantially new, it 
must nevertheless be regarded as an advanced treatise, and its use to the beginner is further 
limited by the unfamiliar notation employed and by the fact that there are no exercises. 

There is a comprehensive 16-page index. A. N. GORDON. 


Introduction to Modern Algebra and Matrix Theory, by O. SCHREIER and 
E. SPERNER. Pp. vili+378. (New York: Chelsea, 1951.) No price given. 


This book is the slightly shortened English version, in one volume, of the authors’ 
work, originally in German, entitled Introduction to Algebra and Analytical Geometry and 
published in two volumes. It deals with the theory of vectors, determinants and the 
concept of field. It introduces briefly the theory of groups and ends with a long and more 
comprehensive chapter on linear transformations and matrices. 

The authors are writing primarily for students of mathematics, but their demands of 
previous knowledge on the part of their readers is slight. Anyone who has mastered the 
chief operations of algebra and who has been introduced to the basic principles of geometry 
and the calculus will be able to follow the course presented here with ease and, indeed, with 
enjoyment. 

Students of physics will find here a sound and very readable presentation of a notation 
which has in modern times become essential in their study. 

They will find here an explanation of the terms and some propositions which are used, 
and often taken for granted, in textbooks of modern physics. 

This is particularly true of the calculus of matrices and linear transformations which 
have come to be of fundamental importance to students of wave mechanics. 

The book can be recommended with confidence to physics students in the final stages 
of their undergraduate course and in early postgraduate years. H. DB 


Oeuvres scientifiques de Paul Langevin. Pp. 687. (Paris: Centre National de 


la Recherche Scientifique; London agents: H. K. Lewis, 1950.) In paper 
44s.; in cloth 53s. 


The appearance of most of the original papers of Paul Langevin collected together in 
one volume must be regarded as an event of great importance. ‘This collection is not quite 
complete, as many papers which are readily accessible are not included; a list of these papers 
is given. In time this collection covers the entire first half of the twentieth century, for 
his first paper appeared in 1902 and his last in 1949 : in substance it covers a range so wide 
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as to include important papers on ionization of gases, kinetic theory, thermodynamics, 
electromagnetic theory, magnetism, relativity, radioactivity, ultrasonics, nuclear physics and 
gravitation. This is indeed a contribution to physics in the grand manner. An obituary 
notice of Langevin has been published in the Proceedings of the Physical Society (1947, 59, 
1041), where tribute to his immense range of achievement was made. His earlier work, 
coming as it did at a vital period in the development of physics, played a most important 
role in effecting the transition from the classical to the modern point of view, and can be 
regarded as forming a bridge between the statistical mechanics of the late 19th century and 
the early quantum ideas of the 20th century. In particular, Langevin’s work is an excellent 
‘example of the power and scope of classical statistics in dealing with the fundamental 
-collisional phenomena involved in the conduction of electricity through gases. ‘Then, 
again, it is interesting to note that he was a pioneer worker in relativity, and also that Einstein 
first learned of de Broglie’s original notions of wave mechanics from Langevin. For these 
reasons this volume of collected papers is, and probably will continue to be, of profound 
interest to the student of physics, and particularly to the young research worker, since, in 
-addition to important theoretical papers, it also gives numerous papers describing his 
-experiments. 

Langevin’s first papers dealt with the properties and motions of ions (the existence of 
‘free electrons in gases had not yet been demonstrated), and his interest in this subject no 
doubt originated during the time he spent at the Cavendish Laboratory under Sir J. J. 
Thomson. An outstanding problem which excited the interest of physicists at the turn of 
the century was that of the explanation of the electrical conductivity produced by x-rays or 
radioactive emanations. It was the application of classical kinetic theory to the experimental 
determinations of conductivity, and of current decay, which enabled the nature of the ion 
carriers to be found, the processes of ionization and recombination to be investigated, and 
the mechanisms of ionic transport, diffusion and recombination to be understood in terms of 
the collisions of elastic sphere atom models. 

By consideration of the polarization of neighbouring neutral atoms produced by a 
charged ion, and also the resulting attraction effecting reduction of the free path, the conse- 
-quent collection of molecules to form large clusters about a single ion was envisaged. ‘These 
charged molecular groups, whether singly or in clusters, were shown to be the carriers of 
electricity. It was in this particular field that Langevin made his first outstanding contri- 
bution. In the series of papers in 1902-1905 (Ann. Chim. Phys., 5, 245), culminating in the 
famous paper Une formule fondamentale de théorie cinétique, he treated the motion of ions 
in gases, and rigorously derived the well-known formula for the drift velocity W : 


W=0°815(Ze/m)(A/C){(m+M)/M}”, 


of an ion of charge e, free path A, r.m.s. velocity of agitation C(S> W), mass m moving under a 
field Z in a gas of molecular mass MV in which the ion exerts no electric forces except at 
collision (low gas pressure). ‘The ions are considered to be nearly in statistical equilibrium 
with the gas atoms : thus restriction is now known to be unnecessary, since the formula also 
applies to the motions of electrons under a regime of elastic conditions when the electrons 
.are not in statistical equilibrium with the gas atoms, provided, however, that their energy 
distribution is Maxwellian. This latter consideration, again, only affects the numerical 
-coefficient, which only varies from about 2/3 to unity over a wide range of electron energy 
-distribution functions. Replacement of the force Ze by that due to the concentration 
gradient leads to Langevin’s formula for the diffusion coefficient of ions. The application of 
his formulae to the measurements of the mobilities of clusters, single ions, or of electrons led 
to determinations of the mass factor and of the collisional cross sections of the ions : and, 
in fact, in the hands of Townsend and Tizard, such application suggested in 1913 the 
possibility of a variation of atomic cross section in collisions with slow electrons. That 
Langevin had peculiar insight in this subject is seen by the fact that his basic formula held 
the field for about 25 years, in fact until the discovery of exchange phenomena and the 
-calculations of Mott and Massey based on quantum scattering, and was used in the inter- 
pretation of the experimental work of Tyndall and Powell. 
The other great contribution with which Langevin’s name is always associated is in the 
‘theory of magnetism, and this was made in 1905. The theory of electrons in the hands of 
Lorentz was making great strides, but doubt was expressed whether the magnetic theories 
.of Ampére and Weber could be expressed adequately in terms of the new electron theory. 
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It was Langevin who set these doubts at rest. Curie’s work on the temperature dependence 
of paramagnetism was almost completely explained by Langevin on the basis of the classical 
kinetic theory on the assumption that paramagnetic susceptibility was due to the reserve of 
permanent magnetic dipoles in e:ch molecule. ‘The results deduced that over a large range 
of temperature paramagnetic susceptibility is inversely proportional to the absolute 
temperatures—known as the Curie—Langevin law. This work today still forms a most 
valuable introduction to the electron theory of magnetism. The full explanation, however, 
had to await the application of the quantum theory of spin, with its restrictions on the 
orientations of the molecular magnets in a magnetic field. 

This book of 687 pages gives a curriculum vitae of Langevin and a complete bibliography 
of his scientific works, and has altogether thirteen chapters. In spite of its size, this volume 
is quite pleasant to handle, and the print is clear and easy to read. The first chapter 
describes his experimental works on ionization of gases by x-rays, his experimental and 
theoretical work on ionic mobility and recombination, including the polarization of neutral 
atoms produced by the ion charge, and the properties of secondary x-rays ; he contributed 
to these subjects throughout his life. 'The second chapter is devoted to considerations of 
ionization in the atmosphere and the properties of ion clusters. 

The next three chapters are of a theoretical character and give his papers on fundamental 
kinetic theory, electromagnetic inertia, relativity, and also his classic work in the theory 
of magnetism. The sixth chapter presents a collection of his papers on relativity, and this. 
is followed by chapters on radioactivity and units. Then follow some twenty pages. 
devoted to a critical review of classical, quantum and statistical mechanics. Chapters ten 
and eleven are concerned with his contribution to ultrasonics, under-water acoustics and 
divers technical problems ; these represent his scientific work war in World War I. He 
originated and developed the piezoelectric under-water detector, and carried out trials 
in Toulon harbour. The two final chapters are devoted to his most recent work on nuclear 
physics and gravitation, and are of a fundamental nature. 

It is not possible within the scope of a short review to discuss the chapters in detail : it: 
is sufficient to stress the scope and quality and importance of the work described, to admire 
the clarity of its exposition, and very strongly to recommend this collection to all interested. 
in the history and development of an important aspect of modern physics. 

F. LLEWELLYN JONES. 


Quantum Mechanics, by ALFRED LANDE. Pp. x+307. (London: Pitman, 
1951.) 40s. 


In recent years there,have appeared a fair number of good introductory treatises on 
Quantum Mechanics, each taking the student from the beginnings of quantum theory up to- 
the description of relativistic electron theory and field quantization. It is typical of the 
many-sided aspects of quantum mechanics that these treatments all have considerable 
points of difference, not only in the detailed problems studied but also in the fundamental 
style of: pproach and in the order of treatment. 

Professor Landé’s book is an excellent addition to this series. As the author himself 
states in the Preface, the emphasis of his approach is placed on “ the physical background of 
quantum mechanics and its close relation to classical experience, rather than on a rigid 
axiomatic formulation of the mathematical method”. The first three chapters are devoted 
to a study of the equivalence of waves and particles and of the older quantum theory. 
Schrédinger’s equation for matter waves is then developed, and the next two chapters are 
devoted to some of the exact and approximate problems of wave mechanics. In chapter VI 
the transition is made to matrix mechanics based on the wave-mechanical representation. 
There follow treatments of some non-steady state problems, and of electronic spin leading 
to the helium atom and hydrogen molecule. There is also a brief discussion of atomic and 
molecular spectra‘and of quantum statistics. The last three chapters are concerned with 
Dirac’s electron theory, and with the quantization of the electromagnetic and meson fields. 

It will be seen that the ground covered in the book is very wide. The discussion is, 
however, generally clear, though at times rather condensed so that in places the mathe- 


matics is left to carry the argument. Th» inclusion of relativistic electron theory and field 
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quantization in an elementary book poses severe problems, which cannot be said to have 
been completely solved by the author. However, it was no doubt right to include these 
topics, and the treatment offered at least forms an introduction to more advanced treatises. 
In general, the book forms a reliable and well presented account of wave mechanics 
suitable for Honours Physics students of reasonable mathematical attainments, and for 
young research workers. An interesting feature is the short summaries at the end of each 
chapter, which preserve the main line of argument against the mathematical details. 
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The Growth from the Melt of Single Crystals of Low Melting Point Metals, by 
A. J. Goss and S. WEINTROUB. 


ABSTRACT. The results of a series of experiments on the preparation of tin, lead, zinc, 
and bismuth single crystals in the form of rods are given. The effects on the crystallization 
of various thermal factors, in particular the temperature gradient and the rate of growth, are 
studied and the results for the four metals are compared. 


Asymmetry of Texture in Cubic Crystals Grown from the Melt, by P. R. ROWLAND. 


ABSTRACT. Careful electro-polishing of single crystal spheres of copper (face-centred. 
cubic system) grown from the melt has revealed that they appear to consist of a number of 
parallel crystallites of square cross section with their lengths parallel to the [001] axis which 
lies nearest in direction to the ‘ direction of growth’ in the melt. 


The Growth of Thallium-Bromo-Iodide Crystals, by K. M. GREENLAND. 


ABSTRACT. Crystals of thallium-bromo-iodide, or KRS-5, are transparent to radiation 
of wavelengths between 0-5 and about 40%. This note briefly describes the growing 
of KRS-5 single crystals and their optical working. The crystals are grown in vertical 
cylindrical crucibles by a progressive freezing of the melt from bottom to top. The physical, 
chemical and mechanical properties of the material are summarized. 


Growth of Crystals of Caesium Bromide, and the Infra-Red Transmission Limit in 
Optical Materials, by A. C. G. MENZIES. 


ABSTRACT. Following the development of methods for growing single crystals of 
materials having good optical properties, the search has been intensified for new substances. 
One such is caesium bromide, and this has been grown by the Kyropoulos—Pohl method 
of drawing from the melt; a modification of the standard process is described. 

In the search for new materials, it is useful to have some way of predicting the useful 
long-wavelength transmission limit, and a rule is stated which forms a useful guide. 


Optics of the Focusing X-Ray Microscope, by J. Dyson. 


ABSTRACT. A brief description is given of Kirkpatrick and Baez’ x-ray reflecting micro— 
scope. The types of aberration possible in optical systems without axial symmetry are 
discussed and the effect of spherical aberration on the diffraction pattern of an object point 
is evaluated. The resolving power of a circular mirror is calculated, and a method of 
correcting field obliquity is described. A method is given for the design of a figured mirror 
to give smaller spherical aberration with zero field obliquity. 

Design equations are given for a system employing two circular mirrors in each plane, by 
which spherical aberration, obliquity and curvature of field canbe corrected. The computed. 
performance of the system is described. 


The Focal Properties and Spherical Aberration Constants of Aperture Electron: 
Lenses, by M. M. MacNauGuton, 


ABSTRACT. The focal properties and spherical aberration constants of a number of 
aperture electron lenses of varying diaphragm thickness and varying strength were deter-- 
mined numerically using potential distributions obtained (a) from a rigorous solution due: 
to Ollendorff, and (6) on a resistor network analogue. The focal lengths of ‘weak’ aperture. 
lenses of finite diaphragm thickness agree closely with those given by Davisson and Calbick’s 
formula for the ‘infinitely’ thin aperture lens; ‘ strong’ aperture lenses of all thicknesses. 
show systematic deviations from this formula. 'The primary spherical aberration constants. 
of diverging aperture lenses are small, and of such sign that they cannot be used for the- 
correction of the spherical aberration of converging electron lenses. 
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. Investigation of Magnetic Lenses having the Axial Field H(0, z)=y/z", by 
U. F. GIANOLA. 


ABSTRACT. Examination of the requirements for reducing the spherical aberration of an 
electron lens by modifying the field shows that the hypothetical lens field H(0, 2)=y/2 
should have small aberrations. 

It is shown that the properties of a certain class of this general field configuration can 
be described in terms of an exact analysis. In particular, the electron optical properties of 
the y/z? field are described in detail. 

The y/2? field is found to have small aberrations, but practical conditions limit appreciable 
improvement in resolving power. The characteristics of the y/z° field are also examined and 
are found to give little further advantage. 


Analysis of Polarized Light Reflected from Absorbing Materials at Normal Incidence, 
by J. Wooprow, B. W. Morr and H. R. Harnes, 


ABSTRACT. Expressions are derived for the intensity of the analysed light reflected from 
absorbing crystals (a) under crossed polarizing units and (6) with the analyser slightly 
rotated from the crossed position. Four intensity maxima per complete rotation of the 
specimen are obtained when the polarizers are crossed, but under certain circumstances, as 
the analyser is rotated from the crossed position, the number of maxima may be reduced to 
two. ‘These relationships are illustrated by experimental curves on specimens of tin, 
uranium, bismuth and stibnite. 


The Magneto-Resistance of Ferromagnetic Al-Si-Fe Alloys, by R. PARKER. 


ABSTRACT. By extending the treatment in an earlier paper an equation is derived for the 
saturation magneto-resistance of some mixed ferromagnetic alloys as a function of composi- 
tion and temperature. Results calculated from this equation are in good agreement with 
experiment for the Al-Si-Fe alloy system. 


Sparking Potentials in a Transverse Magnetic Field, by J. M. SOMERVILLE. 


ABSTRACT. A theory of sparking between parallel plates in a transverse magnetic field, 
put forward by Valle in 1950, is discussed and shown to be in qualitative disagreement with 
observation. By making allowance for the distribution of electron free paths, and for electron 
recapture by the cathode, a modified theory is produced which is in better agreement with 
observation. Further refinement of the theory requires an approximate expression for the 
first Townsend coefficient which is valid when a magnetic field is present. ‘The possibility 
of obtaining such an expression is discussed. 


Metal Transfer during Static Loading and Impacting, by E. RABINOWICZ. 


ABSTRACT. Experiments are described which investigate the metal transfer which 
takes place when a radioactive hemispherically ended slider is pressed normally into a flat 
surface of the same or another metal. 

It is found that metal transfer in the form of a number of small fragments is always 
observed, and this supports the view that metallic junctions are formed between metal 
surfaces in contact. Experiments in which the load is varied suggest that at the higher loads 
the oxide layer on the metal surfaces is broken up to a larger extent than at lower loads, and a 
more than proportional increase in metallic interaction and transfer takes place. An 
analogous effect is observed in the presence of boundary lubricants. 

When the surfaces are impacted together, very similar results are obtained. Somewhat 
less pick-up is observed than for static loading, and the difference is probably due to the fact 
that it takes time for strong junctions to be formed. Impact experiments with surfaces 
covered by lubricants show that a lubricant layer may be trapped between the surfaces, and 
this produces a large reduction in pick-up without greatly reducing the amount of plastic 


deformation. 
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Further Experiments on Evaporation from Small, Saturated, Plane Areas into a 
Turbulent Boundary Layer, by D. R. Davies and T. 5. WALTERs. 


ABSTRACT. ‘This paper describes an experimental investigation designed to evaluate 


ee ee 


the effect of width and length on rates of evaporation from small, saturated, plane, rectangular — 


areas placed in a turbulent boundary layer. An approach is made to the problem of esti- 


mating the distribution of temperature on the evaporating surface by using a number of — 


thermojunctions situated at various points under the surface. The higher rate of evapora- 
tion near the edges of the area is well illustrated by this method, greater temperature 
differences between the surface and air stream being observed near the edges due to.greater 
latent heat effects. Provided the crosswind width of the evaporating area is not less than 
about half the downwind length, Sutton’s two-dimensional solution of this limited area 
problem is found to lead to tolerable accuracy in the circumstances tested. 


Experiments on the Mechanics of Rubber—III: Small Torsions of Stretched Prisms, 
by A. N. Gent and R. S. RIVLIN. 


ABSTRACT. A theorem, relating the torsional modulus for small torsions of a rod of : 


incompressible highly elastic material, held in a state of simple extension, to the extension 
ratio and force necessary to maintain the extension, is verified experimentally for rods of 
vulcanized rubber of rectangular cross section. 
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(6) Rb,(Cu+Zn)(SO,)2 . 6H,O 
Fig. 2. Variation in width of hyperfine lines with temperature. 
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Fig. 3. Absorption lines in KyCu(SO,),.6H,O with magnetic field at 45° to 6 axis. 
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